For  Reference 


NOT  TO  BE  TAKEN  FROM  THIS  ROOM 


©X  IIBB1S 

aramsraw 

MBHOTBUW 


The  University  of  Alberta 
Printing  Department 
Edmonton,  Alberta 


UNIVERSITY  OF  ALBERTA 
LIBRARY 

Regulations  Regarding  Theses  and  Dissertations 


Typescript  copies  of  theses  and  dissertations  for  Master’s  and  Doctor’s 
degrees  deposited  in  the  University  of  Alberta  Library,  as  the  official  Copy  of 
the  Faculty  of  Graduate  Studies,  may  be  consulted  in  the  Reference  Reading  Room 

only. 


A  second  copy  is  on  deposit  in  the  Department  under  vhose  supervision  the 
vork  vas  done.  Some  Departments  are  willing  to  loan  their  copy  to  libraries, 
through  the  inter-library  loan  service  of  the  University  of  Alberta  Library. 

These  theses  and  dissertations  are  to  be  used  only  with  due  regard  to  the 
rights  of  the  author.  Written  permission  of  the  author  and  of  the  Department 
must  be  obtained  through  the  University  of  Alberta  Library  when  extended  passages 
are  copied.  When  permission  has  been  granted,  acknowledgement  must  appear  in  the 
published  work. 

This  thesis  or  dissertation  has  been  used  in  accordance  with  the  above 
regulations  by  the  persons  listed  below.  The  borrowing  'library  is  obligated  to 
secure  the  signature  of  each  user. 


Please  sign  below: 


Date 


Signature 


Institution 


THE  UNIVERSITY  OF  ALBERTA 


MASS  TRANSFER  AT  BUBBLE  FORMATION 


by 


HARISH  KALRA , 


B .  Sc  . 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  FULFILLMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 
OF  MASTER  OF  SCIENCE  IN  CHEMICAL  ENGINEERING 


DEPARTMENT  OF  CHEMICAL  AND  PETROLEUM  ENGINEERING 


EDMONTON,  ALBERTA 
FALL,  1969 


Digitized  by  the  Internet  Archive 
in  2019  with  funding  from 
University  of  Alberta  Libraries 


https://archive.org/details/Kalra1969 


UNIVERSITY  OF  ALBERTA 


FACULTY  OF  GRADUATE  STUDIES 


The  undersigned  certify  that  they  have  read,  and 
recommend  to  the  Faculty  of  Graduate  Studies  for  acceptance, 
a  thesis  entitled  "Mass  Transfer  at  Bubble  Formation" 
submitted  by  Harish  Kalra,  B.Sc.,  in  partial  fulfilment  of 
the  requirements  for  the  degree  of  Master  of  Science  in 
Chemical  Engineering. 


>0nsctq9003  *10*1  e^JEjjfaBTtO  10  V-'  UvS\-  .'nsrninooivr 

- 


I 


ABSTRACT 


High  speed  cine  photography  was  used  to  study  mass 
transfer  at  bubble  formation  for  the  CO^-water  system.  The 
amount  of  mass  transfer  was  determined  by  the  difference  of 
the  bubble  volume  at  formation  as  determined  from  photographs 
and  the  bubble  volume  calculated  from  flow  rate  and  fre¬ 
quency  measurements. 

Point  data  of  the  bubble  periphery  was  taken  from 
photographic  films  and  was  processed  using  computer  pro¬ 
grams  to  determine  the  bubble  surface  area  and  volume.  Two 
programs  were  written  for  the  purpose;  one  of  them  used  a 
numerical  approach  and  the  other  was  based  on  an  empirical 
equation  used  to  describe  the  bubble  shape.  Variation  of 
bubble  surface  and  volume  with  time  could  be  studied  using 
the  numerical  program. 

Bubble  surface  area  was  found  to  vary  according  to 
nl 

the  equation  A=k^t  only  up  to  the  point  where  bubble 
necking  and  oscillations  start.  n^  was  found  to  be  of  the 
order  of  0.850.  Sharp  increases  in  area  were  observed 
during  necking.  The  volume  of  the  bubble  was  found  to 
vary  in  a  substantially  linear  manner  with  time. 

Mass  transfer  coefficients  for  bubble  formation  were 
determined  for  the  CO^-water  system  in  a  column  having  a 
diameter  of  20.2  cm.,  for  a  liquid  velocity  of  1.0  cms/sec., 
bubble  frequencies  of  0.5  to  10.0  bubbles/sec.,  an  orifice 
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internal  diameter  of  0.408  cms. ,  gas  flow  rates  of  0  to  1.5 

cc/sec j  and  at  a  liquid  temperature  of  25.0°c.  The  values 

obtained  were  in  the  range  of  98  to  222  cms/hr  for  the 

frequencies  studied.  The  order  of  magnitude  agrees  well 

with  results  of  previous  investigations  for  mass  transfer 

during  bubble  rise  [3,30,50,52,69,70,71,72],  but  the  value 

of  the  mass  transfer  coefficient  during  bubble  formation  was 

found  to  be  greater  than  that  measured  during  bubble  rise. 

This  difference  was  thought  to  be  due  to  a  liquid  velocity 

superimposed  on  bubble  growth.  The  total  amount  of  mass 

transfer  was  found  to  depend  more  closely  on  Tf0'^^, 

0  5 

where  Tf  is  the  total  bubble  age,  rather  than  Tf  as 
suggested  by  some  previous  investigators  during  studies  of 
mass  transfer  from  drops  in  liquid-liquid  systems  [5,42 
53,54,55,56,57]. 
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CHAPTER  1 


INTRODUCTION 

Gas  liquid  contactors  are  indispensible  in  industry 
and  because  of  their  extensive  application  it  is  desirable 
to  be  able  to  design  such  equipment  reliably.  To  do  this 
it  is  useful  to  have  as  much  information  as  possible  on  the 
fundamental  aspects  of  gas  liquid  contacting.  This  study, 
namely  the  study  of  mass  transfer  during  bubble  formation, 
was  undertaken  with  this  viewpoint  in  mind. 

The  dispersion  of  a  gas  phase  in  the  form  of  bubble 
in  a  liquid  is  one  of  the  modes  of  gas-liquid  contact.  The 
study  of  the  mechanisms  of  bubble  formation  and  of  mass 
transfer  to  or  from  such  bubbles  has  been  in  progress  for 
the  last  thirty  years.  Bubbling  phenomena  has  been  cat¬ 
egorically  subdivided  into  three  distinct  regions  namely: 

(i)  Bubble  formation 

(ii)  Bubble  rise 

(iii)  Bubble  breakup  at  the  free  liquid  surface. 

Mass  transfer  will  occur  during  each  of  the  above  mentioned 
regions . 

Most  of  the  work  done  to  date  in  the  study  of  mass 
transfer  has  been  restricted  to  the  bubble  rise  period. 

Mass  transfer  during- bubble  formation  has  been  treated  as 
an  end  effect  and  any  quantitative  information  obtained  for 
the  amount  of  mass  transfer  occurring- during  these  regions 
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has  therefore  been  indirect.  Calderbank  [1],  Haselden  et. 
al.  [2],  Bowman  et .  al .  [3],  Licht  et.  al .  [4,5]  and  Dixon 
et .  al .  [6,7]  conducted  studies  of  mass  transfer  during 
bubble  rise.  In  order  to  separate  the  amount  of  mass 
transfer  occurring  during  this  period  it  was  necessary 
that  the  amount  of  mass  transfer  occurring  during  bubble 
formation  and  during  bubble  breakup  be  either  known  or 
be  eliminated  by  a  modification  of  the  experimental  tech¬ 
nique.  The  end  effect  due  to  bubble  breakup  or  bubble 
coalescence  could  be  eliminated  through  equipment  modifi¬ 
cations  but  an  estimate  of  the  amount  of  mass  transfer  during 
bubble  formation  was  obtained  by  an  extrapolation  of  a 
plot  of  total  mass  transfer  versus  liquid  height  to  zero 
liquid  height.  Their  conclusions  indicated  that  a  con¬ 
siderable  amount,  varying  from  25  to  50%,  of  the  total  mass 
transfer  occurred  during  bubble  formation  [1 , 2 ,4 ,5 ,6 , 7] • 

From  this  it  was  concluded  that  the  rate  of  mass  transfer 
during  bubble  formation  is  probably  enhanced.  No  attempt 
was  made  to  measure  the  actual  rate  of  mass  transfer 
during  bubble  formation,  even  though  theoretical  postulations 
for  an  enhanced  rate  of  transfer  during  this  region  were 
put  forward  [1,2]. 

Subsequently,  however,  Calderbank  et .  al.  [8] 
reported  a  direct  study  of  mass  transfer  during  bubble 
formation  for  Carbon  dioxide.  Sulphur  dioxide  and  Acetylene 
bubbles  formed  in  a  stagnant  water  phase  at  low  bubble 
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frequencies.  High  speed  cine  photography  was  used  to 
determine  the  change  of  bubble  surface  area  with  time  and 
the  mass  transfer  coefficient  was  represented  by  a  Higbie 
[9]  type  of  equation.  It  was  reported  that  no  enhanced 
rate  of  mass  transfer  occurred  during  bubble  formation 
but  no  quantitative  data  were  given  to  support  this 
conclusion . 

The  presence  of  a  considerable  amount  of  uncertainty 
as  to  the  amount  of  mass  transfer  occurring  in  the  region 
of  bubble  formation  justified  a  further  study  of  this 
phenomena  so  that  concrete  evidence  may  be  obtained  one 
way  or  the  other.  An  experimental  program  similar  in 
approach  to  Calderbank  et .  al.  [8]  was  drawn  up  wherein 
the  range  of  investigation  was  extended  and  the  effect  of 
liquid  velocity  was  also  Imposed  on  the  system.  High 
speed  cine  photography  was  used  to  photograph  the  forming 
bubbles  and  computer  programs  were  written  to  evaluate  the 


bubble  area  and  volume  as  a  function  of  time. 
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LITERATURE  REVIEW 

An  extensive  study  of  dispersions  has  been  done 
with  both  gas-liquid  and  liquid— liquid  systems.  A  dis¬ 
persed  liquid  phase  is  referred  to  as  drops  whereas  a  gas 
is  said  to  form  bubbles.  The  general  mechanisms  of  drop 
and  bubble  behaviour  are  very  similar  in  a  number  of  ways 
and  thus  an  understanding  of  one  throws  considerable  light 
on  the  behaviour  of  the  other. 

The  formation  of  drops  or  bubbles  has  been  found 
to  be  a  function  of  the  physical  properties  of  the  system 
involved,  gas  or  liquid  velocity,  and  the  dimensions  of 
the  apparatus  such  as  the  orifice  diameter. 

Mass  transfer  during  bubble  rise  has  been  studied 
in  considerable  detail  and  the  effect,  on  mass  transfer, 
of  variables  governing  the  formation  of  bubbles  has  been 
established.  The  study  of  mass  transfer  during  bubble 
formation  has  been  somewhat  limited  in  comparison. 

The  literature  review  on  this  subject  has  been 
subdivided  into  the  following  categories: 

(1)  Bubble  behaviour 

(2)  Mass  transfer  during  bubble  (or  drop)  rise  and 
during  bubble  (or  drop)  formation. 

(3)  Effect  of  Surface  Activity  on  Bubble  behaviour 
and  Mass  Transfer. 
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Bubble  Behaviour. 


10 


Effect  of  gas  flow  rate  on  bubble  formation. 

The  phenomena  of  bubble  formation  may  be  categorized 
into  the  following  regions  based  on  increasing  gas  flow 
rate  : 

(a)  Constant  Volume  region 

(b)  Transition  region 

(c)  Constant  frequency  region 

(d)  Turbulent  region. 

As  the  gas  flow  rate  is  increased  the  bubble  be¬ 
haviour  corresponds  to  one  of  the  above  mentioned  regions 
in  the  order  in  which  they  have  been  mentioned.  The 
change  from  one  region  to  the  other  has  been  characterized 
in  terms  of  the  orifice  Reynolds  number  which  is  defined 


by  : 


Re  -  pdUg« 


(2.1) 


where  =  density  of  the  dispersed  phase 

Ug  =  velocity  of  the  gas  at  the  orifice 
6  =  Orifice  diameter 

U  =  Viscosity  of  the  continuous  phase. 

This  criteria  is  suitable  for  analysis  of  bubble  behaviour 
in  industrial  size  equipment  only.  Any  values  of  Re 
mentioned  in  the  subsequent  discussion  bear  this  restriction 
In  addition,  some  overlap  of  Reynolds  numbers  between 
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regions  is  often  found. 

(a )  Constant  Volume  Region 

At  low  gas  flow  rates  corresponding  to  a  value  of 
Re  upto  100  or  150,  the  gas  bubble  forms  discretely  and 
slowly  at  the  given  orifice  and  is  substantially  constant 
in  volume.  The  volume  of  the  bubble  at  detachment  from 
the  orifice  may  be  evaluated  by  a  force  balance  between 
the  surface  tension  and  the  buoyant  forces.  The  frequency 
of  bubble  formation  is  proportional  to  the  gas  flow  rate. 
These  observations  have  been  supported  by  Van  Krevelin  et . 
al.  [10],  Quigley,  Johnson  and  Harris  [11],  Hughes  et .  al 
[12],  Coppock  and  Meiklejohn  [13],  Leibson  et .  al .  [14]  and 
Datta,  Napier  and  Newitt  [15]. 

Davidson  et .  al.  [16],  however,  have  reported  that 
the  bubble  volume  increased  slightly  with  increase  of  gas 
flow  rate  in  this  region  and  thus  that  the  bubble  freq¬ 
uency  is  slightly  less  than  proportional.  Benzing  and  Myers 
[17]  concluded  from  their  investigations  that  the  bubble 
volume  had  a  tendency  to  increase  with  bubble  frequency. 

In  effect  this  substantiates  the  statement  made  by  Davidson 
et .  al .  [ 16 ] . 

(b )  Transition  Region 

The  transition  region  is  the  changeover  from  the 
constant  volume  to  the  constant  frequency  region.  A  range 
of  Re  of  100  -  300  has  been  suggested  for  this  region  by 
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Valentin  [18],  Neither  the  bubble  volume  nor  the  bubble 
frequency  are  constant  in  this  region  of  bubble  behavior. 

The  bubble  formation  though  discrete  is  characterized  by  a 
tendency  of  the  bubbles  to  proceed  up  the  column  in  pairs. 

A  bubble  forms  at  the  orifice,  breaks  from  it  and  appears 
to  hover  a  few  centimeters  above  the  orifice  till  it  is 
joined  by  the  next  bubble.  Both  of  them  then  move  up  the 
column  together.  This  phenomenon  persists  for  short  inter¬ 
vals  of  time.  This  implies  that  the  bubble  formation  often 
reverts  to  the  formation  of  single  discrete  bubbles  as 
observed  in  the  constant  volume  region.  This  observation 
has  been  supported  by  Davidson  et .  al.  [16],  Benzing  and 
Myers  [17]  and  Davidson  and  Schuler  [19]  have  reported 
an  increasing  predominance  of  the  inertial  effects  during 
this  region  of  bubble  formation.  Benzing  and  Myers  [17] 
refer  to  it  as  a  kinetic  impingement.  This  effect,  which 
is  due  to  a  higher  velocity  of  gas  flowing  into  the  bubble, 
tends  to  dislodge  the  bubble  at  smaller  volumes,  but  this 
effect  is  more  than  compensated  by  a  large  amount  of  gas 
flowing  into  the  bubble  at  breakoff.  As  a  result  of  this 
the  bubble  volume  increases  with  increase  of  frequency. 

( c  )  Constant  Frequency  Region 

With  a  further  increase  in  gas  flow  rate  the  bubble 
frequency  becomes  independent  of  gas  flow  rate  and  approaches 
a  constant  maximum  value,  which  is  a  function  of  the  orifice 
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diameter.  The  bubble  volume,  however,  increases  proport¬ 
ionally  with  gas  flow  rate.  The  orifice  Reynolds  numbers 
characterizing  this  region  have  been  suggested  as  1000  to 
2100 . 

Due  to  a  high  gas  flow  rate,  the  gas  bubbles  tend 
to  merge  into  one  another  resulting  in  coalescence  after 
detachment  from  the  orifice.  The  size  of  the  bubbles  at 
detachment,  for  a  given  gas  flow  rate  in  this  region,  is 
reproducible,  indicating  that  the  bubble  formation  is 
descrete.  Though  evidence  to  support  the  above  statements 
has  been  provided  by  Davidson  et .  al .  [ 1 6  ]  ,  Benzing  and 
Myers  [17],  Eversole  et .  al.  [20],  and  Calderbank  [21], 
Leibson  et .  al .  [14]  have  indicated  that  some  small  bubbles 
are  formed  in  this  region  by  the  rupture  and  collapse  of 
the  elongated  bubble  neck.  Thus  some  small  bubbles 
accompany  larger  ones  and  their  fraction  gradually  in¬ 
creases  with  increasing  gas  flow  rate. 

(d )  Turbulent  Region 

As  the  gas  flow  rate  is  increased  further  from  Re 
of  about  2100  to  40,000,  bubble  coalescence  begins  to  take 
place  closer  and  closer  to  the  orifice;  it  becomes  difficult 
to  see  discretised  bubbles  forming  at  the  orifice.  A 
variety  of  sizes  of  bubbles  is  formed.  Gas  bubbles  break 
up  on  formation  into  smaller  bubbles  and  the  size  dis¬ 
tribution  of  this  break  up  differs  in  different  regions  of 
the  fluid.  There  are  larger  bubbles  near  the  orifice 
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and  much  smaller  ones  higher  up  the  column.  Bubbles 
during  this  region  have  been  termed  as  being  spherical 
capped  or  toroidal  in  shape  [ 1 8 ] . 

A  spherical  bubble  has  a  relatively  stable  configur¬ 
ation  while  a  toroidal  bubble  breaks  up  on  formation.  At 
higher  rates  of  gas  flow  the  proportion  of  toroidal  bubbles 
increases.  Leibson  et .  al .  [14],  Seimes  [22]  and  Calder- 
bank  [21]  have  reported  agreement  to  the  above  mentioned 
facts.  Leibson  et .  al.  [14]  have  classified  the  bubble 
formation  above  a  Reynolds  number  of  10,000  where  turbulence 
is  considered  fully  developed,  as  being  similar  to  a 
disintegrat ing  liquid  jet  in  air.  The  gas  ensuing  from  the 
orifice  does  so  in  the  form  of  a  jet  and  disintegrates  into 
coarse  and  fine  bubbles  near  the  office.  Rennie  and  Evans 
[23],  however,  indicate  that  even  at  Reynolds  numbers  as 
high  as  40,000,  there  is  still  bubble  formation  followed 
by  bubble  breakup,  the  bubbles  formed  being  all  toroidal. 

Effect  of  Chamber  Volume 

The  chamber  volume  is  defined  as  the  volume  upstream 

of  the  orifice  to  the  nearest  restriction  in  the  gas  line. 

Variation  of  the  chamber  volume  serves  as  a  very  effective 

method  to  control  bubble  stability  and  results  in  a  change 

of  bubble  volume.  In  order  to  study  the  effect  of  the 

chamber  volume  on  bubble  behaviour  it  is  convenient  to 

define  a  capacitance  number  N  given  as  follows : 

c 
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N 


c 


4g<Pc-Pd)Vc 

2  2 

n  6  pdc 


(2.2) 


where 

p  =  density  of  continuous  phase 
c 

Pd  =  density  of  dispersed  phase 

V  =  chamber  volume 
c 

6  =  orifice  diameter 

c  =  acoustic  velocity  in  the  gas. 

The  formation  of  stable  descrete  bubbles  occurs 
when  Nc  is  below  a  certain  critical  value.  For  large 
chamber  volumes  where  is  greater  than  critical  the 
bubble  formation  becomes  unstable  and  a  phenomena  referred 
to  as  chain  bubbling  occurs.  During  chain  bubbling  a  large 
number  of  gas  bubbles  form  in  succession  followed  by  a 
pause  which  is  subsequently  broken  by  another  chain  of 
bubbles.  This  process  repeats  itself.  When  the  chamber 
volume  is  small  the  formation  of  one  gas  bubble  is 
enough  to  decrease  the  gas  pressure  in  the  chamber  up¬ 
stream  of  the  orifice  sufficiently  so  as  to  prevent  the 
next  bubble  from  forming  till  such  time  that  the  gas  flow 
restores  the  chamber  pressure.  The  formation  of  a  chain 
of  bubbles  at  values  of  Nc  greater  than  critical  is  due 
to  the  fact  that  one  bubble  does  not  decrease  the  chamber 
pressure  sufficiently  and  as  a  result  a  chain  of  bubbles 
is  formed  till  the  chamber  pressure  drops.  The  pause 
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following  a  chain  of  bubbles  is  used  to  restore  the  chamber 


pressure  to  the  pressure  of  the  system.  The  formation  of 
bubbles  at  very  large  chamber  volumes  is  said  to  occur  at 
constant  pressure  conditions. 

When  the  chamber  volume  is  very  much  smaller  than 
the  critical  value,  the  volume  of  the  bubbles  formed  is 
independant  of  the  chamber  volume.  On  increasing  the 
chamber  volume,  the  bubble  volume  increases  till  the  value 
of  Nc  approaches  the  critical.  At  the  critical  value  very 
sharp  changes  in  bubble  size  take  place.  If  the  bubbling 
phenomena  favours  the  formation  of  multiple  bubbles  the 
bubble  volume  drops  sharply. 

The  critical  chamber  volume  is  itself  a  function  of 
gas  flow  rate  for  a  given  orifice  size.  At  high  gas  flow 
rates  the  critical  value  of  Nc  drops  sharply.  This  implies 
that  at  high  gas  flow  rates  small  changes  of  chamber  volume 
could  upset  the  stability  of  bubble  formation.  Hughes  et . 
al.  [12],  Davidson  et .  al .  [ 1 6 ] ,  Bowman  et .  al.  [3],  Hayes 
et .  al.  [24],  Leibson  et .  al .  [14]  and  Mahoney  et .  al.  [25] 
have  conducted  investigations  of  the  effect  of  chamber 
volume  on  bubble  behaviour  and  their  conclusions  are  in 
agreement  with  the  above  given  facts.  Bowman  et .  al .  [3] 
proposed  a  dimensionless  correlation  to  predict  the  effect 
of  this  variable  on  bubble  behaviour. 
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Most  workers  agree  on  the  effect  of  orifice  dia¬ 
meter  on  bubble  formation.  During  the  constant  volume 
region  the  bubble  size  is  proportional  to  the  cube  root 
of  the  orifice  diameter.  This  is  readily  seen  by  making 
a  force  balance  between  the  surface  tension  and  the 
buoyancy  forces.,  at  bubble  breakoff.  As  the  orifice  size 
is  increased  the  size  of  the  constant  volume  bubble  will 
eventually  become  smaller  than  the  orifice  diameter.  The 
gas-liquid  interface  will  then  become  horizontal  at  the 
orifice  perimeter  and  this  causes  unstable  bubble  formation. 
Liquid  weeping  through  the  orifice  is  sometimes  caused  by 
the  use  of  too  large  an  orifice.  Davidson  et .  al .  [16] 
suggested  that  stable  bubbles  in  the  constant  volume  region 
cannot  be  formed  at  an  orifice  of  radius  0.79  cms .  but 
found  that  such  bubbles  could  be  formed  at  orifices  of 
radius  0.64  cms.  They  then  concluded  for  the  systems 
they  investigated,  that  orifices  of  radius  larger  than  0.7 
cms.  cannot  produce  stable  bubbles  in  the  constant  volume 
region . 

At  high  gas  flow  rates,  that  is  during  the  constant 
frequency  region,  Van  Krevelin  et .  al.  [10]  have  reported 
that  the  bubble  volume  is  independent  of  orifice  diameter. 
The  maximum  frequency  of  bubble  formation  varies  inversely 
with  the  orifice  radius.  Davidson  et .  al .  [ 1 6 ]  have  pro¬ 
posed  correlations  relating  the  maximum  bubble  frequency 
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with  the  orifice  radius.  They  suggested  that  the  maximum 
frequency  varies  inversely  as  the  cube  root  of  the  orifice 
radius.  Other  correlations  relating  bubble  volume  with  gas 
flow  rate  and  orifice  diameter  were  also  suggested  by  them 
and  a  data  fit  within  1%  for  high  gas  flow  rates  and  within 
5%  for  the  transition  region  has  been  claimed. 

Effect  of  Liquid  Seal  Height 

The  effect  of  liquid  seal  height  has  been  studied 
by  several  authors  but  common  observations  for  all  of  them 
do  not  exist . 

Quigley  et .  al.  [11],  Davidson  et .  al.  [16],  Coppock 
and  Meiklejohn  [13],  Hayes  et .  al .  [24]  and  Towell  et .  al . 
[26],  have  all  reported  a  negligible  effect  of  seal  height 
on  bubble  formation.  Bowman  et .  al.  [3]  have  shown  that 
this  is  not  so.  A  range  of  seal  heights  from  2.5^  to  30.5 
cms .  was  investigated  and  their  conclusions  indicate  that 
bubble  frequency  increases  with  increase  of  seal  height, 
particularly  for  seal  heights  greater  than  20  cms.  (CO^  - 
water  system) .  The  bubble  size  therefore  decreases  with 
increase  of  seal  height  for  the  same  gas  flow  conditions. 
They  developed  an  exponential  equation  in  an  attempt  to 
describe  the  submergence  effect  quantitatively.  They 
suspected  that  the  submergence  effect  may  diminish  with 
increasing  orifice  throat  length.  Leonard  et .  al.  [27] 
also  investigated  this  variable  and  their  conclusions 
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reiterate  the  above  facts. 


Effect  of  Orifice  Throat  Length 

Hughes  et .  al.  [12],  Leibson  et .  al .  [14]  and  Haye 
et .  al.  [24]  have  reported  that  there  is  no  effect  of 
varying  orifice  throat  length  on  bubble  behaviour.  Benzing 
and  Myers  [17]  disagreed  with  this  and  reported  that  for 
an  orifice  of  size  0.438  cms .  as  the  throat  length  was 
decreased  from  4.6  cms.  to  0.1  cms.,  no  significant 
change  in  bubble  volume  occurred  till  the  throat  length 
reached  a  value  of  0.33  cms.  when  the  bubble  volume 
dropped  sharply.  Bowman  and  Johnson  [3]  have  reported  that 
where  chamber  volume  cannot  be  sufficiently  reduced  to  get 
a  descrete  stable  bubble,  reduction  of  throat  length  helps 
to  obtain  stable  bubble  formation. 

Effect  of  Column  walls  and  wettability 

The  effect  of  the  column  wall  on  bubble  formation 
becomes  evident  only  when  the  column  diameter  becomes 
comparable  with  the  bubble  diameter.  The  gas  bubble  then 
gets  elongated  as  in  gas  lifts.  If  the  ratio  of  column 
to  orifice  diameter  is  10  or  greater,  the  wall  effect  is 
negligible . 

Orifice  surface  wettability  is  important  for  pre¬ 
venting  liquid  leakages  through  the  orifice.  Stable  bubble 
formation  requires  that  the  orifice  surface  be  non-wettable 
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by  the  liquid.  Teflon  covers  have  been  successfully  used 
to  enhance  non-wettability  where  orifice  diameters  were  too 
large  and  liquid  leakages  occurred  [3]. 

2 . 2  Mass  Transfer  during  Bubble  Rise  and  Bubble  Formation 

Bubble  Rise 

Mass  transfer  from  or  to  bubbles  or  drops  rising 
singly  in  a  continuous  phase  is  a  function  of  the  fluid 
properties,  the  bubble  size  and  the  flow  conditions  of  the 
dispersed  phase.  The  analysis  of  the  rate  of  mass  transfer 
in  the  continuous  phase  may  be  categorized  as  follows: 

1)  Mass  transfer  to  or  from  immobile  surfaces 

2)  Mass  transfer  to  or  from  mobile  surfaces 

Each  of  the  above  mentioned  categories  may  be  further  sub¬ 
divided  according  to  bubble  size  as  follows : 


(a) 

Small 

spherical 

bubbles 

(  d  <0.2 

e 

cms .  ) 

(b) 

Fluid 

oblate  spheroids 

(0.2<d  <1.8 

cms .  ) 

(c ) 

Fluid 

spherical 

caps 

(  de,1.8 

cms .  ) 

Mass  transfer  from  each  of  the  above  mentioned  bubble  sizes 
can  be  analysed  according  to  the  relative  importance  of 
viscous  and  inertial  forces  in  a  given  fluid  system.  That 
is,  the  analysis  differs  with  the  magnitude  of  the  Reynolds 
number . 

The  rates  of  mass  transfer  from  small  drops  or 
bubbles  in  a  continuous  phase  can  be  analysed  rigorously 
because  the  bubbles  are  spherical  in  shape  and  rise  or  fall 
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in  straight  paths.  Larger  drops  or  bubbles  are  generally 
irregular  both  in  shape  and  path.  An  exact  theoretical 
treatment  to  determine  the  mass  transfer  rates  from  such 
bubbles  is  therefore  very  difficult.  It  is  customary  to 
assume  the  larger  bubbles  to  have  well  defined  shapes  and 
to  rise  in  straight  lines.  In  addition  bubble  shapes  are 
assumed  to  be  axisymmetric 3  that  is,  symmetric  about  the 
axis  in  the  direction  of  rise. 

By  assuming  that  the  Peclet  and  Schmidt  numbers  are 
large  Calderbank  [28]  has  derived  an  equation  in  dimensionless 
coordinates  relating  the  tangential  velocity  of  the  bubble 
surface  with  the  concentration  of  the  solute  in  the  liquid 
phase.  The  restriction  of  large  Peclet  and  Schmidt  numbers 
is  based  on  the  fact  that  the  concentration  boundary  layer 
is  assumed  to  be  far  thinner  than  the  velocity  boundary 
layer.  The  equation  derived  is  given  as 
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The  significance  of  the  coordinates  u3v,x,y,r  for  an 
axisymmetric  body  of  revolution  is  shown  in  Figure  1. 

The  boundary  conditions  applicable  to  equation  2.3  are 

C1  =  ls  y-L  =  0  all  x1 

C1  =  °5  yl  =  °°  a11  X1 

Using  these  boundary  conditions  equation  2.3  can 
be  solved  for  the  mobile  and  immobile  surfaces  mentioned 
above . 

Case  1:  Immobile  Interface 

This  implies  that  u1Q=0.  Using  this  fact  equation 
2.3  has  been  solved  by  Locheil  [29]  under  the  above  given 
boundary  conditions  to  give  the  following  general  mass 
transfer  equation. 

Sh  =  0.641  (i,nre)[/1(u'  r  )1/2r  dx1]2/3Pe1/3 - (2.4) 

T  0  10  1  11 


Where  A  is  the  total  surface  area  of  the  body. 
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Figure  1:  Coordinates  For  an  Axisymmetric 
Body  of  Revolution 
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Case  2 :  Mobile  Interface 

For  a  moving  fluid  boundary  u]_o>>uio^ls  and  equation 
2.3  may  be  solved  under  the  same  boundary  conditions  as 
before  to  give  the  following  equation  [29]. 

Sh  =  (|)V2(!]ji)[/1Uior2aXi]l/2peV2 - (2.5) 

These  equations  can  be  easily  reduced  to  more  familiar  forms 
for  special  conditions  of  flow  and  bubble  size. 

(a)  Solid-Spheres  (Immobile  Interface) 

Depending  on  the  system  under  investigation,  bubbles 
or  drops  below  a  certain  size  behave  essentially  like  rigid 
spheres.  For  purposes  of  mass  transfer  analysis  the  inter¬ 
face  may  thus  be  considered  immobile.  Hammerton  and  Garner 
[30]  have  indicated  that  bubbles  formed  in  water  behave  as 
rigid  spheres  below  a  size  of  0.02  cms .  and  when  formed  in 
glycerol  this  limit  goes  up  to  0.2  cms. 

For  a  sphere  r^  =  sin  0,  A  =  4llr  and  dx^  =  d0. 
Substituting  in  equation  2.4,  the  equation  for  mass  transfer 
becomes : 

Sh  =  0.64l[/(u'  sin  0)1//2sin  0  d©]2//^Pe - (2.6) 

0  10 

This  equation  has  been  derived  by  Baird  and  Hamielec  [31] 
also.  For  cases  where  the  inertial  terms  are  very  small. 
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that  is  Re<l,  Friedlander  [32]  obtained  the  following 
simplification  of  equation  2.6. 
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(2.7) 


Equation  2.7  has  been  verified  for  dissolution  of  small 
bubbles  by  Calderbank  et .  al.  [33],  for  liquid  drops  by 
Ward  et .  al .  [3^]  and  for  solid  shperes  by  Aksel’  Rud  [35]. 

For  cases  where  the  inertial  forces  cannot  be 
neglected,  that  is  Re>l,  a  boundary  layer  type  of  equation 
can  be  obtained  from  the  Navier-St okes  equation  and  using 
a  quartic  profile  for  the  velocity  distribution  in  the 
boundary  layer,  assuming  no  separation  at  the  sphere’s 
surface,  the  following  equation  can  be  obtained. 


Sh  =  0.62  Re1/2  Sc1/3 


(2.8) 


Separation  of  the  boundary  layer,  however,  occurs  at 
0-108°  [36]  and  the  equation  for  mass  transfer  occurring 


over  the  front  surface  of  the  sphere  upto  the  point  of 
separation  is  given  as 


Sh  =  0.56  Re1/2  Sc1/3 


(2.9) 


This  equation  was  also  obtained  by  Froessling  [37]  and 
verified  for  dissolution  of  small  drops  by  Calderbank  et.al. 


' 

o  :  '  : ’!  ‘ 

'■  w  no.'  B  :  ■■  1  ■' 


[33]  and  Hammerton  and  Garner  [30]. 


(b )  Fluid-Spheres  (Mobile  Interface) 

If  the  values  of  r^  A  and  dx^  for  a  sphere  be 
substituted  in  equation  2.5,  the  equation  for  a  fluid 


sphere  becomes 


0 


l/2pel/2 


(2 .10) 


Once  again  depending  on  whether  Re<l  or  >1,  specific  sub¬ 
stitutions  for  u^q  may  be  made  in  equation  2  .10  to  obtain 
particular  solutions  for  mass  transfer  from  fluid  spheres. 

For  Re<l  Hadamard  [38]  and  Rybczynsky  [39]  developed 
velocity  profile  about  a  fluid  sphere.  Using  their  relations 
for  the  tangential  surface  velocity  the  equation  for 

mass  transfer  becomes 

Sh  =  0.65  ( — ^3- — )1/2pe1/2 - (2.11) 

Vud 

Where  y  and  y  ,  are  the  viscosities  of  the  continuous  and 
dispersed  phases  respectively.  Part  of  the  work  of  Ham¬ 
merton  and  Garner  [30]  and  Ward  et .  al.  [34]  correlates  to 
this  equation  reasonably  well. 

For  cases  where  Re>l  Chao  [40]  has  proposed  equations 
for  the  velocity  distribution  both  inside  and  outside  the 
sphere  as  dimensionless  perturbations  from  the  potential 
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flow  fields.  Using  his  equations  for  the  velocity  profiles 
Locheil  [29]  has  derived  the  following  equation  for  mobile 
spheres . 


Sh  =  1.13  [1- 


2+(3Md/yc) 


;  .1/2  0  1/2 

l+(p,y,/p  y  )  Re 

'a  d  c 


1.45  jl/2  pel/2_ 


(2.12) 


Griffith  [ 4 1 ]  and  Herrtjes  et .  al .  [42]  have  used  this 
equation  to  correlate  data  for  water  drops  in  isobutanol 
and  ethyl  acetate  drops  in  water. 

For  gas  bubbles  in  a  liquid  phase ^equation  2.12  can 
be  reduced  to 


Sh  =  1.13  [l-2,^p]1/2Pe1/2 - (2.13) 

Re  7 

This  equation  applies  for  Re  in  the  range  100<Re<400.  For 
very  high  Re  i.e.  potential  flow  equation  2.13  reduces  to 
the  very  familiar  Higbie  [9]  equation. 

Sh  =  1.13  Pe1/2— - (2.14) 


This  equation  was  first  derived  by  Stefan  [43]  and  later  by 
Boussinesq  [44]  and  Higbie  [9].  The  rates  of  dissolution 
of  large  bubbles  has  been  satisfactorily  correlated  with 
this  equation  by  Hammerton  and  Garner  [30],  Higbie  [9],  West 
[45]  and  Bowman  and  Johnson  [3]. 

Ruckenstein  [46]  has  also  developed  equations  for 
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fluid  spheres  both  in  the  range  Re>>l  and  Re<<l  by  using 
an  interpolation  method.  This  allows  an  evaluation  of 
mass  transfer  in  the  intermediate  range  also. 
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Fluid  Oblate  Spheroids 

When  dealing  with  irregular  but  axisymmetric  shapes 
such  as  oblate  spheroids,  it  is  difficult  to  specify  precise 
relations  for  surface  area.  Use  is  made,  therefore,  of  a 
shape  factor  called  eccentricity  which  is  defined  as  the 
ratio  of  the  bubble  width  to  bubble  height.  Oblate  spheriods 
are  in  general  ellipsoidal  shapes  and  the  eccentricity  is 
given  by  the  ratio  of  the  major  to  minor  axes. 

The  equivalent  diameter  of  these  bubbles  lies  in 
the  range  0.2<d  <1.8  cms .  and  the  Reynolds  numbers  are  very 
high.  The  flow  is  therefore,  almost  always  assumed  as 
potential.  Zahm  [47]  has  developed  velocity  profiles  for 
flow  around  oblate  spheroids  and  using  his  relations  in 
combination  with  equation  2.5  the  following  equation  for 
mass  transfer  can  be  obtained. 


Sh  =  {|(l+K*)}1/2{ 


2 . 26e  1/3(e2-l)1/2 
c  c 


•}Pe1/2 - (2.15) 


,  2  ..  ,1/2X1  r  ,,  2  ,  x  1/2  -i 

ec(ec_1)  +ln[ec+(ec-l)  ] 


Where  K*  is  a  complex  function  of  the  eccentricity 
e  .  This  equation  has  been  successfully  used  to  interpret 
the  data  of  Hammerton  and  Garner  [30],  Higbie  [9]  and  West 

[45]. 


' 

4 
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Fluid  Spherical  Caps 
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When  the  bubble  diameter  exceeds  about  1.8  cms .  the 
bubbles  assume  mushroom  like  shapes  that  are  termed  as 
spherical  caps.  Potential  flow  is  assumed  to  describe  the 
velocity  profiles  of  the  fluid  around  such  bubbles.  The 
interfacial  velocity  is  given  by  the  potential  flow  value 
for  a  sphere  namely 


u1Q  =  3/2  sin  0 - (2.16) 

Using  equation  2.16  and  neglecting  mass  transfer  from  the 
rear  of  the  bubble,  Locheil  [29]  obtained  an  equation 
given  as 


(3eUl))2/3 

Sh  =  1.79  [-,£ — ] 

e2+4 

c 


(2.17) 


Rosenberg  [48]  and  Takadi  and  Maeda  [49]  have  suggested  a 
value  for  e  =  3.5  for  all  spherical  capped  bubbles.  Using 
this  value  equation  2.17  reduces  to 

Sh  =  1.28  Pe1/2 - (2.18) 


which  is  very  nearly  the  same  as  the  Higbie  equation  2.14. 

Most  mass  transfer  studies  from  bubbles  involve  an 
analysis  of  the  variation  of  the  mass  transfer  coefficient 
with  equivalent  bubble  diameter.  This  is  generally  pre¬ 
sented  in  the  form  of  a  plot.  If  the  range  of  bubble 
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diameters  covered  is  wide,  the  mass  transfer  coefficient 
conforms  to  the  Proessling  [37]  type  of  equation  for  rigid 
spheres  i.e.  small  bubble  diameters  and  then  increases  as 
the  bubble  diameter  increases.  For  larger  bubbles  the 
Higbie  type  of  equation  is  more  closely  followed.  The 
curve  for  mass  transfer  coefficient  vs.  equivalent  bubble 
diameter  as  determined  by  experiment  thus  lies  in  between 
the  values  predicted  by  the  equations  of  Froessling  [37] 
and  Higbie  [9].  The  experimental  curve  passes  through  a 
maximum  at  a  diameter  which  is  a  function  of  the  properties 
of  the  system  under  study.  These  type  of  curves  have  been 
obtained  by  Hammerton  and  Garner  [30],  Leonard  and  Houghton 
[27],  Zieminski  et .  al .  [50]  and  Siddique  [51]. 

The  mass  transfer  studies  so  far  reviewed  have  all 
been  done  in  stagnant  liquid  systems.  Himmelblau  et .  al . 

[52]  have  reported  that  the  mass  transfer  coefficient  from 
gas  bubbles  in  a  moving  liquid  is  higher  than  that  in  the 
stagnant  system.  During  a  study  of  multiple  dispersions, 
however,  Calderbank  [28]  has  stated  that  no  such  effect  of 
liquid  velocity  is  observed. 

Bubble  Formation 

In  contrast  to  the  bubble  rise  period,  the  study  of 
mass  transfer  during  bubble  formation  has  been  rather  limited. 
As  mentioned  earlier,  indications  of  a  significant  amount 
of  mass  transfer  occurring  during  bubble  formation,  were 
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obtained  by  extrapolating  the  bubble  rise  data  to  zero  liquid 
height.  This  was  done  after  the  coalescence  end  effect  had 
been  eliminated.  To  investigate  these  findings  in  greater 
detail,  direct  studies  of  mass  transfer  during  drop  or 
bubble  formation  have  been  undertaken  by  Licht  and  Pansing 
[5] 5  Heertjes  et.  al.  [42],  Groothius  and  Kramer  [53], 

Ilkovic  [54,  55],  Michels  [56],  Popovich  et .  al .  [57]  and 
Calderbank  and  Patra  [8]. 


Five  different  models  for  mass  transfer  during  drop 


or  bubble  formation  have  been  reported  in  the  literature 
[5 , 42 , 53 ,  ( 54 , 55 ) , 56 ] .  All  of  them  have  been  based  on 
the  unsteady-state  diffusion  of  the  solute  in  the  cont¬ 
inuous  phase.  Different  assumptions  describing  the  hydro¬ 
dynamics  of  the  system  were  superimposed  to  obtain  explicit 
equations  for  the  mass  transfer  occurring  up  to  a  time  ’ tT. 
Popovich  et.  al.  [57]  assumed  that  the  variation  of  drop 
surface  area  with  time  could  be  suitably  described  by  the 
equation 


n 


1 


(2.19) 


and  reduced  all  the  final  equations  developed  in  the  five 
models  mentioned  above  to  two  common  equations  of  the  form 


=  const (k1) (C|-C^) 


D  [(2n,-l)/2] 

H  Z  1 


(2.20) 


. 
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Where  is  the  mass  transfer  rate  at  time  ft'.  The 
equation  for  total  mass  transfer  in  time  ft’  is  given  by 
integration  of  equation  2.20  and  is  of  the  form 

s  b  /d~  [(2n  +1)/ ,] 

Na  =  const  (kj^HC^-C^)/ j  t  1  2 - (2.21) 

The  constant  in  equations  2.20  and  2.21  is  a  function  of  the 

exponent  n^  and  assumes  a  different  value  for  each  model 

proposed.  In  all  cases  the  value  of  n1  has  been  assumed 

to  be  close  to  2/3  which  is  the  value  for  a  dispersed 

spherical  particle  growing  uniformly  with  time.  For  a 

2  -2/3 

spherical  drop  or  bubble  =  TId^Tf  and  equations  2  .20 
and  2  .21  take  the  form 

=  const  (k  )(C®-C^)/  DII  d2  T"2/3  t1/6 - (2.22) 

and 

Na  =  const  (Ca-Ca)/"dI  d2  T”2/3  t7/6 - (2.23) 

To  obtain  the  mass  transfer  during  drop  formation  to  the 
point  of  release  substitution  of  t  =  Tf  reduces  the  above 
equations  to 


const 


(°X) 


m-l/2 

f 


(2.24) 
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and 

Na  =  const  (CA-CA)/“Dlf  d2  T^/2 - (2-25) 

Popovich  et .  al.  [57]  found  that  their  data  cor¬ 
responded  most  closely  to  the  model  proposed  by  Ilkovic  [54, 
55],  wherein  the  value  of  the  constants  for  equations  2.20 
and  2.21  were  predicted  to  be  1.53  and  1*31  respectively. 

The  study  of  Popovich  et .  al.  [57]  was  confined  to 
liquid  drops  of  approximately  0.2  to  0.4  mm.  in  diameter. 

The  assumption  of  spherical  drop  growth  under  these  con¬ 
ditions  is  reasonable.  However,  bubbles  forming  at  an 
orifice  are  generally  less  nearly  spherical  than  liquid 
drops,  especially  if  the  bubble  size  is  of  the  order  of 
0.5  to  1.0  cms .  Under  these  circumstances  the  analogy  of 
liquid  drop  behaviour  to  gas  bubbles  and  the  applicability 
of  the  above  given  theory  is  uncertain. 

Calderbank  and  Patra  [8]  studied  the  mass  transfer 
during  bubble  formation  for  carbon  dioxide,  acetylene  and 
sulphur  dioxide  bubbles  in  water.  They  proposed  an  empir¬ 
ical  relation  for  the  variation  of  bubble  surface  area  with 
time  which  is  given  as 


A  _  ,  ,Tf  Un2 
Am  "  (  Tf  } 


(2  .26) 


where  (Tf-t)  is  the  life  expectancy  of  the  element 
of  surface  generated  at  time  t.  For  low  bubble  frequencies 
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ranging  from  0.2  to  2.0  sec”^.  the  value  of  the  exponent 
n^  was  reported  as  1.86  +0.06.  Assuming  that  each  fresh 
element  of  surface  formed  does  not  mix  with  the  older  ones, 
and  using  the  theory  of  diffusion  into  a  semi-infinite  body, 
a  Higbie  type  of  equation  was  proposed  [9].  This  equation 
is  valid  for  short  exposure  times  only  and  is  given  as 


d  N .  ,  _ _ 

xrrrr  =  (C®-C„)/d  - (2.2 

dtdA  a  a  ✓  Trirpyi 

Substituting  for  dA  from  equation  2.26  and  integrating  be¬ 
tween  time  limits  0  and  the  following  equation  may  be 
obtained . 


4n, 


N 


(C®-Cfl)  f TT  a  t 
A  A  /  m 


A  2n2+l  V~A  "A  .  n 


1/2 


m  f 


(2. 28) 


With  the  value  of  n^  =  1.86  the  above  expression  becomes 


NA  ■  i-56  <c]tca>  /§“  AmTf1/2' 


(2.28a) 


Assuming  that  the  quantity  of  material  transferred  as  given 
by  equation  2.28a  is  the  same  as  would  be  dissolved  if  the 
bubble  were  formed  instantaneously  and  remained  at  its 
maximum  size  for  a  time  Tf,  Calderbank  et .  al .  [8]  suggested 
the  following  equation  for  the  mass  transfer  coefficient. 


t 
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kT  = 
Lav 


(2.29) 


They  have  reported  that  a  good  agreement  was  obtained  with 
equation  2.28  for  the  carbon-dioxide-water  and  acetylene- 
water  systems.  No  data  to  this  effect  has  been  reported. 

As  mentioned  earlier,  mass  transfer  during  bubble  rise  may 
be  given  by  equation  2.14  another  form  of  which  is  given  as 


(2.30) 


Where  t  is  the  time  taken  by  the  bubble  to  rise  a  distance 
equal  to  its  diameter  at  its  vertical  rise  velocity.  Based 
on  this  fact,  Calderbank  and  Patra  [8]  have  suggested  that 
the  mass  transfer  rate  during  bubble  formation  is  not  en¬ 
hanced  but  is  lower  than  that  during  bubble  rise. 


2. 3  Effect  of  Surface  Activity  on  Bubble  Behaviour  and 
Mass  Transfer 

A  substance  whose  addition,  even  in  extremely  small 
quantities,  changes  the  surface  properties  of  the  system 
under  study,  is  termed  as  being  surface  active.  The  effect 
of  surface  active  agents  on  bubble  behaviour  and  on  mass 
transfer  from  bubbles  has  been  and  still  is  a  subject  of 
extensive  controversy,  even  though  work  on  this  variable 
has  been  in  progress  for  some  time. 
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There  is  more  agreement  amongst  authors  on  the 
effect  of  surface  active  agents  on  bubble  behaviour  than 
there  is  on  its  effect  on  mass  transfer. 

Addition  of  surface  active  agents  to  the  system 
being  investigated  reduces  the  surface  tension.  As  a 
result  of  this,  gas  bubbles  formed  at  a  given  orifice  for 
identical  conditions  of  gas  flow  rate  are  smaller  in  size 
than  a  system  where  no  surface  active  agents  have  been 
added.  Reduced  surface  tension  and  reduced  bubble  sizes 
have  been  observed  by  Kintner  et .  al .  [58],  Leibson  et .  al. 
[14],  Towell  et .  al.  [26],  Coppock  and  Meiklejohn  [13], 
Benzing  and  Myers  [17],  Eversole  et .  al .  [20],  Calderbank 
et .  al.  [33]  and  Zieminski  et .  al .  [59].  However,  reduction 
of  bubble  size  is  not  indefinite  but  is  a  function  of 
concentration  of  the  surfactant  in  the  system.  Coppock 
and  Meiklejohn  [13]  and  Benzing  and  Myers  [17],  have  reported 
that  at  low  concentrations  of  the  surfactant,  the  bubble 
size  drops  sharply  ,but  as  the  concentration  of  the  surface 
active  agent  is  continuously  increased  beyond  a  certain 
limit jthe  bubble  size  remains  unaffected.  Benzing  and 
Myers  [17]  attributed  this  to  the  fact  that  at  higher 
concentrations  the  bubble  surface  is  rapidly  oriented  to 
its  new  state,  that  is  to  say  that  the  static  surface 
tension  value  is  reached  in  a  very  short  period  of  time. 
Zieminski  et .  al .  [59]  have  also  reported  a  similar 
observation,  but  in  addition  they  state  that  the  structural 
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character  of  the  surfactant  is  also  a  function  in  deciding 
how  much  a  bubble  size  may  be  reduced.  The  long  carbon 
chain  compounds  produce  smaller  bubbles,  and  hence  a  larger 
interfacial  area  for  a  given  gas  holdup. 

The  rate  of  coalescence  of  gas  bubbles,  after 
formation  at  the  orifice,  is  also  significantly  reduced  by 
the  addition  of  surfactants.  The  use  of  surfactants  is, 
therefore,  sometimes  advised  for  systems  where  coalescence 
becomes  a  serious  problem.  Experimental  data  on  coalescence 
reduction  has  been  presented  by  Zieminski  et .  al .  [59]  and 
Towell  et .  al.  [26]. 

The  mass  transfer  coefficient  in  a  gas-liquid 
system  is  a  function  of  bubble  size  and  is  reported  to  in¬ 
crease  with  increase  of  bubble  size  [27,33].  Since  surface 
active  agents  reduce  bubble  size  Calderbank  et.al.  [33] 
concluded  that  the  mass  transfer  coefficient  drops  with 
addition  of  surface  active  agents.  This  phenomena  has  been 
supported  by  Hammerton  and  Garner  [30],  Baird  and  Davidson 
[60]  and  Leonard  and  Houghton  [27].  On  formation,  a  gas 
bubble  has  some  surface  disturbances  referred  to  as  rip¬ 
pling;  this  is  considerably  reduced  by  the  action  of 
surfactants  in  the  system  and  the  surface  stagnation  thus 
created  reduces  the  mass  transfer  coefficient  [60]. 

In  contradiction  to  the  above  observations  Quigley 
et .  al.  [11]  have  reported  that  no  effect  on  mass 
transfer  was  observed  by  the  addition  of  the  surfactant 
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Alkanoxto  the  system  but  their  data  indicates  that  the 
surfactant  was  added  in  fairly  large  concentrations ,  which 
may  have  been  the  reason  for  their  conclusions.  However, 
Houghton  et .  al .  [6l]  have  reported  an  increase  in  mass 
transfer  coefficient  in  a  system  containing  0.05%  byweight 
of  Lubrol  W. 

The  above  facts  indicate  that  both  the  concentrat¬ 
ion  of  the  surfactant  and  its  chemical  nature  govern  the 
effect  that  its  addition  would  have  on  the  mass  transfer 
coefficient.  This  was  indeed  found  to  be  so.  Zieminski 
et.al.  [59]  have  reported  that  for  a  given  surfactant, 
the  mass  transfer  coefficient  increases  with  concentration 
of  the  surfactant  in  the  low  concentration  range  and  then 
drops  at  higher  concentration  finally  levelling  off  at  a 
fairly  constant  value  at  very  high  concentrations.  The 
maxima  of  the  mass  transfer  coefficient  was  found  to  occur 
at  lower  concentrations  of  the  surfactant  for  higher  carbon 
chain  compounds. 

Further  work  in  this  field  is  necessary  to  really 
establish  the  role  played  by  surface  active  agents  in 
dispersed  systems. 
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THEORY 

3 . 1  Evaluation  of  Mass  Transfer  Coefficient 

The  liquid  phase  mass  transfer  coefficient  can  be 
defined  by  the  following  equation. 


^A 

dt 


kLA(t)«j*-cJ) 


(3.1) 


Where  N.  Is  the  number  of  moles  of  solute  transferred, 
s  b 

C^  and  C^  are  the  concentrations  of  the  solute  in  the 
liquid  phase  at  the  interface  and  in  the  bulk  respectively. 
The  area  available  for  mass  transfer  is  given  by  A  which 
is  a  function  of  the  bubble  age  t. 

The  equation  may  be  rearranged  to  give 

dNA  =  kL(C®-CA)A(t)dt  - (3.2) 


and  integrating  between  the  appropriate  limits 


T 


na  - 


kLav.(CA-CA>  /  dt 


(3.3) 


Where  ^Lav  is  the  averaged  mass  transfer  coefficient  during 
the  time  of  formation  Tf  of  one  bubble.  The  derivation  of 


equation  3-3  assumes  that  the  solute  concentrat ions  at  the 


' 


If 
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interface  and  in  the  bulk  are  constant  over  the  time  Tf. 

If  the  resistance  to  mass  transfer  at  the  interface  is 
negligible  then  represents  the  solubility  of  the  solute 
in  the  liquid  phase. 

Integration  of  equation  3.3  requires  the  value  of 

the  interfacial  area  A  as  a  function  of  time.  To  determine 

it  is  necessary  that  a  suitable  analytical  procedure  be 

available.  Then  if  N.  is  known.  kT  can  be  evaluated. 

A  a  Lav 


3 . 2  Evaluation  of  Bubble  Area  and  Volume 

A  record  of  bubble  growth  can  be  obtained  by  taking 
motion  pictures  of  the  bubble  during  formation.  The 
evaluation  of  bubble  area  and  volume  from  photographs 
requires  that  the  bubble  periphery  be  described  either  by 
a  continuous  function  or  by  discrete  data  points.  Then  if 
the  bubbles  are  assumed  symmetrical  about  the  axial  direction, 
the  surface  and  volume  can  be  determined  by  the  following 
integrals : 


b 


1 


s  =  2n  /  f ( v ) 

ai 


/  '2 

/l+f(v)  dv 


(3.4) 


and 


2 

f  (v) 


dv 


(3.5) 


where  f(v)  is  the  function  describing  the  periphery  and 
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f'(v)  is  the  derivative  of  f(v)  with  respect  to  v.  Thus, 
two  approaches  may  be  used  to  describe  the  bubble  periphery. 

Continuous  Function  Approach  -  Computer  Program  I 

Figure  2  represents  a  cardioid  whose  equation  in  the 
parametric  form  is 

r  =  a(l+Cos  0)  - (3.6) 

If  the  portion  of  the  curve  below  AB  is  sliced  off  the  re¬ 
maining  part  of  the  curve  takes  the  shape  of  a  gas  bubble 
in  its  early  stages  of  formation  and  immediately  after 
detachment  from  the  orifice.  A  continuous  function  which 
will  describe  the  periphery  of  such  a  flat  bottomed  bubble 
can  therefore  be  obtained  by  modifying  the  equation  for  a 
cardioid.  Three  modifications  were  attempted: 


(i)  r  =  a(k+Cos  0)  - (3.7) 

(ii)  r  =  a(k+Cosn  0)  - (3.8) 

(iii)  vn  =  a(k+Cos  0)  - (3.9) 


Equation  3.9  was  found  not  to  converge  to  a  definite 
value  of  ’  n'  and  was  therefore  discarded.  The  method  used 
for  testing  whether  or  not  equations  3.7  or  3.8  could  be 
used  to  describe  the  bubble  periphery  was  as  follows : 

Considering  equation  3*7 

0=O=a (k+1 ) 


r  =  a(k+Cos  0) | 


(3.10) 


v  agnriw  K6UU9  6  eJneaaiqai  £  !  I 
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v 


Figure  2 : 


The  Cardioid:  r=a(l+cos0) 


' 


- 

i 


r  =  a ( k+Cos  0 )  I 
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(3.11) 


e=n/2  ak 

These  limiting  values  can  be  measured  from  a  photographic 
film  i.e.  r  @  0=0  and  Jl/^  is  known.  Therefore  'a*  and 
’k’  are  known  by  simultaneous  solution  of  equation  3«10 
and  3.11. 

Substituting  the  values  of  ’a'  and  1 k T  in  equation 
3.7,  'rT  is  obtained  as  an  explicit  function  of  0  only. 

By  varying  0  between  0  and  n/2  at  small  intervals,  a  plot 
of  Trf  vs.  0  can  be  generated.  If  the  equation  3*7  is  a 
good  fit  of  the  bubble  periphery,  the  plot  of  fr’  vs.  0 
should  match  it . 

Considering  equation  3*8,  the  limiting  values  of 
'r'  @  0=0  and  0=11/ 2  are  once  again  given  by  equations  3.10 
and  3.11.  Rearranging  equation  3.8 

,  n~ 

r-ak  =  a  cos  0 

or  r-ak  =  cosn  0 

a 

taking  logs  of  both  sides 

log  ( /— a/ ) =  n  log  Cos  0 
a 

or  log  (~k)  =  n  log  Cos  0 - (3.12) 

a 

’a’  and  T k T  are  known  through  equations  3.10  and  3.11.  A 

data  table  of  log  (— k)  vs.  log(Cos  0)  at  values  of  O<0<n/9 

a 
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can  therefore  be  constructed  using  a  photographic  film. 

The  value  of  the  exponent  ’n’  can  be  determined  by  a  least 

square  fit  through  the  linear  plot  of  log( — k)  and  log(Cos  0). 

a 

TrT  is  then  known  explicity  as  a  function  of  0  only.  The 
generated  data  of  Tr’  vs  0  may  be  plotted  to  see  if  the 
fit  obtained  is  satisfactory.  Figure  3  shows  both  equations 
3^  7  and  3.8  plotted  in  comparison  to  the  actual  bubble 
periphery  obtained  from  a  film.  The  match  shows  clearly 
that  equation  3*8  gives  a  satisfactory  fit  and  can  be 
used  to  develop  the  surface  area  and  volume  of  flat-bottomed 
bubbles . 

From  Figure  2 


v  =  r  Cos  0 - (3.13) 

and  f(v)  =  r  sin  0 - ( 3 •  1 4 ) 

f’(v)  =  d  f ( v )  =  d  f ( v )  .  d0 - (3*15) 


dv  d0  dv 


Using  equations  3#13,  3.14  and  3*15  with  equation  3*8  and 
substituting  in  equations  3/4  and  3/5,  the  following  eq¬ 
uations  may  be  obtained. 


S-2nf-a2Sin0(k+Cosn0)[Sin20{k+(n+l)Cosn0}2+ 

n/2 


{nSin20Cosn  ^0-kCos0-Cosn+20}2 ]2//2d0 - (3*16) 


' 

o  $  id  pwori  i  C  :■?  'gf'L 

br  '  ■  i.i  fo  ;  qoleveo  o) 


.  •  C  ,  ■<■: 

•  * 


o  Actual  bubble  periphery 
a  Equation  r  =  a(k+cos0) 


v  □  Equation  r=a(k+cosno) 


Figure  3:  Plot  of  Equations  3-7  and  3.8  in  comparison 
to  actual  Bubble  Periphery 
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and 


V=n/-a3Sin30(k+Cosn0)2{k+(n+l)Cosn0}d0 


(3.17) 


The  derivation  of  equations  3.16  and  3*17  is  given 


in  Appendix  A.  The  complete  analysis  using  equation  3.8 
was  programmed  for  the  IBM  360/67  computer  and  the  program 
was  tested  on  a  flat-bottomed  sphere  of  radius  3-03  cms . 
The  surface  area  as  evaluated  by  equation  3*16  was  within 
0.5%  and  the  volume  within  1.0%.  This  includes  the  errors 
inherant  in  the  procedure  used  to  obtain  data  from  the 
photographic  film.  The  program  documentation  is  included 
in  Appendix  B. 

Numerical  Approach-Computer  Program  II 


The  bubble  periphery  can  also  be  described  by 


descrete  points.  The  surface  and  volume  can  then  be  eval¬ 
uated  using  numerical  techniques.  The  various  steps  used 
in  generating  the  surface  and  volume  follow  the  sequence 
given  below: 


(i)  Obtain  raw  data  of  (f(v),v)  from  bubble  periphery 


on  film. 


(ii)  Smooth  the  raw  data  using  Forsythe  Orthogonal 
Polynomials  [62]. 


(iii)  Find  the  derivative  f  (v)  at  the  smoothed  points 
using  Lagrangian  Interpolation  [63]. 


(iv)  Obtain  g(v) 


■ 

. 
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Use  Lagrangian  Interpolation  to  generate  equispaced 
data  for  g(v)  [64]. 

Integrate  using  Newton  Cotes  fourth  order  formula 
with  m  subintervals  to  obtain  surface  area  [65]. 

Use  Lagrangian  Interpolation  to  generate  equispaced 
data  for  f (v)  [64] . 

Integrate  using  method  (vi)  above  to  obtain  bubble 
volume . 

The  equations  used  for  the  Forsythe  Polynomials, 
Lagrangian  Interpolation  and  Newton  Cotes  Integration  have 
been  included  in  Appendix  A. 

This  approach  can  be  used  to  evaluate  the  surface 
and  volume  for  any  symmetrical  shape  and  thus  can  be  used 
for  bubbles  during  the  entire  period  of  formation  and 
during  bubble  rise  provided  the  symmetry  is  not  destroyed. 

Program  II  was  used  for  the  analysis  of  data  obtained 
during  this  investigation.  The  program  was  tested  with  the 
same  test  data  as  for  program  I  and  the  errors  are  within 
0.1/S  for  surface  and  1.0 %  for  volume.  The  reliance  on  the 
surface  area  and  volume  is  however  limited  to  1  to  2% 
because  the  errors  involved  in  obtaining  data  from  the 
film  are  of  this  magnitude. 

3 . 3  Evaluation  of  the  Sample  Size  Required  For  Bubble 
Analysis  -  Computer  Program  III 


(v) 

(vi) 

(vii) 

( viii ) 


The  procedure  described  here  can  be  used  to  determine 
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the  sample  size  of  bubbles  that  is  representative  of  a 
given  normally  distributed  population.  Such  an  analysis 
is  necessary  so  as  to  determine  the  number  of  bubbles  that 
must  be  analyzed  in  order  to  obtain  reliable  values  for 
bubble  volume  and  surface  area. 

If  the  size  of  the  given  population  is  and  it 
consists  of  individual  readings  x^; 

Then  the  required  sample  size  N  is  given  by 

s 


N 


s 


x  y 

- ] 

-x) 


2 


(3.18) 


Where  y  is  the  standard  deviation  from  the  mean  and  is 
given  by 


(3.19) 


x  is  the  mean  of  the  population  and  is  given  by 


x 


N 

IP  X. 
i=l  ■ 
N 


(3.  20) 


is  a  factor  which  indicates  the  confidence  that  an  event 

will  fall  within  the  sample  size  chosen. 

Depending  on  the  accuracy  required  the  term  (x  -x) 

s 

may  be  made  as  small  as  desirable.  The  smaller  this  term 
is  the  more  stringent  are  conditions  of  tolerance  and  hence 
the  larger  the  sample  size  required.  For  this  investigation 
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equation  3.18  has  been  used  in  the  following  form 
On  x  y  9 

N  =  [  — - r]  - (3.21) 

0.01  x 

A  program  to  carry  out  the  above  analysis  was  written  and 
is  given  in  Appendix  B. 

3 . 4  Liquid  Velocity  Profile 

Liquid  velocity  profiles  can  be  used  to  characterize 
the  hydrodynamics  of  the  system  when  mass  transfer  during 
bubble  formation  is  measured  under  conditions  of  a  super¬ 
imposed  liquid  velocity.  The  principle  of  polarography  as 
applied  to  electrochemical  reactions  can  be  used  to  measure 
velocity  profiles.  This  method  is  also  used  to  measure  the 
instantaneous  mass  transfer  coefficient  in  electrochemical 
reactions  [ 66 ] . 

When  a  current  is  passed  through  an  aqueous 
solution  of  a  supporting  electrolyte,  transfer  of  dissolved 
oxygen  in  the  water  takes  place  towards  the  cathode  surface. 
This  current  which  is  called  cathodic  current,  increases 
with  applied  voltage  across  the  electrodes  until  the  cat¬ 
hode  surface  becomes  completely  polarized.  Thereafter  the 
current  reaches  a  limiting  value  called  the  diffusion 
current .  In  this  range  of  applied  voltage  the  transfer  of 
oxygen  is  diffusion  controlled.  The  magnitude  of  this 
current  increases  with  increasing  velocity  of  the  liquid. 


- 

' 

' 

1 

-  ;  -  ;  joi  3  "  -  a  >  -  ■  ~J  -  0  '  n ' 


50 


If  the  applied  voltage  is  maintained  at  a  value  correspond¬ 
ing  to  the  diffusion  current,  point  liquid  velocities  may 
be  measured  by  following  the  change  in  the  magnitude  of 
the  diffusion  current.  A  sample  set  of  curves  for  dif¬ 
ferent  liquid  velocities  in  the  full  applied  voltage  range 
is  given  in  Figure  4.  The  chemical  reactions  that  take 
place  at  the  cathode  are  as  follows: 

The  cathode  reaction  is 


02  +  2H20  +  2e  -*  H202  +  20H  - (3-22) 

H202  +  2e  +2  OH" - (3.23) 


The  overall  reaction  is 

02  +  2H20  +  4e  ->  40H~ - (3.24) 

According  to  Faraday’s  Law  a  general  relationship  may  be 
written  for  the  electrochemical  current  at  a  dropping 
mercury  or  solid  stationary  electrode. 

* 

i  =  iiF  | - (3.25) 

Where  i  is  the  diffusion  current,  n  the  number  of  electrons 
taken  up  (in  this  case  4)  and  F  the  Faraday  Number.  Since 
the  current  is  the  diffusion  current,  the  rate  of  transfer 
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(Vrt)  1U3J.JL  no 


Applied  Voltage  (volts) 

Figure  ^ :  Plot  of  Polarographic  Current  vs .  Applied  Voltage 


52 


dN* 


is  the 


dt 

cathode  by 


given  by  A 


transfer  of  the 
diffusion  only, 
then  by  Fick's 


depolarisor  per  unit  area  of 
If  the  cathode  surface  area  is 
first  law  of  diffusion 


* 

dN  = 


D  A 


3C 

3x 


dt 


(3.26) 


Where  x  is  the  distance  from  the  electrode. 

Substituting  in  3.25  the  following  equation  for 
diffusion  current  may  be  obtained. 


i  =  n  F  A 


D  (— ) 
v3x;x=0 


(3.27) 


tt— .the  concentration  gradient  of  oxygen  in  water,is  evaluated 

d  X  3 

from  Fick’s  second  law  given  by  an  expression  of  the  type 


3C 

at 


D 


S2C 

3x2 


(3.28) 


Equation  3.28  is  subject  to  change  depending  on  the  shape 
of  the  electrode. 

The  initial  boundary  conditions  for  solution 
of  the  above  given  equations  are 

C  ,  =  C*  for  t  =  0 

X  ^  13 


c 


x  3 1 


C 

o 


for  t  >  0 . 
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where  C*  is  the  concentration  of  the  oxygen  in  the 
bulk  of  the  liquid  and  CQ  is  the  concentration  of  oxygen 
on  the  cathode  surface  and  is  zero  for  the  limiting  dif¬ 
fusion  current.  The  final  solution  is  of  the  form 

i  =  n  F  Ae  =  n  F  Ae  k£eC* - (3.29) 

/[Dt 

where  k„  is  the  mass  transfer  coefficient. 

Joe 

The  mass  transfer  coefficient  k„  is  a  function  of 
the  dimensions  of  the  electrode,  the  properties  of  the 
electrode  material  and  of  the  liquid  velocity.  For  a 
given  electrode  the  mass  transfer  coefficient  is  determined 
by  the  value  of  the  liquid  velocity,  and  the  limiting 
diffusion  current  may  therefore  be  measured  as  a  function 
of  liquid  velocity. 
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CHAPTER  4 


EXPERIMENTAL 
4 . 1  General  Approach 

The  purpose  of  this  investigation  was  to  measure 
the  liquid  phase  mass  transfer  coefficient  of  carbon- 
dioxide  in  water  at  bubble  formation.  The  method  employed 
was  to  determine  the  difference  between  the  volume  of  gas 
fed  to  an  orifice  and  the  product  of  the  number  of  bubbles 
and  their  mean  detached  volume.  Due  to  the  fact  that  the 
mass  transfer  from  air  to  water  in  an  air-water  system  is 
negligible  the  difference  mentioned  above  should  be  es¬ 
sentially  zero.  The  validity  and  accuracy  of  the  exper¬ 
imental  technique  was  thus  determined  using  the  air-water 
system. 

The  volumetric  flow  rate  of  a  gas  through  the 
orifice  can  be  related  to  the  bubble  volume  and  frequency 
by  the  following  equation. 

G  =  vV - (4.1) 

Where  G  is  the  gas  flow  rate  In  cc/sec. ,  v  the  bubble 
frequency  in  sec"'*'  and  V  the  volume  equivalent  of  one 
bubble  in  cc . 


For  both  the  carbon-dioxide-water  and  the  air-water 
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systems,  the  gas  flow  rate  and  bubble  frequency  were  adjusted 
to  the  desired  value  so  that  stable  bubbles  were  formed  at 
the  orifice.  High  speed  motion  pictures  of  the  bubbles 
were  then  taken  during  formation.  This  was  done  for  a 
range  of  frequencies.  The  motion  pictures  were  then 
analyzed  with  a  digitizer  and  the  surface  area  and  volume 
determined . 

If  the  volume  of  the  detached  bubble  evaluated  from 

*  * 

photographs  be  designated  as  V  then  (V-V  )  should  be  zero 

for  the  air-water  system  assuming  that  the  solubility  of 

air  in  water  is  negligible.  This  was  the  criteria  used  to 

ascertain  whether  or  not  the  experimental  techniques  used 

* 

are  feasible.  For  the  carbon-dioxide-water  system  (V-V  ) 
gives  a  measure  of  the  mass  transfer  to  the  liquid  phase 
and  using  this  information  the  mass  transfer  coefficient 
can  be  evaluated. 

4-2  Experimental  Equipment 

A  schematic  diagram  of  the  experimental  set  up  is 
given  in  Figure  5 . 

Carbon  dioxide  gas  from  a  high  pressure  cylinder 
(1)  is  saturated  with  water  In  a  packed  saturator  (2)  and 
then  metered  through  a  soap  film  meter  (3).  The  gas 
pressure  at  the  soap  film  meter  is  regulated  by  a  pressure 
regulating  valve  and  recorded  on  a  pressure  gauge  (5).  A 
micrometer  needle  valve  (4)  is  installed  in  the  gas  line 
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so  that  small  gas  flow  rates  may  be  accurately  metered. 

The  gas  then  passes  through  an  adjustable  chamber  volume 
(6)  to  the  orifice  (7)  installed  in  the  bubble  column  (8) 
where  it  forms  bubbles  in  demineralized  water  which  is  the 
continuous  liquid  phase. 

The  demineralized  water  is  stored  in  a  reservoir 
tank  (11)  from  where  it  is  pumped  with  a  centrifugal 
pump  (10)  through  a  set  of  rotameters  (9)  into  the  column. 
The  flow  of  water  through  the  column  is  cocurrent  with 
the  gas . 

The  analysis  of  the  system  requires  that  the  temp¬ 
erature  of  the  liquid  be  controlled  accurately.  A  temp¬ 
erature  control  thermometer  (13)  is  used  to  activate  a 
control  circuit  when  the  water  temperature  rises  above 
the  desired  value.  When  this  happens  cooling  water  is 
circulated  from  a  refrigerated  storage  (15)  to  a  cooling 
coil  (12)  immersed  in  the  demineralized  water.  A  small 
capacity  centrifugal  pump  (14)  is  used  for  this  purpose. 

Individual  pieces  of  equipment  that  were  used  in 
this  investigation  are  described  below. 

Column  and  Accesories 

Details  of  the  bubble  column  are  given  in  Figure  6. 

The  column  was  constructed  in  two  sections  from  a 
lucite  pipe  of  20.2  cms .  I.D.  and  0.635  cms .  wall  thickness 
with  a  total  height  of  100  cms.  The  two  sections  of  the 
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Overflow  Tank- 
Bubble  Column 
Orifice  Cover 
(Teflon ) 

Orifice  Stem  (SS) 
Wire  Mesh  No.  3 
Flanges  and  Gaskets 
Glass  Beads  (5mm<}>) 
Wire  Mesh  No.  2 
Perforated  Plate 
( Lucite ) 

Tripod  Stand  (SS) 
Wire  Mesh  No.  1 
Chamber  Volume 


Figure  6 :  Details  of  the  Experimental  Column 
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column  were  flanged  together  with  stainless  steel  bolts. 
Nylon  reinforced  rubber  gaskets  were  used  between  the 
flanges  to  prevent  water  leakage.  The  bottom  section  of 
the  column  was  provided  with  four  liquid  inlets  and  a 
suitable  attachment  to  secure  the  orifice  in  position. 

An  overflow  tank  with  four  liquid  outlets  was  constructed 
on  the  top  of  the  column. 

The  orifice  was  made  of  316  stainless  steel  tubing 
0.408  cms .  inside  diameter  and  23.0  cms .  total  length.  To 
hold  the  orifice  a  suitable  holder  was  designed  to  match 
the  corresponding  socket  provided  for  it  in  the  bottom 
section  of  the  column.  A  teflon  cover  was  made  for  the 
top  of  the  orifice  so  that  it  could  be  made  non-wettable 
with  water. 

Tygon  tubing  of  about  1.0  cms.  inside  diameter 
and  15.0  cms.  long  was  connected  to  the  bottom  end  of  the 
orifice  through  a  glass  T.  This  served  as  the  chamber 
volume  for  the  system  and  could  be  varied  by  adjusting 
the  position  of  a  pinch  cock. 

Spherical  glass  beads  5  mm.  in  diameter  were  used 
to  pack  the  bottom  section  of  the  column.  A  total  packing 
depth  of  16.5  cms.  was  provided.  To  support  the  packing, 
a  circular,  perforated  lucite  plate,  three  stainless  steel 
wire  meshes  and  a  stainless  steel  wire  stand  were  used. 

The  manner  in  which  the  glass  beads  were  supported  is 
shown  in  Figure  6 . 
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Eight  sample  points  to  sample  the  liquid  in  the 
column  were  constructed  by  drilling  equispaced  holes 
through  the  column  wall  at  different  heights  from  the 
flanged  end.  A  liquid  height  of  about  80  cms .  could  thus 
be  sampled  if  desired.  Special  sample  points  in  level 
with  the  top  edge  of  the  orifice  were  also  constructed. 

Flow  Measurement 

Rotameters  were  used  to  measure  liquid  flow  rates. 
These  were 

(i)  Brooks  model  R13M-25-3  for  0  to  2 
litres/sec 

(ii)  Brooks  model  R10M-25-3  for  0  to  0.42 
litres/sec . 

The  floats  in  both  rotameters  were  made  of  stainless  steel. 

A  soap  film  meter  was  used  to  measure  gas  flow 
rates.  It  was  constructed  from  a  precision  50  cc  glass 
burette.  The  time  taken  to  displace  a  given  volume  of  gas 
through  the  burette  was  a  measure  of  the  gas  flow  rate. 

Feed  Water  Tank 

The  feed  water  tank  was  constructed  from  1/2  inch 
lucite  plate  to  hold  a  maximum  capacity  of  100  gallons. 

Six  3/4  inch  openings  were  provided  for  lines. 
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Connecting  Lines 

Liquid  lines  were  made  from  3/4  inch  inside  diameter 
pvc  pipe.  Three  types  of  lines  were  used  for  gas  flow. 

These  were,  1/4"  inside  diameter  copper  tubing,  1/4M  inside 
diameter  polyflow  tubing  and  3/8  inch  diameter  tygon  tubing. 

Valves  and  Pumps 

Rubber  diaphram  valves  were  used  to  control  liquid 
flow  rates  and  a  Whitey  Needle  Valve  Model  22RS4  316  SS  was 
used  to  control  the  gas  flow  rate. 

A  centrifugal  pump  model  882E  "Dynapump"  manufactured 
by  Crane  Co.  was  used  for  liquid  flow.  The  pump  specific¬ 
ations  were  stainless  steel  impellor,  carbon  steel  body, 

1/2  HP,  3/4  inch  suction  and  discharge  ports. 

Gas  Saturators 

Two  columns  were  provided  for  saturating  air  and 
carbon-dioxide  with  water.  They  were  constructed  from 
glass  tubes  2.5  inch  in  diameter  and  were  packed  with  5  mm 
spherical  glass  beads.  The  top  part  of  the  columns  was 
left  unpacked  to  dampen  fluctuations  in  gas  flow.  The 
columns  were  about  three  fourths  full  of  water. 

Temperature  Control  System 

In  order  to  maintain  a  constant  temperature  of  the 
water  used  in  the  bubble  column,  a  temperature  control 
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scheme  was  used.  The  layout  of  this  scheme  Is  presented 
schematically  in  Figure  5. 

The  various  pieces  of  equipment  that  were  used  to 
control  the  water  temperature  were: 

(1)  A  refrigerator  to  maintain  the  cooling  water 
temperature  to  about  0.5°c. 

(2)  A  centrifugal  pump  for  circulation  of  the  cooling 
water . 

(3)  A  cooling  coil  made  from  316  stainless  steel. 

This  was  immersed  in  the  water  reservoir. 

(4)  A  solenoid  valve  to  actuate  the  pump  when  the 
water  temperature  rose  above  the  set  point. 

(5)  A  Beckmann  thermometer  to  specify  the  set  point 
and  to  actuate  the  solenoid  valve  when  the  temp¬ 
erature  rose  above  the  desired  value. 

The  temperature  of  the  demineralized  water  could  be  controlled 
within  +  0.05°c.  by  this  control  scheme. 

Titration  Equipment 

All  titrations,  to  standardize  solutions  and  to 
measure  the  concentrations  of  carbon-dioxide  in  water,  were 
done  electrometrically  in  an  atmosphere  of  nitrogen.  The 
following  pieces  of  equipment  were  used  for  the  purpose. 

(1)  A  Coleman  Metrion  II  model  28A  AC  pH  meter 

(2)  A  Coleman  3-450,  shielded  glass,  heavy  duty,  screw 
base,  glass  electrode 
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(3)  A  reference  electrode  model  Coleman,  Calomel 

3-510,  screw  base,  capillary  type  liquid  function. 
A  photograph  of  the  titration  equipment  as  set  up  is  given 
in  Figure  7 . 

Water  Demineralisers 

Demineralised  water  was  used  throughout  the  experi¬ 
mental  investigation.  The  quality  of  the  water  was  main¬ 
tained  at  <0.1  ppm  impurity  as  NaCl.  Mixed  bed  resin 
cartridges  no.  0808  were  used  in  a  Double  Chamber  Type, 

Model  BD-22  Pressure  Bantam  Demineralizer  manufactured  by 
Barnstead  still  and  sterlizer  Co » 

The  demineralizer  is  shown  in  Figure  8. 

Equipment  for  Measurement  of  Liquid  Velocity  Profiles 

The  electrical  circuit  required  for  the  measurement 
of  velocity  profiles  in  the  column  is  given  in  Figure  9. 

The  magnitude  of  the  diffusion  current  in  electro¬ 
chemical  reactions  of  the  type  that  occur  in  this  system, 
is  quite  small  and  therefore  difficult  to  measure.  In 
order  to  facilitate  this  measurement,  a  high  resistance 
was  installed  in  series  with  the  cathode  and  the  voltage 
across  this  resistance  was  measured  by  using  a  digital 
voltmeter.  The  magnitude  of  this  voltage  was  taken  to  be 
a  measure  of  the  point  liquid  velocity  in  the  column.  In 
order  to  obtain  the  absolute  value  of  the  latter  the  plat- 
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Figure  7:  Photograph  of  Titration  Equipment 
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Photograph  of  Demineralizer 
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Figure  9:  Circuit  Diagram  for  Measurement 
of  Liquid  Velocity  Profile 
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inum  probes  were  calibrated  to  relate  the  digital  voltage 
with  liquid  velocity.  Details  of  the  equipment  used  for 
the  measurement  of  velocity  profiles  are  given  below: 

Digital  voltmeter  -  Series  4800  digital  voltmeter,  four 

figure  display  made  by  Non-Linear 
Systems  Inc. 

Platinum  probe  (cathode)  -  A  platinum  wire  of  diameter 

0.02184  inches  shielded  and 
insulated  from  outer  casing. 

This  was  prepared  in  the  lab¬ 
oratory.  Details  of  its  con¬ 
struction  are  given  in  Figure  10. 
Platinum  foil  (Anode)  -  A  thin  one  inch  square  platinum 

foil . 

Dry  cells  -  6  volts  equivalent. 

Vacuum  tube  voltmeter  (VTVM)  -  Heathkit  model  MM1  -  The 

Heathkit  Company. 

Set  of  Variable  Resistances  (potent ionmeters )  -  100  ohms, 

1000  ohms  and  10  kilo  ohms. 

The  platinum  probes  were  calibrated  by  rotating 
them  in  a  tank  of  water  at  a  known  velocity  and  measuring 
the  digital  voltage  created.  The  electrical  circuit  used 
is  given  in  Figure  9  and  details  of  the  tank  are  shown  in 
Figure  11. 
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LEGEND 

1  Stainless  Steel 
Tubing 

2  Insulation  between 
Tubing  and  Lead 
Wire 

3  Copper  Lead  Wire 

4  Solder  between 
Copper  and  Platinum 
Wires 

5  Platinum  Wire 

6  Platinum  Wire  Tip 

7  Stainless  Steel 
Casing 

8  Holes  in  Casing 
for  Liquid  Flow- 
3  in  Number 


Figure  10:  Details  of  the  Platinum  Probe 
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Figure  11:  Photograph  of  the  Calibration  Tank 


70 

Equipment  for  Measurement  of  Bubble  Frequency. 

The  frequency  of  bubble  formation  was  measured  by 
using  an  electrical  circuit  as  shown  in  Figure  12.  A  laser 
beam  was  used  to  intercept  the  path  of  bubbles  and  a  photo¬ 
cell  was  used  to  generate  an  electrical  signal.  The  equip¬ 
ment  used  included  the  following. 

Laser  Beam  -  Spectra  Physics  model  130B  (ruby  laser) 
Photocell  -  Clairex  photocell  No.  CL  707HM 
Digital  Counter  -  Universal  counter  model  361-R, 

Transistor  Specialities  Inc. 
Oscilloscope  -  Heathkit  model  EUW-25-  The  Heathkit 

Company . 

Standard  Voltage  Source  -  Heathkit  model  EUW-16 
Vacuum  tube  voltmeter  -  Heathkit  model  MM1 
Variable  resistance  -  10  Kilo  ohms. 

Photographic  Equipment . 

A  schematic  diagram  of  the  photographic  set  up  is 
given  in  Figure  13. 

The  gas  bubbles  were  photographed  at  camera  speeds 
ranging  from  120  to  1200  frames  per  second  using  a  Hycam 
high  speed  motion  picture  camera  model  K2004E.  This  is  a 
16  mm  camera  utilizing  a  rotating  prism  for  optical  comp¬ 
ensation.  It  has  a  film  capacity  of  400  feet  of  standard 
thickness  film  but  adjustments  for  accomodating  100  feet 
and  200  feet  film  spools  are  also  available.  The  frame 
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STD  Voltage  Source 


Figure  12:  Circuit  Diagram  for  Measurement 
of  Bubble  Frequency 
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Figure  13:  Schematic  Layout  of  the 
Photographic  Equipment 
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rate  is  variable  from  0  to  11000  frames  per  second. 

A  Milli-Mite  timing  light  generator  model  TLG-3 
was  connected  to  the  camera  and  was  used  to  determine  the 
exact  speeds  of  the  film.  This  is  a  battery  operated 
transistorized  unit  designed  to  put  a  10,  100  to  1000  cycle 
time  base  on  the  moving  film. 

Illumination  of  the  bubbles  was  provided  with  two 
1500  watt  high  spot  beams.  These  lights  were  placed  to 
provide  front  and  back  lighting  of  the  object  as  shown  in 
Figure  13.  Suitable  light  diffusers  were  provided  to  dif¬ 
fuse  the  light  falling  on  the  object,  and  heat  sheilds 
were  used  to  protect  the  column  against  thermal  expansion. 
An  aluminum  foil  reflector  was  used  to  facilitate  lighting 
of  the  bottom  surface  of  the  bubble. 

In  order  to  prevent  distortion  of  the  bubble  due 
to  the  curvature  of  the  column,  a  correction  lens  was 
designed  and  constructed.  Details  of  the  ray  diagram  and 
the  correction  lens  are  given  in  Figure  14 .  The  lens  con¬ 
sists  essentially  of  a  hollow  plano-concave  lucite  casing 
filled  with  water. 

Eastman  Tri-X  reversal  film  type  7278  was  used 
throughout  the  investigation. 

The  Digitizer. 

Figure  15  is  a  photograph  of  the  equipment  used  to 
measure  the  bubble  dimensions  from  photographic  film. 
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Lucite  Column  Wall 
Circular  Lucite  Wall 
Water  in  Correction  Lens 
Flat  Lucite  Plate 
Air  Atmosphere 
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Figure  15:  Photograph  of  the  Digitizer 
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The  film  is  projected  on  a  ground  glass  screen 
using  a  model  224-A  LW  Photo  Optical  Data  Analyzer.  This 
is  a  16  mm  stop-start  projector  with  a  variable  film  speed. 
Its  projection  speed  and  stop-start  action  can  be  controlled 
with  a  remote  control  panel. 

The  projected  image  of  the  bubble  can  be  scanned 
manually  using  a  cross  wire  and  any  point  (x,y)  on  the 
bubble  periphery  can  be  represented  by  two  voltages  for 
the  x  and  y  coordinates  respectively.  These  voltages  can 
be  read  and  punched  on  computer  cards  in  a  pre-determined 
format . 

4 . 3  Experimental  Procedure 

The  experimental  procedure  may  be  discussed  under 
the  following  subheadings : 

(1)  Cleaning  and  assembly  of  equipment 

(2)  Photographic  procedure 

(3)  Experimental  runs 

(4)  Analysis  of  film 

(5)  Precautions 

1 .  Cleaning  and  assembly  of  equipment 

In  order  to  obtain  reliable  and  reproducible  results 
for  both  the  bubble  behaviour  and  mass  transfer  it  was 
important  that  the  system  be  cleaned  carefully  and  thoroughly 
to  avoid  contamination .  Particular  care  was  taken  to  avoid 
any  oily  material  from  getting  into  the  system.  To  do  this. 
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rubber  gloves  were  used  to  handle  all  the  equipment  during 
the  process  of  cleaning  and  assembly. 

Glass  beads  used  for  packing  the  bottom  section  of 
the  column  were  cleaned  by  immersing  them  in  concentrated 
chromic  acid  for  at  least  a  period  of  forty  eight  hours. 

They  were  then  rinsed  repeatedly  with  demineralized  water 
and  left  to  dry.  The  orifice  tube,  its  teflon  cover,  the 
rubber  gaskets  and  the  wire  meshes  were  washed  thoroughly 
using  a  detergent  and  were  then  rinsed  with  water  and  dried 
with  air.  In  order  to  improve  the  non-wetting  characteristics 
of  the  teflon  cover,  it  was  left  immersed  in  acetone  for  a 
short  period  of  time.  The  main  column,  reservoir  tank  and 
connecting  lines  were  cleaned  by  recycling  demineralized 
water  with  detergent  through  the  system.  A  period  of  about 
three  to  four  hours  was  considered  to  be  sufficient  for  this 
purpose.  Fresh  batches  of  demineralized  water  were  used  to 
rinse  the  system  free  of  detergent.  This  was  done  for  as 
many  as  ten  or  fifteen  times,  and  each  batch  of  water  was 
drained  after  a  maximum  of  one  hour. 

The  importance  of  careful  cleaning  cannot  be  over 
emphasized.  Some  runs  were  conducted  where  adequate  care 
was  not  devoted  for  this  purpose  and  it  was  found  that 
unstable  bubbles  were  formed  at  the  orifice  giving  a  poor 
reproducibility  of  bubble  frequency. 

Details  of  the  assembly  of  the  column  and  its 
accessories  are  given  in  Figure  6.  Experience  with  the 
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system  has  shown  that  particular  attention  should  be  paid 
to  certain  aspects  of  column  assembly.  These  include  the 
following : 

(a)  In  order  to  avoid  air  entrappment  in  the  packing,  the 
glass  beads  should  be  introduced  into  the  column  by  dis¬ 
placement  of  water.  Entrapped  air  interferes  with  measure¬ 
ment  of  bubble  frequency  and  tends  to  give  a  non-uniform 
liquid  velocity  profile. 

(b)  For  the  same  reason  as  in  (a)  the  wire  meshes  should 
be  wetted  from  the  bottom  up  when  the  column  is  filled  with 
water . 

(c)  The  glass  beads  should  be  uniformly  distributed  through 
the  entire  cross-section  of  the  column.  This  is  done  to 
ensure  a  reasonably  flat  liquid  velocity  profile. 

(d)  Before  the  rubber  gaskets  are  placed  in  position,  both 
flanges  should  be  dried.  This  prevents  water  leakage  from 
the  flanges. 

(e)  A  torque  wrench  should  be  used  to  tighten  the  nuts 
and  bolts  holding  the  two  flanges  together.  Diametrically 
opposite  bolts  should  be  tightened  in  sequence.  This 
ensures  an  even  pressure  on  the  rubber  gaskets  and  prevents 
water  leakage. 

2 .  Photographic  procedure. 

Figure  13  gives  a  schematic  arrangement  of  the 
photographic  set  up. 

The  correction  lens  (9)  was  secured  tightly  with 
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the  wall  of  the  column  (7)  eliminating  any  air  gap  between 
the  lens  and  the  column  wall  as  far  as  possible.  The 
aluminum  foil  reflector  (12)  was  then  wrapped  around  the 
base  of  the  column  in  such  a  way  that  the  square  hole  (11) 
was  in  line  with  the  top  of  the  orifice  and  the  gas  bubble. 
Two  high  spot  lights  (1,  2)  were  placed  in  position  so  as 
to  provide  back  and  front  lighting  of  the  gas  bubble  as 
shown  in  Figure  13.  A  light  collimater  (4)  was  adjusted 
between  the  high-speed  camera  (3)  and  the  correction  lens 
so  that  extraneous  light  from  the  environments  could  be 
prevented  from  getting  into  the  camera.  The  camera  was 
mounted  on  a  suitable  stand  (not  shown  in  the  figure)  and 
its  level  was  adjusted  so  that  the  camera  lens  was  in 
direct  line  with  the  object.  The  camera  stand  was  then 
locked  in  that  position.  Light  shields  (5,6)  and  the  light 
diffuser  (13)  were  then  provided  so  as  to  completely  mask 
the  column  from  direct  lighting. 

The  following  procedure  was  used  to  photograph  the 
bubbles  during  the  course  of  an  experimental  run. 

The  camera  was  threaded  with  a  400  ft.  roll  of 
16  mm  film  and  the  millimite  cycle  timer  was  hooked  into 
the  receptacle  provided  for  it  in  the  camera.  The  desired 
timing  frequency  and  the  photographing  speed  were  then  set 
on  the  cycle  timer  and  camera  respectively.  Both  high  spot 
lights  were  turned  on  and  the  light  intensity  was  measured 
through  the  collimater  by  using  a  suitable  exposuremeter . 
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The  camera  was  then  accurately  focussed  on  the  orifice 
keeping  the  lens  aperture  wide  open  and  making  sure  that 
the  correction  lens  is  full  of  water.  The  lens  aperture 
was  then  adjusted  according  to  the  reading  of  the  exposure- 
meter.  To  prevent  any  motion  of  the  camera  during  filming, 
it  was  locked  in  position  after  the  above  mentioned  setting. 
After  adjusting  the  variables  pertinent  to  the  experiment 
to  their  desired  values  the  camera  switch  was  turned  on 
and  the  desired  footage  of  film  was  exposed. 

3 .  Experimental  runs. 

To  carry  out  an  experimental  run  with  air  or  carbon 
dioxide  the  following  steps  were  followed: 

(a)  Admit  gas  (air  or  CC^)  into  the  system  through  the 
gas  line  at  a  desired  gas  pressure. 

(b)  Admit  water  into  the  column  and  adjust  the  rate  of 
flow  of  water. 

(c)  Control  liquid  temperature  to  the  desired  value. 

(d)  Adjust  the  gas  flow  rate  and  obtain  a  constant  bubble 
frequency . 

(e)  Make  final  adjustments  of  the  gas  pressure. 

(f)  Analyse  liquid  samples  of  CO^-water  runs  only. 

(g)  Take  photographs. 

The  supply  gas  pressure  was  adjusted  to  the  desired 
value  and  the  gas  was  admitted  into  the  saturator  (see 
Figure  5)  slowly  and  intermittently  to  prevent  a  surge  of 
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liquid  through  the  gas  line.  The  control  valve  in  the 
gas  line  was  set  to  give  the  desired  pressure  of  the  gas 
at  the  soapfilm  meter,  and  the  micrometer  needle  valve  was 
opened  to  let  the  gas  into  the  column  through  the  orifice. 

Any  air  in  the  water  lines  and  rotameter  tubes  was 
then  expelled  from  the  system  through  the  air  vent.  To  do 
this,  the  water  inlet  ports  to  the  column  were  closed  and 
air  was  displaced  through  the  vent  by  pumping  water  through 
the  lines.  The  column  was  then  filled  with  water  at  a  very 
low  flow  rate  and  the  latter  was  adjusted  to  the  desired 
value  once  the  column  was  full. 

The  set  point  of  the  Beckman  Thermometer  was  fixed 
at  the  desired  temperature  and  the  temperature  control 
circuit  was  turned  on.  The  flow  of  cooling  water  through 
the  cooling  coil  was  adjusted  till  the  on-off  control  was 
set  in  operation.  The  liquid  temperature  was  thus  controlled 
within  a  limit  of  +  0.05°c. 

In  order  to  measure  bubble  frequencies,  a  laser 
beam  and  photocell  were  used.  The  electrical  circuit  for 
this  measurement  is  given  in  Figure  12.  The  position  of 
the  laser  beam  was  adjusted  such  that  the  bubble  intercepted 
the  beam  after  detachment  from  the  orifice.  In  order  to 
avoid  erroneous  measurements  at  high  gas  flow  rates  due  to 
incipient  coalescence,  the  laser  beam  position  was  adjusted 
away  from  the  central  axis  of  the  bubble.  Every  interruption 
of  the  laser  beam  was  measured  as  a  descrete  signal  by 
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the  photocell  and  was  recorded  on  a  universal  counter.  An 
oscilloscope  was  used  to  display  a  waveform  which  was 
indicative  of  the  stability  of  bubble  formation. 

The  gas  needle  valve  and  the  gas  pressure  were 
adjusted  alternately  till  the  desired  bubble  frequency  and 
gas  pressure  were  obtained. 

For  mass  transfer  runs  with  CO^-water  system,  the 
continuous  water  phase  has  to  be  analysed  for  the  amount  of 
dissolved  CO^  .  This  is  necessary  to  determine  the  driving 
force  for  mass  transfer. 

The  analysis  of  CO^  in  water  was  done  by  back 
titrations.  A  blank  titration  with  Ba(OH)2  against  HC1  in 
an  atmosphere  of  nitrogen  was  done  for  each  experimental 
run.  To  titrate  the  water  sample,  a  known  volume  of  it  was 
added  to  a  known  volume  of  Ba(OH)2  solution  and  the  mixture 
was  titrated  with  HC1  to  a  green  end  point  of  Bromo  Thymol 
Blue  indicator.  A  nitrogen  atmosphere  was  used  throughout 
the  titration.  The  concentration  of  CO^  in  water  could 
then  be  determined  by  calculation  as  indicated  in  Appendix 
C.  The  water  sample  was  titrated  either  immediately  before 
or  immediately  after  photographing  the  bubbles. 

The  bubbles  were  then  photographed  by  using  the 
procedure  outlined  in  section  2.  The  values  of  all  the 
variables  pertinent  to  the  experiment  were  recorded  before 
and  after  filming.  In  order  to  determine  the  number  of 
bubbles  to  be  analysed  to  obtain  a  reliable  value  of  bubble 
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surface  and  volume,  it  was  necessary  to  record  a  represent¬ 
ative  population  of  bubble  frequency  during  each  experimental 
run. 

^ .  Analysis  of  films. 

The  films  of  gas  bubbles  taken  during  the  experimental 
runs  were  analysed  in  the  following  fashion. 

(a)  Film  Editing 

(b)  Determination  of  the  number  of  data  points  per  bubble. 

(c)  Analysis  of  detached  bubbles. 

(d)  Analysis  of  bubbles  during  formation. 

Before  any  data  could  be  taken  from  the  films,  they 
were  screened  to  determine  the  quality  of  the  pictures 
obtained.  A  given  film  was  accepted  only  if  it  conformed 
to  the  following  specifications. 

(i)  The  orifice  should  have  sharp  well  defined  edges, 

(ii)  The  bubble  should  be  sharp  and  completely  in  focus 
at  each  stage  of  formation. 

(iii)  There  should  be  a  deep  contrast  of  the  orifice 
and  bubble  with  its  background. 

(iv)  There  should  be  no  light  reflections  on  the  bubble 
edges  at  any  time  during  formation  and  detachment. 

If  the  pictures  failed  to  meet  any  one  of  the  above 
conditions,  it  was  discarded  and  the  experimental  run  was 
repeated . 

The  total  number  of  bubbles  photographed  during  a 
run  were  then  counted  and  a  given  number  of  them,  as  required 
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by  the  statistical  sample  size,  were  isolated.  A  reference 
mark  was  made  on  the  film  at  a  position  corresponding  to 
the  first  bubble  to  be  analysed.  This  was  done  so  that  the 
frame  counter  on  the  film  projector  could  be  set  at  zero. 

In  order  to  obtain  reliable  results  of  bubble  sur¬ 
face  area  and  volume  it  was  necessary  to  decide  on  the 
number  of  data  points  to  be  taken  from  the  periphery  of 
one  bubble.  This  was  done  by  a  trial  and  error  procedure. 

The  number  of  data  points  taken  from  a  given  bubble  were 
increased  for  each  trial  till  the  bubble  surface  area  and 
volume  evaluated  for  two  consecutive  trials  did  not  differ 
by  more  than  an  acceptable  tolerance.  This  analysis  was 
done  for  detached  bubbles  only.  For  forming  bubbles  the 
number  of  data  points  had  to  be  increased  as  the  bubble 
grew  larger.  The  maximum  number  of  data  points  were  taken 
just  prior  to  detachment. 

Analysis  of  detached  bubbles  involved  taking  data 
from  a  single  frame  corresponding  to  the  bubble  in  its 
detached  position.  This  was  done  for  each  bubble  of  a 
given  sample  size.  In  order  to  obtain  surface  area  and 
volume  as  a  function  of  time,  however,  different  frames  of 
the  same  bubble  were  analysed  at  various  stages  of  formation. 
Analysis  of  all  bubbles  in  a  given  sample  size  would  there¬ 
fore  require  taking  data  from  a  prohibitively  large  number 
of  frames.  To  avoid  this, a  representative  bubble  from  a 
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sample  was  selected  and  the  surface  area  as  a  function  of 
time  evaluated  from  it  was  used  to  determine  the  mass  tran¬ 
sfer  coefficient. 

Data  from  the  bubble  periphery  was  taken  with  the 
digitizer.  Taking  the  data  consisted  of  reading  the  co- 
odinates  of  points  on  the  periphery  and  punching  these  on 
computer  cards  in  a  pre-determined  format. 

5  .  Precautions . 

Careful  experimental  work  is  prerequisite  to  obtaining 
reliable  data  from  an  experimental  program.  Success  of  the 
experimental  work  invariably  depends  on  the  care  taken  in 
observing  a  number  of  precautions  pertinent  to  a  particular 
system.  Some  precautions  that  should  be  observed  while 
performing  an  experiment  such  as  outlined  here,  are  re¬ 
iterated  below: 

(a )  Cleaning  and  assembly  of  equipment. 

(1)  The  quality  of  the  demineralised  water  should  not  fall 
below  the  desired  value. 

(2)  Handling  the  components  inside  the  column  with  bare 
hands  should  be  avoided. 

(3)  If  the  system  is  washed  with  detergent,  adequate  care 
should  be  taken  to  wash  all  the  detergent  away. 

(4)  Some  extraneous  material  can  sometimes  collect  in  the 
orifice  passage.  This  situation  should  be  avoided.  The 
orifice  cover  should  be  finally  cleaned  with  acetone. 
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(5)  The  rubber  gaskets  and  wire  meshes  should  be  dried  with 


air  before  being  used. 

(6)  Glass  beads  should  be  evenly  packed  in  the  base  of  the 
column  taking  care  to  avoid  entrapped  air  bubbles. 

(7)  The  column  flanges  should  be  evenly  tightened  to  avoid 
liquid  leakage . 

(b )  Photgraphic  Procedure. 

(1)  The  column  should  be  kept  well  shielded  from  direct  light. 

(2)  The  camera  should  be  focussed  on  the  object,  with  the 
lens  aperture  fully  open  to  take  full  advantage  of  the  light. 
The  desired  aperture  for  filming  should  subsequently  be 

adj  usted . 

(3)  Care  should  be  taken  to  ensure  that  there  is  water  in 
the  correction  lens  before  filming  the  bubbles. 

(c )  Experimental  Run. 

(1)  The  experimental  run  should  be  performed  as  soon  as 
the  system  is  cleaned  and  assembled  to  avoid  deviations  in 
data  resulting  from  a  change  in  the  surface  age  of  the  orifice. 

(2)  All  gas  leaks  in  the  gas  line  should  be  completely 
eliminated  before  starting  a  run.  Any  gas  leak  could  be 
interpreted  as  additional  mass  transfer  and  this  could 
jeopardise  the  purpose  of  the  investigation. 

(3)  Thermal  effects  to  the  column  from  its  surroundings 
should  be  minimised  by  adequate  shielding. 
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(4)  The  flow  of  gas  through  the  orifice  should  be  started 


before  the  column  is  filled  with  the  liquid. 

(5)  Care  should  be  taken  to  avoid  moving  any  gas  lines 
during  the  experiment.  This  affects  bubble  frequency. 

(6)  To  prevent  desorption  of  CO^  from  the  water  sample, 
the  pipette  tip  should  dip  into  the  BaCOfOp  solution  before 
the  sample  is  discharged. 

(d )  Analysis  of  films. 

(1)  Film  editing  should  be  done  very  critically  before 
deciding  to  take  data. 

(2)  The  number  of  points  read  from  the  bubble  periphery 
should  be  evenly  distributed  for  good  results. 

(3)  All  instructions  for  proper  use  of  the  digitizer 
should  be  closely  followed. 
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RESULTS 

5 . 1  General 

The  results  of  this  investigation  are  subdivided 
as  follows: 

(i)  Liquid  Velocity  Profiles 

(ii)  Air-Water  System 

(iii)  Mass  transfer  during  bubble  formation  with  the  CO^ 
water  system. 

5 . 2  Liquid  Velocity  Profiles 

The  glass  bead  packing  and  the  stainless  steel  wire 
meshes  in  the  bottom  section  of  the  column  were  intended  to 
provide  a  flat  liquid  velocity  profile  in  the  column.  It 
was  found,  however,  that  due  to  varying  degrees  of  chan¬ 
nelling  from  different  sections  of  the  packing,  the  velocity 
profile  was  different  across  different  diameters  of  the 
column.  Measurements  of  point  velocities  were  thus  made  to 
select  the  diameter  along  which  minimum  deviation  from  a 
flat  profile  could  be  obtained.  Thereafter,  velocity 
measurement  were  made  along  the  selected  diameter  at  dif¬ 
ferent  average  liquid  velocities,  so  as  to  select  the 
liquid  velocity  which  gave  the  most  flat  profile.  These 
measurements  were  made  for  an  orifice  size  of  2.1  mm,  and 
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bubble  frequencies  ranging  from  0  to  25  per  second.  All 
measurements  were  made  using  the  air-water  system.  The 
velocity  profiles  obtained  for  average  velocities  of  1.0 
and  2.0  cms/sec  are  given  in  Figures  16,  17  &  18.  The 
data  is  reported  in  Appendix  C. 

5 . 3  Air-Water  System 

The  air-water  system  was  used  to  investigate  the 
following : 

(i)  The  inter-relationship  of  orifice  diameter,  bubble 
frequency,  gas  flow  rate  and  bubble  diameter.  This  may  also 
be  referred  to  as  a  study  of  bubble  phenomena. 

(ii)  An  evaluation  of  the  technique  used  to  measure 
mass  transfer  during  bubble  formation. 

(iii)  An  observation  of  bubble  behaviour  during  bubble 
formation . 

Investigations  of  the  effect  of  orifice  diameter, 
bubble  frequency,  gas  flow  rate  and  bubble  diameter  on 
bubble  phenomena  have  been  carried  out  by  a  number  of 
previous  investigators.  A  further  investigation  was  war¬ 
ranted  here  in  order  to  establish  a  conformity  of  funda¬ 
mental  bubble  behaviour  with  previous  work  in  the  field. 

The  range  of  experimental  conditions  for  the  air- 
water  system  are  given  in  Table  1. 
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Figure  17 :  Liquid  Velocity  Profiles  at  Different 

Bubble  Frequencies  for  an  average  liquid 
velocity  of  1.0  cms/sec. 
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Figure  18 :  Liquid  Velocity  Profiles  at  different 

Bubble  Frequencies  for  an  average  liquid 
velocity  of  2.0  cms/sec . 
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Inside  Orifice 
Diameter 
(cms ) 

Gas  Flow  Rate 
(cc/sec ) 

Bubble  Frequency 
( 1/sec ) 
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i — i 
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o 
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^r 

i — 1 

* 
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0 . 408 

0  +  50.0 

0  +  25.0 

0.210 

o 

• 

I — 1 

i — 1 

* 

o 

0  -*>  26.0 

0.0936 

0  +  2.0 

o 

CO 

C\J 

+ 

o 

0.0485 

0  -+  1.0 

0  -+  45.0 

Table  1:  Range  of  Experimental  Conditions  for 


The  Air-Water  System 

An  equivalent  spherical  bubble  diameter  was 
evaluated  from  the  bubble  volume  at  detachment  and  this  is 
plotted  against  the  Reynolds  number  based  on  the  orifice 
diameter  in  Figure  19.  This  investigation  was  conducted 
in  a  stagnant  liquid  medium.  The  following  significant 
conclusions  may  be  drawn  from  the  data  of  Figure  19. 

(i)  For  the  smallest  orifice  size  the  bubble  formation 
occurred  completely  in  the  constant  volume  region. 

(ii)  For  the  larger  orifice  sizes,  the  bubbles  were 
formed  in  both  the  constant  volume  and  the  transition 
regions . 

(iii)  At  a  given  Reynolds  number  the  bubble  diameter 
increased  with  increasing  orifice  size. 

The  existence  of  the  transition  region  cannot  be 
established  through  Figure  19  only  but  must  be  confirmed 
by  reference  to  Figure  20  also.  This  is  due  to  the  fact 
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Figure  19:  Plot  of  Bubble  Diameter  vs.  Orifice  Reynolds  Number 
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Figure  20:  Plot  of  Bubble  Frequency  vs.  Orifice  Reynolds 
Number 
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that  bubble  volume  increases  with  increasing  gas  flow  rate 
in  the  constant  frequency  and  the  transition  regions. 

Figure  20  is  a  plot  of  bubble  frequency  and  Reynolds  number 
and  for  the  range  of  gas  flow  rates  studied  it  indicates 
that  the  constant  frequency  region  was  not  reached.  This 
confirms  the  existence  of  the  transition  zone.  It  is  also 
evident  from  Figure  20  that  the  correlation  is  independent 
of  orifice  diameter. 

These  observations  are  in  agreement  with  the  work 
of  previous  authors. 

The  air-water  system  was  used  to  characterize  the 
orifice  and  to  establish  the  effectiveness  of  the  technique 
adopted  to  evaluate  bubble  surface  areas  and  volumes.  It 
was  assumed  that  the  mass  transfer  from  the  air  bubble  to 
the  water  phase  is  negligible.  The  bubble  volume  before 
and  after  formation  should  thus  remain  unchanged.  A  com¬ 
parison  was  made  between  the  volume  of  the  air  bubbles  as 
evaluated  from  the  gas  flow  rate  and  bubble  frequency  with 
the  volume  evaluated  after  bubble  formation  from  photographic 
films.  Figure  21  shows  this  comparison.  Though  there  is 
some  scatter  about  the  45°  line,  the  variation  is  distributed 
on  either  side  and  the  accuracy  of  the  data  is  within 
experimental  limitation  of  -  2%. 

The  behaviour  of  gas  bubbles  during  formation  was 
observed  from  photographic  films.  The  photographs  of 
Figure  22  show  a  typical  bubble  during  its  various  stages 
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Figure  21:  Comparison  of  Bubble  Volume 

from  Gas  flow  rate  and  bubble 
frequency  with  Bubble  volume 
from  Photographs 
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Figure  22:  Photographs  of  Bubble  Behaviour 
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of  formation. 

In  its  initial  stages  of  formation,  the  gas  bubble 
experiences  a  mechanical  oscillation  in  the  axial  direction 
resulting  from  an  equal  and  opposite  reaction  due  to  the 
break  off  of  the  previous  bubble.  These  oscillations  cease 
fairly  quickly  and  the  bubble  grows  steadily  and  uniformly 
in  an  almost  spherical  pattern.  This  growth  occurs  till 
the  volume  of  the  bubble  approaches  a  limiting  value  det¬ 
ermined  by  an  equilibrium  between  the  surface  and  gravitat¬ 
ional  forces.  A  neck  then  develops  and  the  bubble  begins 
oscillating  in  its  axial  direction.  The  neck  becomes  thinner 
with  each  oscillation  but  the  bubble  volume  keeps  on 
increasing.  The  bubble  eventually  breaks  off  from  the 
orifice  and  this  creates  an  imbalance  in  the  forces  acting 
on  the  bubble.  This  effect,  which  is  sometimes  referred 
to  as  a  kinetic  impingement,  is  responsible  for  the  initial 
oscillations  of  the  oncoming  bubble  and  for  severe  rippling 
on  the  surface  of  the  detached  bubble.  On  detachment,  the 
bubble  undergoes  violent  deformations  in  addition  to  surface 
rippling  till  finally  it  assumes  an  oblate  spheroidal  shape 
which  is  deformed  only  by  its  spiral  motion  up  the  liquid 
column . 

9 . 4  Mass  Transfer  During  Bubble  Formation:  CQ^-Water  System 


Experiments  were  carried  out  to  measure  the  amount 
of  mass  transfer  during  bubble  formation  for  the  absorption 
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of  CO^  in  water.  The  following  conditions  were  used: 

Orifice  Inside  Diameter  =  0.408  cms . 

Liquid  Average  Velocity  =1.0  cms/sec. 

Bubble  Frequencies  =  0.5,  1.0,  5.0  and  10.0  bubbles/sec. 
The  0.408  cms.  diameter  orifice  was  chosen  in  preference  to 
the  others  from  its  overall  performance  during  studies  with 
the  air-water  system.  Bubble  formation  was  more  stable  for 
this  orifice  size.  The  sets  of  experimental  conditions  used 
for  individual  runs  are  given  in  Table  2. 


Run  No . 

Gas  Rate 
cc/sec 

Bubble  Frequency 
( 1/sec ) 

Average 
Liquid 
Velocity 
( cms/sec ) 

Ml 

0.1011 

1.0190 

1.0 

M2 

0.5008 

5.008 

1.0 

M3 

1.3298 

10 . 002 

1.0 

M4 

0 . 0850 

1 . 0269 

1.0 

M5 

0 » 4846 

4.996 

1.0 

M6 

1.2441 

9.883 

1.0 

M7 

0.0483 

0.5014 

1.0 

Table  2:  Experimental  Conditions  used  for  CO^-Water  Runs 

The  surface  area  and  volume  of  bubbles  determined 
as  a  function  of  time  from  photographs  is  shown  in  Figures 
23  to  28.  In  evaluating  bubble  area  and  volume,  corrections 
were  incorporated  to  account  for  the  rest  bubble  which  is 
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Figure  23:  Plot  of  Bubble  Surface  Area  vs. 

Time  for  Runs  M-l,  M-2  and  M-3 
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Time  (sec) 

Figure  24:  Plot  of  Bubble  Surface  Area  vs.  Time  for 
Runs  M-4,  M-5  and  M-6 
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Figure  25:  Plot  of  Bubble  Surface  Area 
vs .  Time  for  Run  M-7 
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Figure  26:  Bubble  Volume  vs.  Time  for  Runs 
M-l,  M-2  and  M-3 
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Figure  27:  Bubble  Volume  vs.  Time  for  Runs 
M-4 ,  M-5  and  M-6 
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Figure  28: 


Bubble  Volume  vs.  Time  for 
Run  M-7 


V  -  1  1  -■ 


.  J 


107 

that  part  of  the  bubble  left  behind,  on  the  orifice  after 
detachment.  The  surface  area  corresponding  to  the  rest 
bubble  was  included  in  the  evaluation  of  the  area  vs.  time 
integral  and  the  volume  of  the  rest  bubble  was  deducted 
from  the  volume  of  the  oncoming  bubble  at  each  time  interval. 
The  necessity  of  this  correction  to  obtain  reliable  results 
was  pointed  out  by  Heertjes  et .  al .  [67]. 

The  mass  transfer  coefficient  kT  evaluated  from 
the  data  is  an  average  value  determined  during  the  time  of 
formation  0<t<T^.  The  general  expression  for  evaluation 
of  the  mass  transfer  coefficient  is  given  by  equation  3  .3 
expressed  as 

Tf 

NA  ■  kLav(CI  -  CA}  j  A<t)dt - (3.3) 

0 

where  the  significance  of  the  symbols  has  been  previously 

defined.  The  value  of  the  concentration  of  CO  in  the  bulk 

2 

liquid  phase  i.e.  was  found  to  be  essentially  zero. 

g 

Details  of  the  calculations  involved  for  evaluation  of 
and  are  given  in  Appendix  D.  The  results  obtained  for 
amount  of  mass  transfer  during  bubble  formation  and  the 
mass  transfer  coefficient  are  given  in  Tables  35  4  and  5. 
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CHAPTER  6 


DISCUSSION 

6 . 1  Air-Water  System. 

The  variation  of  bubble  frequency  with  Orifice 
Reynolds  number  as  shown  in  Figure  20  indicates  that  the 
correlation  is  independent  of  orifice  diameter  in  the 
constant  volume  and  transition  regions.  At  high  Reynolds 
numbers  it  is  seen  that  the  curve  levels  off  to  some  extent 
indicating  the  approach  of  the  constant  frequency  region. 
Davidson  et .  al.  [16]  indicated  that  the  maximum  bubble 
frequency  corresponding  to  the  constant  frequency  region  is 
a  function  of  the  orifice  diameter  and  is  given  as 


where  n  is  the  maximum  frequency  attained.  This  effect 
can  be  seen  in  Figure  20  at  high  Reynolds  numbers  where 
the  frequency  curve  tends  to  level  off  at  different  values 
for  different  orifice  diameters. 

6 . 2  CO^-Water  System 

The  interpretation  of  the  results  obtained  for 
the  CO^-water  system  may  be  subdivided  into  the  following 


sections : 
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(1)  Variation  of  Surface  Area  and  Volume  with  time  of  bub¬ 
ble  formation. 

(ii)  Mass  transfer  during  bubble  formation. 

During  a  study  of  mass  transfer  at  bubble  formation 
Popovich  et .  al .  [57]  assumed  that  the  surface  area  of  the 
growing  bubble  or  drop  may  be  represented  by  the  equation 


A  =  k  t  - (2.  19) 

or  In  A  =  In  +  n^  In  t  - (6.1) 

If  the  bubble  is  assumed  to  grow  as  a  sphere  then  n^  =  2/3. 


The  surface  area  of  the  bubbles  obtained  in  this 
investigation  is  plotted  according  to  equation  6,1  in 
Figures  29,  30  and  31.  It  is  seen  from  these  figures  that 
(In  A)  is  not  linear  with  (In  t)  throughout  the  period  of 
formation  as  indicated  by  equation  6.1.  A  departure  from 
linearity  is  observed  in  each  case,  resulting  in  a  sharp 
rise  of  surface  area  and  examination  of  photographic  films 
indicated  that  the  non-linear  portion  of  the  curve  cor¬ 
responds  to  bubble  necking  and  subsequent  bubble  break-off. 
It  may  be  concluded  from  this  that  during  neck  formation 
and  bubble  breakoff  the  surface  area  of  the  bubbles  rises 
at  a  faster  rate  than  predicted  by  equations  2.19  or  6.1. 

Calderbank  et .  al.  [8]  suggested  that  the  bubble 
area  as  a  function  of  time  could  be  represented  by  equation 
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Figure  29:  Log-Log  Plot  of  Bubble  Area  vs 
Time  for  Runs  M-l ,  M-2  and  M-3 
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Figure  31:  Log-Log  Plot  of  Bubble  Area  vs. 
Time  for  Run  M-7 


E  rpp  -  y-LGZ  (fed’  c.  is ) 


2.26  which  is  given  as  follows: 
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1- 


(T^-t)  n0 

[  — - - ]  2 

T,, 


(2.26) 


The  value  of  the  exponent  n^  was  reported  as  1.86  for  bub¬ 
ble  frequencies  between  0.2  and  2.0  sec  ^  and  equivalent 
spherical  bubble  diameters  between  0.5  and  1.0  cms . 

Using  the  values  of  A  and  T„  as  obtained  from  this 

m  f 

investigation,  equation  2-26  was  plotted  as  (In  A)  vs. 

(In  t)  and  is  given  in  Figures  29,  30  and  31.  The  mean 
slope  n^  of  the  linear  portion  of  the  plot  as  obtained 
from  this  investigation  was  found  to  be  0.850.  This  value 
compares  well  with  a  mean  slope  of  0.875  as  obtained  from 
the  plot  of  Calderbank’s  equation  [8].  It  is  possible  to 
conclude  from  this  that  the  bubble  growth  was  not  spherical 
for  the  range  of  bubble  sizes  and  bubble  frequencies 
studied.  Non  spherical  growth  was  also  observed  by 
Poutanen  and  Johnson  [68]  who  found  that  a  value  of  n^ 
varying  between  2/9  and  5/9  best  suited  their  data. 

When  plotted  on  linear  coordinates  the  sharp  rise 
in  bubble  area  at  necking  is  represented  by  a  pronounced 
inflection  as  shown  in  Figures  23  to  25.  This  effect 
though  present  at  ail  frequencies  cannot  be  clearly  seen 
at  higher  frequencies  from  Figures  23  and  24  because  of  the 
scale  chosen.  It  is  interesting  to  plot  bubble  area  vs. 
time  on  transformed  time  coordinates  taking  the  time  axis 
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to  be  (Tf-t)  where  Tf  is  the  total  formation  time.  This 
affords  a  comparison  with  the  work  of  Poutanen  and  Johnson 
[68].  The  inflection  region  is  fairly  conspicuous  from  this 
plot  also  as  is  shown  in  Figure  32,  but  the  work  of  Poutanen 
and  Johnson  [68]  indicates  no  such  observation.  This  dif¬ 
ference  may  be  due  to  the  different  techniques  used  to 
evaluate  bubble  surface  area  during  formation.  Poutanen 
et .  al.  represented  the  bubble  profile  by  an  emperical 
equation  of  the  form 

x 

rn  0  =  1 - (6  .2) 

X 

where  r  and  0  are  polar  coordinates  and  n  is  an  exponent 
determined  from  some  shape  factors  evaluated  using  specific 
bubble  dimensions.  The  expressions  for  bubble  surface  area 
and  volume  were  determined  explicity  in  a  manner  similar  to 
the  one  used  in  this  study  to  derive  expressions  given  by 
equations  3  .16  and  3.17.  The  area  and  volume  of  the  elongated 
neck  during  bubble  necking  was  evaluated  assuming  the  neck 
to  be  made  up  of  frustra  of  right  cones.  In  this  study, 
however,  the  numerical  approach  was  used  in  preference  to 
the  continuous  function  approach  because  the  neck  could  be 
accurately  represented  by  actual  data  points  before  evaluating 
the  surface  and  volume  and  also  because  this  procedure  was 
thought  to  give  more  reliable  results. 


The  variation  of  bubble  volume  with  time  was  found 


' 

3  ~  ■  ..  t  if-  :  ■■■  -  • 

" 

i 


Bubble  Area  (sq.  cms. 


118 


Figure  32:  Plot  of  Bubble  Area  vs.  Time  on 
transformed  Time  Coordinates 
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to  be  substantially  linear  throughout  the  period  of  formation. 
This  is  in  agreement  with  the  results  of  Poutanen  et .  al. 

[68]  for  the  Nitrogen-Water  System.  The  linearity  of  this 
plot  merely  indicates  that  a  steady  gas  flow  rate  was 
maintained  during  the  period  of  formation. 

The  general  expressions  for  mass  transfer  during 
bubble  formation  assuming  a  spherical  growth  of  the  bubble 
were  suggested  by  Popovich  et .  al.  [57]  and  are  given  by 
equations  2.22  to  2.25.  The  mass  transfer  flux  at  time  t 
given  by  and  the  total  mass  transfer  from  time  0  to  t 
given  by  are 


=  const  (C®  -  C^)  /Dn  d2  Tf  2/3  t1/6 - (2.22) 

Na  =  const  (Ca-Ca)  /dIT  d2  Tf“2/3  t7/6 - (2.23) 


When  t  =  Tf  the  above  expressions  reduce  to 

NAf  -  const  (Ca-Ca)  /DF  df2  Tf“1/2 - (2.24) 

and 

Na  =  const  (Ca-Ca)  /dF  df2  Tf±7^ - (2.25) 

Inspection  of  equation  2.25  suggests  that  the  total  mass 
transfer  from  the  gas  bubble  upto  the  point  of  detachment 
varies  as  the  square  root  of  Tf.  The  results  of  this 
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study  indicate  that  this  is  not  so.  A  plot  of  total  mass 
transfer  vs.  time  of  formation  is  shown  in  Figure  33  and 
the  slope  of  a  least  square  line  through  the  data  was  found 
to  be  0.685  indicating  that  the  exponent  on  T^.  is  0.685  and 
not  0.50  as  suggested  by  equation  2.25. 

This  difference  can  be  explained  by  the  fact  that 
the  growth  of  the  bubble  in  this  study  was  not  spherical 
whereas  this  is  implicity  assumed  in  equation  2=25.  Popovich 
et .  al.  [57]  confined  their  work  to  liquid  drops  of  approxi¬ 
mately  0.2  to  0.4  mm  diameter  and  spherical  growth  of  the 
drop  for  small  drop  sizes  is  reasonable.  The  equivalent 
bubble  diameters  in  this  investigation,  however,  were  in 
the  range  of  5.2  to  6.3  mm  and  an  assumption  of  spherical 
growth  for  these  sizes  is  not  justified.  In  addition  to 
this  the  variation  of  bubble  surface  area  with  time  was  not 
represented  by  equation  2  ,19  throughout  the  period  of  bubble 
growth  and  this  too  has  been  assumed  in  deriving  equation 
2  .25. 

Bowman  et .  al .  [3]  conducted  some  work  with  CO^- 
water  and  reported  the  variation  of  mass  transfer  coefficient 
with  bubble  frequency  during  the  bubble  rise  period.  Their 
results  indicated  that  mass  transfer  coefficient  increases 
with  bubble  frequency  and  the  results  of  this  study  in¬ 
dicated  that  this  was  true  during  bubble  formation  also, 

A  comparison  of  the  results  of  Bowman  et .  al .  [3]  and  this 
work  is  presented  in  Figure  34. 
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Figure  33:  Plot  of  Total  Mass  Transfer  vs. 
Total  Bubble  Age 
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An  Increase  of  mass  transfer  coefficient  with  in¬ 
creasing  bubble  frequency  can  be  explained  on  the  basis  of 
increased  disturbances  in  the  liquid  phase  at  higher  freq¬ 
uencies.  The  development  of  a  neck  during  bubble  growth 
followed  by  a  bubble  breakup  causes  considerable  turbulence 
in  the  liquid  in  the  region  of  the  orifice  mouth  and  this 
tends  to  increase  the  mass  transfer  coefficient.  For  lower 
frequencies  the  fraction  of  neck  formation  time  to  total 
bubble  growth  time  is  smaller  than  that  for  higher  frequencies. 
In  addition,  there  are  a  larger  number  of  bubble  breakoffs 
from  the  orifice  in  a  given  time  for  higher  frequencies. 

These  statements  are  supported  by  the  results  shown  in 
Table  6,  and  the  increase  of  the  mass  transfer  coefficient 
with  an  increase  of  bubble  frequency  is  therefore  justified. 


Run  No. 

Bubble  Frequency  % 

( 1/sec  ) 

Neck  Formation 
Time 

Mass  Transfer 
Coefficient 
cms/hr . 

Ml 

1.019 

11.2 

161.6 

M2 

5.008 

19-6 

180.0 

M3 

10.002 

24.6 

222 . 4 

M4 

1.027 

11.7 

158.6 

M5 

4.996 

15.6 

149.4 

M6 

9.883 

i — i 

in 

C\J 

214.5 

M7 

0.5014 

6.7 

98.7 

Table 

6 :  Percentage  Neck 

Formation  Time 

in  Relation 

to  Bubble  Frequency  and  Mass  Transfer  Coefficient 
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As  seen  from  Figure  3^,  the  mass  transfer  coefficient 
obtained  during  bubble  formation  is  of  the  same  order  of 
magnitude  as  that  obtained  during  bubble  rise.  The  former 
is,  however,  greater  in  absolute  value.  Calderbank  and 
Patra  [8]  reported  that  the  mass  transfer  coefficient 
during  bubble  formation  is  lower  than  that  during  bubble 
rise  but  their  investigation  was  limited  to  a  stagnant 
liquid  phase.  The  higher  value  of  mass  transfer  coefficient 
during  bubble  formation  obtained  in  this  study  may  be 
attributed  to  a  superimposed  liquid  velocity  on  the  forming 
bubble.  This  is  supported  by  the  fact  that  higher  mass 
transfer  coefficients  during  bubble  rise  due  to  a  moving 
liquid  phase  have  also  been  reported  by  Himmelbau  et .  al. 
[52]. 

Figure  35  presents  the  results  of  the  mass  transfer 
coefficients  during  bubble  formation  obtained  in  this  study 
as  a  comparison  with  the  mass  transfer  coefficients  obtained 
during  free  bubble  rise  for  the  CCi-water  system  by  other 
authors.  The  plot  presents  the  mass  transfer  coefficient 
vs.  the  equivalent  spherical  bubble  diameter  at  bubble 
detachment  from  the  orifice. 

It  is  interesting  to  note  the  wide  disparity  in  the 
results  obtained  by  different  authors  for  the  same  system 
during  bubble  rise.  Bogandy  et .  al.  [69]  and  Leonard  [70] 
indicate  an  increase  in  the  mass  transfer  coefficient  with 
increasing  bubble  diameter  while  Himmelbau  et .  al„  [52], 
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Figure  35:  Plot  of  Mass  Transfer  Coefficient  vs.  Equivalent  Bubble  Diameter 
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Hammerton  and  Garner  [30],  Bowman  et .  al .  [3]  and  Calderbank 
et .  al.  [71]  Indicate  a  decrease  in  mass  transfer  coefficient 
with  increasing  bubble  size  in  the  same  range  of  bubble 
diameters.  The  results  of  Fricke  [72]  together  with  those 
of  Hammerton  and  Garner  [30]  and  Calderbank  et .  al.  [71] 
suggest  that  over  the  complete  range  of  bubble  sizes  the 
curve  of  the  mass  transfer  coefficient  during  bubble  rise 
should  pass  through  a  maximum.  This  fact  has  been  sub¬ 
stantiated  by  Siddique  [51].  The  work  of  Zieminski  et .  al. 
[50]  indicates  a  similar  trend  over  a  relatively  narrow 
range  of  bubble  sizes. 

The  results  obtained  for  mass  transfer  coefficients 
during  bubble  formation  in  this  investigation  indicate 
that-  over  the  range  of  bubble  sizes  studied  the  mass  trans¬ 
fer  coefficient  increases  with  increasing  bubble  diameter. 

The  values  of  mass  transfer  coefficient  obtained  are  of  the 
same  order  of  magnitude,  but  higher  than  those  obtained 
by  other  authors  during  bubble  rise  studies  and  are  in  the 
same  range  as  those  of  Zieminski  et ,  al.  [50]. 
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CHAPTER  7 

CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  study  of  mass  transfer  during  bubble  for¬ 
mation  and  the  analysis  of  the  data  obtained  in  this  in¬ 
vestigation  using  the  CO^-water  system,  the  following 
conclusions  may  be  made. 

(1)  The  total  mass  transfer  during  formation  depends  on 
T^O'685  w^ere  ip  j_s  the  time  of  formation  of  the  bubble. 

This  fact  does  not  correspond  to  the  equations  proposed 
by  Popovich  et .  al .  [57]  wherein  the  mass  transfer  varies 
as  Tf0'5. 

(2)  The  mass  transfer  coefficient  increases  with  increasing 
bubble  frequency  in  the  range  of  frequencies  investigated. 
This  is  supported  by  the  work  of  Bowman  et .  al.  [3]. 

(3)  The  mass  transfer  coefficient  during  bubble  formation 
with  a  superimposed  liquid  velocity  is  higher  than  the  mass 
transfer  coefficient  (i)  during  bubble  rise  in  a  stagnant 
liquid  (ii)  during  bubble  rise  in  a  moving  liquid  (iii) 
during  bubble  formation  in  a  stagnant  liquid.  The  mass 
transfer  coefficients  in  all  the  above  given  cases  are  of 
the  same  order  of  magnitude. 

(4)  With  a  superimposed  liquid  velocity  the  mass  transfer 
coefficient  during  formation  increases  with  increase  of 
the  equivalent  spherical  bubble  diameter  at  detachment 
over  the  range  of  bubble  sizes  examined. 
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(5)  The  bubble  surface  area  during  growth  follows  the 

nl 

relationship  A  =  k^t  up  to  the  point  where  the  bubble 
develops  a  neck.  Subsequent  to  this  the  surface  area 
increases  sharply  with  time  till  bubble  break-off.  The 
average  value  of  n^  found  in  the  investigation  was  0.850. 

The  growth  of  the  bubble  is  therefore  not  spherical. 

(6)  The  volume  of  the  gas  bubble  increases  linearly  with 
time  throughout  the  period  of  formation. 

(7)  The  surface  area  and  volume  of  the  bubble  could  be 
evaluated  at  detachment  with  sufficient  accuracy  using  the 
computer  program  based  on  the  continuous  function  approach. 
The  continuous  function  was  a  modification  of  the  equation 
of  a  cardioid. 

(8)  The  computer  program  based  on  the  numferical  approach 
gives  much  more  flexibility  of  analysis  in  that  it  can  be 
used  to  evaluate  surface  area  and  volume  of  the  bubble  as 
a  function  of  time  during  bubble  formation.  It  can  also 
be  used  to  evaluate  the  area  and  volume  of  the  detached 
bubble . 

The  following  conclusions  could  be  made  from  the 
air-water  runs  conducted  to  characterize  the  system. 

(1)  The  correlation  of  the  frequency  of  bubble  formation 
with  orifice  Reynolds  number  is  independent  of  orifice 
size.  The  frequency  increases  with  gas  flow  rate  up  to  an 
orifice  Reynolds  number  of  about  500  after  which  it  increases 
with  gas  flow  rate  at  a  decreasing  rate.  The  maximum 
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bubble  frequency  is  an  Inverse  function  of  orifice  size. 

(2)  Bubble  volume  increases  with  orifice  diameter  at  a 
given  gas  flow  rate. 

(3)  The  bubble  volume  is  essentially  constant  at  low  gas 
flow  rates  but  increases  during  the  transition  region  with 
increase  of  gas  flow  rate.  The  extent  of  the  constant 
volume  region  and  transition  region  from  a  point  of  view 

of  orifice  Reynold’s  numbers  is  a  function  of  orifice  size. 
Essentially  constant  bubble  volume  is  obtained  below 
Reynolds  numbers  of  30  for  all  orifice  sizes  investigated. 
This  limit  is  higher  for  smaller  orifice  sizes. 

On  the  basis  of  the  results  obtained  during  this 
investigation  the  following  recommendations  for  future 
studies  in  the  field  can  be  made. 

(1)  Additional  data  should  be  obtained  using  the  same 
orifice  size  but  different  liquid  velocities  to  sub¬ 
stantiate  the  effect  of  liquid  velocity  on  mass  transfer 
coefficient.  This  should  be  extrapolated  to  zero  liquid 
velocity  to  compare  the  results  with  mass  transfer  studies 
in  stagnant  liquids  conducted  by  previous  authors. 

(2)  Mass  transfer  data  could  be  obtained  using  different 
size  orifices  to  study  the  effect  of  orifice  size  on  mass 
transfer  coefficient.  This  would  help  in  optimum  system 
design , 

(3)  More  information  is  necessary  on  the  effect  of  small 
amounts  of  surfactants  in  the  system  being  investigated. 
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(4)  Immediately  after  detachment,  the  gas  bubble  undergoes 
rapid  changes  In  shape  constantly  renewing  its  surface  area. 
The  mass  transfer  coefficient  in  this  region  is  likely  to 
be  high.  It  is  suggested  that  attempts  be  made  to  measure 
the  mass  transfer  coefficient  during  this  period. 
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NOMENCLATURE 


a 


a 


1 


A 


A 


e 


b 


1 


c 


C 


c* 


c 

o 


c 


A 


C 


b 

A 


C 


A 


d 


e 


d 


f 


D 


=  Constant  in  equation  of  the  Cardioid  and  subsequent 

equations  presented  in  Chapter  3S  and  Appendix  A. 

=  Lower  limit  of  the  surface  and  volume  integrals  given 

by  equations  3.4  and  3.5. 

=  Surface  area  of  the  bubble  or  drop  (cm  ). 

2 

=  Cathode  surface  area  (cm  ). 

2 

=  Maximum  surface  area  of  bubble  (cm  ). 


Upper  limit  of  the  surface  and  volume  integrals. 
Acoustic  velocity  in  the  gas  (cm/sec). 

Concentration  of  oxygen  in  water  at  a  given  position  in 

the  column.  (6 ms  or  Sm  moles)  or  (SiS.  .P.rm  moles , 

cc  litre 


=  Concentration  of  oxygen  in  the  bulk  water  phase 
/gms  or  gm  moles,  „  /gms  or  gm  moles N 

I  /  o 27  l  jo 

cc  '  litre  ' 

=  Concentration  of  oxygen  on  the  cathode  surface 

,gms  or  gm  moles,  ,gms  or  gm  moles 

cc  or  "  litre  ' 6 

=  Concentration  of  diffusing  solute  TAT. 

/gms  or  gm  moles N  /gms  or  gm  moles N 
cc  j  or  litre 

=  Concentration  of  the  solute  'A'  in  the  bulk  liquid 

,  gms  or  gm  moles N  ^gms  or  gm  moles ^ 

Pnase  *  cc  '  or  1  litre 

=  Concentration  of  solute  ’A’  in  the  liquid  phase  at 
saturation.  (gms  or,  gm  moles.  or  (gms  or  gm  moles 

0  C  _L  1  u  27  0 


=  Bubble  equivalent  spherical  diameter  (cms. ). 

=  Final  equivalent  spherical  diameter  of  the  bubble  at 
formation  (cms.). 

=  Diffusivity  of  the  solute  in  the  liquid  phase  (cm  /sec). 
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=  Function  in  expression  for  derivative  by  Lagrangian 


D., 

lk 


Interpolation.  Given  by  equation  A. 19. 

e„  =  Eccentricity  of  bubble  (  =  7 — - - rj — : — — ,  )  . 

c  J  Maximum  Height 

f(v)  =  Function  describing  bubble  periphery  i.e.  abscissa 
of  a  point  on  the  bubble  periphery  (cms). 
f'(v)  =  Derivative  of  f(v)  with  respect  to  V. 

F  =  Faraday  Number. 

g(v)  =  Function  defined  as  [f(v)  /l+f'(v)^  ]. 

c  c 

G  =  Gas  flow  rate  through  the  orifice  ( - ) . 

s  0  c 

h  =  Spacing  interval  in  Newton  Cotes  integration  formula. 
H  =  Henry's  Law  constant, 
i  =  Diffusion  current  (yA). 


k  = 


k 


1 


Constant  in  equations  for  bubble  periphery  given  as  3.7 
through  3.12,  3.16  and  3.17. 

=  Constant  in  equation  2.19  for  surface  area  of  a  bubble 
as  a  function  of  time. 


K 


^  =  Instantaneous  liquid  phase  mass  transfer  coefficient. 

/emsx  /CmsN 

( - )  or  (r- — )  . 

sec  hr 

kLav=  Average  liquid  phase  mass  transfer  coefficient  during 

bubble  formation.  (^L±L)  or  (Oms_) 

sec  hr 

k’Lav  =  Average  liquid  phase  mass  transfer  coefficient  during 

bubble  rise  period.  ( - )  or  {r — ). 

s  0  c  Jnir3 

kfe  =  Mass  transfer  coefficient  of  oxygen  in  water. 

,-cmsN  /CmsN 

t — -)  or  (r- — )  . 
sec  hr 

K*  =  a  function  of  bubble  eccentricity  used  in  equation  2-15 
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n  =  Exponent  for  equations  representing  models  of  the  bub¬ 
ble  periphery  (equations  3.8,  3-9,  3.16  and  3.17). 
n,  =  Exponent  of  time  in  equation  2,19  for  variation  of 
bubble  area  with  time. 

n^  =  Exponent  of  dimensionless  time  in  surface  area  equations 

2.26. 


n*  =  Exponent  of  polar  coordinate  ' r '  in  equation  6.2. 
n  =  Number  of  electrons  taken  up  during  the  electrochemical 
reaction . 


N  =  Normality  of  Ba(OH)2  solution. 

N-^  =  Normality  of  HC£  solution. 

Ncq  =  Normality  of  carbon-dioxide  in  solution 


N*  =  Mass  transfer  of  oxygen  in  water  to  the  cathode, 
(gms  or  gm  moles). 

N^  =  Total  mass  transfer  in  time  ’tT  (gms  or  gm  moles). 

N’  -  Mass  Transfer  rate  at  time  Tt’  (^^ — or  gm  moles ^ 

A  sec 

N  =  Capacitance  Number  (  =  ^ &-{ ~ P P- i-V  £ )  e 

n6"pd  c“ 

N^  =  Total  number  of  data  points  in  a  given  population 
size  of  the  population. 

N^  =  Representative  sample  size  for  a  given  population. 

pfl  -  Partial  pressure  of  'A’  in  the  gas  phase  (mm.Hg). 

dp  U 

Pe  =  Pec let  Number  of  bubble  or  drop  (  =  — )  . 

P^(x)  -  i^^  Lagrangian  polynomial,  in  x. 

P.(x)  =  Polynomial  of  degree  j ( j =0 , 1 , . . . ,m) . 

J 

r  =  Polar  coordinate  of  bubble  periphery. 


i 
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r  =  Radius  of  the  column  (cms). 
o 

r-^  =  Any  radial  position  in  the  column  from  the  centre  line 
( cms ) . 

rg  =  Radius  of  a  sphere  of  volume  equivalent  to  the  axi- 

symmetric  body  of  revolution  (cms). 

Re  =  Orifice  Reynolds  Number  (  =  — — )  . 

°  2 

S  =  Surface  area  of  revolution  of  a  body  (cm  ). 

Sc  =  Schmidt  Number  (  =  - =r)  . 

P  A  D 

Sh  =  Sherwood  Number  (  =  L  e)„ 

D 


t  =  Time  ( sec  )  . 

t  =  Time  taken  by  bubble  to  rise  a  distance  equal  to  its 
equivalent  spherical  diameter  at  its  vertical  rise 
velocity .  ( sec ) . 

T^  =  Total  time  of  formation  of  bubble  (sec). 

=  Average  liquid  velocity  in  the  column  (cms/sec), 

Ug  =  Velocity  of  the  gas  at  the  orifice  (cms/sec) . 

U  =  Free  rising  or  falling  velocity  of  bubble  or  drop  (cms/sec) 
v  =  Vertical  height  of  the  bubble  or  the  ordinate  of  a  point- 
on  the  bubble  periphery  (cms.). 

V  =  Volume  of  revolution  of  a  body  or  volume  of  a  bubble  (cc). 
V^  =  Volume  of  HC£  used  to  neutralize  a  given  volume  of 
Ba(OH)p  solution  in  a  blank  titration  (cc), 

V^  =  Volume  of  HC£  used  at  the  end  of  the  first  stage  of 
the  back  titration  (cc), 

V*  =  Volume  of  a  detached  bubble  estimated  from  photographs 
(cc). 
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Vc  =  Chamber  Volume  (cc). 
x  =  Distance  from  the  electrode  (cms). 
x  =  Mean  of  the  sample . 

x^  =  Value  of  the  individual  data  point  in  a  given  population. 
xg  =  Mean  of  the  sample. 

X*  =  Characteristic  length  for  platinum  probe  1  (mm) . 

=  Moles  of  solute  'A'  per  mole  of  CO^-Water  Solution. 

y.  =  Value  of  a  function  of  Tx'  at  x=x0. 

J  J 

Y*  =  Characteristic  length  for  platinum  probe  2  (mm) . 

Greek  Symbols 

a o  =  Factor  in  Forsythe  polynomial  given  by  equation  A.l6. 

J 

3c  =  Factor  in  Forsythe  polynomial  given  by  equation  A . 17 . 

J 

6  =  Orifice  diameter  (cms). 

0  =  Polar  coordinate  for  bubble  periphery  and  an  axi- 

symmetrie  body  of  revolution, 
li  =  Standard  deviation  of  the  individual  data  point  from 
the  mean  of  the  population. 

u  =  Viscosity  of  the  continuous  phase  (^ -  )  . 

p  =  Density  of  the  continuous  phase  (gms/cc). 

=  Density  of  the  dispersed  phase  (gms/cc). 
v  =  Frequency  of  bubble  formation  (sec  ^). 
nm  =  Maximum  frequency  of  bubble  formation  attained  in  the 

/  -In 

constant  frequency  region  (sec  ). 
c^  =  Cut  off  limits  or  sigma  limits  on  a  given  distribution 


. 
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of  Individual  data  points  for  a  given  level  of  confidence, 
(p  (x)  =  Function  for  a  polynomial  fit. 

Q  (x)  =  Function  given  by  equation  A. 21. 
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APPENDICES 


Appendix  A 

Evaluation  of  Bubble  Surface  and  Volume  for  Program  I 
With  reference  to  Figure  2 


b-L  _ _ 

S  =  2n  /  f  (v)  /  1+f  ’  (v)  dv - 

al 

r  =  a(k+Cosn0)  - 

f(v)  =  rSinO  =  a (k+Cosn0 ) Sin0 

f '  (v)  =  ^££0  - 

d0  dv 

v  =  r  Cos  0  =  a(k+Cosn0)  Cos0 
dv  =  aSin0 (k+ (n+1 ) Cosn0 )d0  -- 


A-i 


- - - A- 2 

akSin0+aSin0Cosn0  - A-3 

- A-4 

a.kCos0+aCosn+10 - A-5 

- ___ - A-6 


d0 

dv 


_ 1 _ 

aSin0 (k+ (n+1 ) Cosn0 ) 


also 


df  (v) 
d0 


ak  Cos0 


2  n— 1  n+1 

anSin  0Cosx  0  +  aCos  0 


A-7 


A-8 


substituting  in  A-4  we  get 


I  ■ 

■ 

\ 

\ 


A-2 


2  n— T  n  + 1 

ft(v)  =  nSin  0  Cos  0  -  k  Cos9  -Cos  0 _ A_9 

SinG  (k+(n+l )Cosn0) 

Substituting  for  f(v)  and  f'(v)  and  dv  in  A-l  we  get 
S  =  211  /  -a2Sin0  (k+Cosn0 )  [  Sin20  (k+  (n+1 )  Cosn0 )  2  + 

n/2 

(nSin20Cosn“10-kCos0-Cosn+10)2]1/2d0  - A-10 

The  volume  is  evaluated  as  shown  below: 

b  2 

V  =  n  /  f(v)  dv - - - - - A-ll 

a 

substituting  for  f(v)  from  A-3  and  for  dv  from  A-6  we  get 

V  =  n  f  -  a3Sin30(k+Cosn0)2(k+(n+l)Cosn0)d0  - A-12 

n/2 

Equations  Used  for  development  of  Program  II 

Forsythe  Orthogonal  Polynomial 

The  Forsythe  Orthogonal  Polynomial  combined  with 
the  least  square  criterion  is  used  in  this  investigation  to 
smcothen  the  raw  data  obtained  from  photographic  films. 

The  maximum  order  permitted  is  12.  The  standard  deviation 

-4 

is  of  the  order  of  10  .  The  program  terminates  when  either 

of  the  two  limits  specified  above  is  first  achieved. 
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The  equations  are  given  as  follows: 


j=m 

<P  (x)  =  l  b  P  (x)  - A-13 

j=o  J  J 


Where  0m(x)  is  the  representation  of  the  descrete 

data  points  in  the  form  of  a  polynomial  fit  where  P.(x)  is 

J 

a  polynomial  of  degree  j .  The  Forsythe  orthogonal  poly¬ 
nomials  are  given  by 


P_l(x) 

P0(x)  = 

P  (x)  = 
1 


t 

I 

I 


p . 
J 


+1 


(x) 


0 

1 

(x-a1)Po(x)-3QP_1(x) 


=  (x-a  )P  (x)  -  3  P  (x) 

tJ  - J  J 


A-14 


where  3Q~0.  The  values  of  a  and  3  are  chosen  in  such  a  way 
that 


In  general 


A-15 


A-16 


and 


- 


' 


i? 


Evaluation  of  the  Derivative  at  descrete  points  Using 


Lagrangian  Interpolation 

The  equation  used  to  evaluate  the  derivative 

dy  at  a  point  (x  ,y,  )  is  given  by 
dx  *  K 


?’L  =  {  l  ll  (n  (x  -x  ))+Y  l  - A-18 

Xk  i=0  D  3=0  J  Xj=0  Xk  xj 

i+k  K  j  =f=k  j  =t=k 

Where  is  given  by 

Dik =  (xk-xi)pi(xi) - A-19 

and 

n 

P.  (x.  )  =  II  (x.-x.)  - A-20 

11  3=0  1  J 

j+i 

Pi (xi)  is  the  Lagrangian  Polynomial  at  x^ . 

Lagrangian  Interpolation  Formula 


This  is  implicitly  stated  in  the  previous  section 


a a1  _  sdnb  j_s>Tt)ioasb  3£  :  iJ-rilL  -  — ■ ~  — '  — ^ 


. 


%  • 


A-5 


n 


$..(x)  =  l  pi(x)  Y4  - A-21 

n  i=o  qrrp  1 

n 

where  P.(x)  =  n  (x-x.)  - A-22 

1  j =0  J 
J+i 


Newton  Cotes  4th  order  integration  Formula  with  1  m*  subintervals 

The  integration  formula  for  five  evaluation  points 
xq}  x13  .  ..,  x is  given  by 


x  i 


/  f (x)dx  = 


7 

(7y0+32y1+l2y2+32y3+7y1))-^5-f[6](e) 


- A-23 


■o 


Where  Th’  is  the  interval  spacing.  This  formula 
may  be  applied  repeatedly  to  cover  more  than  5  points,  e.g. 


x, 


•A.  9  -9 

f  f(x)dx  =  J  f(x)dx  +  /  f(x)dx  and  so  on. 


x 


0 


x 


0 


x , 


In  general,  if  this  is  applied  ’m’  times  in  succession  the 
equation  is  given  by 


‘4m 


/  f(x)dx  =  ||  [7yQ  +  32(y1+y3+y5+.  .  .yi|m_1)  + 


x 


12(y2+y6+.  .  •y„m_2)+l1*(y1,+  .  •  •7l)m_4^+7yi)m] 


A-24 


V 


. 


+  (;  V  i)  ^V"  ■  ’  ^  1 
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Appendix  B 

Program  I 

The  program  is  entitled  PROG  1  and  consists  of  the 
mainline  program  five  subroutines  and  two  function  sub¬ 
programs.  The  subroutines  are  entitled  PRELIM,  LSQFIT, 
GAUSS,  STD  and  GAQ10;  the  function  subprograms  are  called 
FCT1  and  FCT2 . 


Input  to  the  program  is  done  completely  in  the 
mainline  and  includes  the  following  sequence  of  cards. 

(1) NN,  Ml  where  NN  is  the  total  number  of  data  points  minus 
one  and  Ml  Is  the  number  of  functions  of  x  involved  in  the 
least  square  fit  routine  LSQFIT. 

(2)  XL,  SV : -  The  lower  and  upper  limits  of  integration 
respectively.  These  are  used  in  the  routine  GAQ10 .  For 
this  investigation  XL=n/2  and  SV=0. 

(3)  NUMRUN  which  indicates  the  total  number  of  separate 
runs  being  analysed. 

(4)  LI:-  The  Total  number  of  bubbles  being  analysed  in 
any  given  run. 

(5)  G,U,SP,D  where  G  is  the  gas  flow  rate  in  cc/sec,  U  is 
the  liquid  velocity  in  cms/sec,  SP  is  the  film  speed  in 
frames/sec  and  D  is  the  outside  orifice  diameter  in  cms . 

(6)  M  which  indicates  the  frame  number  of  the  bubble  being 
analysed.  This  is  indicated  on  the  film  projector. 

(7)  0R1  (I),  1=1,4  :-  These  are  four  numbers  representing 
the  x  and  y  coordinates  of  the  two  edges  of  the  orifice  to 
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give  the  orifice  size  on  the  projected  film.  These  numbers 
are  obtained  from  the  photographic  film  of  the  bubble  with 
the  aid  of  the  digitizer.  The  Y  coordinate  should  precede 
the  X  coordinate  in  all  cases. 

(8)  Data  points  from  the  bubble  periphery:-  These  are 
obtained  with  the  digitizer.  A  format  of  (1615)  is  specified 
indicating  a  maximum  of  16  numbers  or  8  points  on  the  bubble 
periphery  per  card.  The  sequence  in  which  the  bubble 
periphery  should  be  traversed  is  shown  in  Figure  36 .  The 
total  number  of  cards  used  up  with  this  data  depends  on 

the  total  number  of  points  taken  on  the  bubble  periphery. 

(9)  NAPEXV,  NAPEXF : -  These  are  the  Y  and  X  coordinates 
respectively  of  the  apex  of  the  bubble  i.e.  corresponding  to 
point  1?  in  Figure  36.  This  is  read  separately  from  the 
rest  of  the  data  so  that  the  same  data  deck  could  be  used 
for  Program  II . 

For  each  bubble  belonging  to  the  same  run,  cards 
corresponding  to  sections  6  through  9  above  must  be  in¬ 
cluded.  For  each  run, cards  corresponding  to  sections  4 
through  9  should  be  included. 

The  operations  performed  by  different  parts  of  the 
program  are  as  follows : 

Mainline  Program:-  (i)  Inputs  all  the  data  mentioned  above. 

(ii)  Transforms  the  data  of  the  bubble 
periphery  so  that  the  bubble  is  symmetrical  above  the 
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Figure  36:  Sequence  of  Reading  the 
Bubble  Periphery 


'■ 
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central  axis . 

(iii)  Converts  the  digitized  data  to 
actual  bubble  size  by  comparison  with  the  actual  orifice 
diameter  which  is  supplied  as  data. 

Subroutine  PRELIM:-  (i)  Determines  the  values  of  ’a’  and 

' k ’  which  are  parameters  of  the  equation  r=a(k+Cos0n) 

2? 

(ii)  Determines  ln( — k)  and  ln(Cos0) 

a 

for  each  data  point  on  the  bubble  and  stores  them  as  vectors 
Y  and  X  respectively. 

Subroutine  LSQFIT:-  (i)  Determines  the  set  of  simultaneous 
equations  through  the  least  square  formulation 
Subroutine  GAUSS:-  (i)  Solves  the  set  of  simultaneous 
equations  determined  in  LSQFIT  by  Gaussian  Elimination 
Method  to  get  a  value  of  the  parameter  'n'. 

Subroutine  STD:-  (i)  Determines  the  standard  deviation  of 
the  least  square  fit . 

Subroutine  GAQ10:-  (i)  Integrates  the  functions  for  bubble 

volume  and  surface  derived  in  Appendix  "A'  by  a  ten  point 
Gaussian  Quadrature  formula. 

Function  Subprogram  FCT1:-  (i)  Defines  the  function  for 
bubble  surface. 

Function  Subprogram  FCT2:-  (i)  Defines  the  function  for 
bubble  volume. 

The  output  from  the  program  includes  a  listing  of 
the  following: 

(i)  The  run  number  being  considered. 

(ii)  The  number  of  the  bubble  being  analysed. 
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(iii)  The  frame  number  the  bubble  corresponds  to 

(iv)  The  gas  and  liquid  flow  rates,  the  film  speed  and 
the  orifice  outside  diameter  . 

(v)  The  standard  deviation  of  the  least  square  fit  to 
determine  the  parameter  fn’. 

(vi)  The  bubble  surface  in  cm  . 

(vii)  The  bubble  volume  in  cm  . 

Program  II 

The  program  is  entitled  "FORMING"  and  is  used  to 
evaluate  the  bubble  surface  and  volume  using  a  numerical 
procedure.  It  consists  of  the  mainline  program  and  four 
subroutines  which  are  entitled  FORSYT,  LADERI,  LAGINT  and 
NCOTES . 

The  input  to  the  program  is  done  in  the  mainline 
program  only  and  consists  of  the  following  sequence  of 
cards : 

(1)  ERROR:-  This  is  a  specification  of  the  upper  bound 

on  the  variance  in  the  Forsythe  Orthogonal 
Polynomial  fit  used  to  smoothen  the  raw  data. 

(2)  NUMRUN:-  Indicates  the  total  number  of  runs  analysed. 

(3)  G, U,SP,D, TIMER, (C(I) ,1=1,4) The  variables  are 

defined  as  (i)  G  =  Gas  flow  rate  in  cc/sec 

(ii)  U  =  Liquid  velocity  in  cm/sec 

(iii)  SP=  Film  speed  in  frames/sec 

(iv)  D  =  Outside  orifice  diameter 
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(v)  TIMER  =  The  cycle  timer  fre¬ 
quency  in  cycles/sec.  This 
variable  is  not  used  in  the 
program  but  could  be  output 
if  required  and  is  useful 
to  check  on  the  film  speed 
SP. 

(vi)  C(I),I=1,4  =  These  are  four 
coefficients  used  in  the 
fourth  order  Newton  Cotes 
integration  formula.  The 
coefficients  are  7.0,  32.0, 
12.0  and  14.0. 


(4) 

Ll:- 

The 

total 

number 

of 

bubbles  to 

be  analysed. 

(5) 

Ml:- 

The 

total 

number 

of 

frames  per 

bubble.  This  is 

used  for  cases  where  surface  area  is  desired  as 
a  function  of  time  for  any  one  given  bubble 
during  the  period  of  formation. 

(6)  N,N1:-  N  is  the  total  number  of  points  on  the  bubble 

periphery  excluding  the  apex  of  the  bubble  e.g. 
with  reference  to  Figure  36,  N=l6.  N1  is  the 
number  of  intervals  the  total  bubble  height  has 
been  divided  into  upto  the  value  of  N.  i.e. 

-  1].  With  reference  to  Figure  36,  Nl=7 . 

(7)  L2(M):-  The  frame  number  of  the  frame  being  analysed. 
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(8)  The  data  subsequent  to  section  (7)  above  includes  the 
same  deck  as  indicated  in  sections  (7),  (8)  and  (9) 
of  Program  I.  Even  though  both  coordinates  of  the 
apex  point  are  supplied3  the  program  reads  only  the  Y 
coordinate  and  generates  the  X  coordinate  itself. 

If  the  program  is  being  used  for  analysis  of  bubble 
formation,  then  M1>1.  Cards  corresponding  to  sections  6, 

7  and  8  above  must  therefore  be  supplied  for  each  frame 
analysed.  If  more  than  one  bubble  is  being  analysed  in  a 
given  run,  cards  corresponding  to  sections  5  through  8  must 
be  included  for  each  bubble. 

The  complete  set  of  data  excluding  the  cards  1  and 
2  for  the  ERROR  and  NUMRUN  respectively,  must  be  included 
for  each  run  being  analysed  if  the  total  number  of  runs  is 
greater  than  1. 

The  function  of  the  subroutines  used  in  the  program 
is  as  follows : 

Main  Line  Program:-  The  operation  of  the  mainline  in  this 
program  is  identical  to  that  in  Program  I. 

Subroutine  FORSYT : -  Smoothens  the  data  obtained  from 
photographs  by  using  the  Forsythe  Orthogonal  Polynomials. 
Subrout ine  LADERI : -  Evaluates  the  derivative  f  ’  (v)  at 
each  of  the  smoothed  data  points.  The  Lagrangian  inter¬ 
polation  formula  is  used. 

Subroutine  LAGINT:-  Generates  equispaced  data  for  f(v)  and 

/ - 

for  f(v)  /1+f 1 (v)  using  a  third  order  Lagrangian  Inter- 
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polation  formula. 

2 

Subroutine  NCQTES:-  Integrates  the  functions  f(v)  and 
f  (v)^L+f  f (v)2  with  respect  to  'v'  by  a  fourth  order'K' 
panel  Newton  Cotes  formula.  The  number  of  panels  K  are 
determined  by  the  total  number  of  points  on  the  bubble 
periphery.  This  number  is  generated  in  the  main  line 
program. 


If  detached  bubbles  are  being  analysed  the  output 
from  the  program  includes  the  following: 

(i)  The  run  number  being  considered 

(ii)  The  number  of  bubble  being  considered 

(iii)  The  gas  flow  rate  in  cc/sec  and  the  liquid  velocity 
in  cms/sec 

(iv)  The  frame  number  and  the  cumulative  time  in  seconds 

Frame  Number 


determined  by  the  ratio 


Film  Speed 


(v) 


The  bubble  surface  in  cm  . 

,3 


(vi)  The  bubble  volume  in  cm  . 

(vii)  The  number  of  panels  used  in  the  Newton  Cotes  inte¬ 
gration  formula. 

If  different  frames  from  a  single  bubble  are 
analysed  to  obtain  surface  area  and  volume  as  a  function 
of  time  of  formation,  the  output  in  sections  (iv)  through 
(vii)  appears  as  a  matrix  in  which  each  row  corresponds  to  a 


frame  analysed. 
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Program  III 


The  program  is  called  "SAMPLE"  and  is  used  to 
estimate  the  sample  size  form  a  given  population  of 
bubbles  so  that  the  bubble  area  and  volume  evaluated  from 
the  sample  is  representative  of  the  total  population. 

The  program  input  consists  of  the  following 
sequence  of  cards. 

(1)  M1,K1:-  These  are  identification  numbers  for  the 

run  being  investigated. 

(2)  N,UP:-  N  is  the  total  population  size  and  UP  is  the 

sigma  limit  for  a  given  level  of  confidence. 

(3)  X(I),I=1,N:-  These  are  the  N  actual  data  points  from 

which  the  sample  size  is  to  be  evaluated. 

(4)  LI:-  This  indicates  tht  total  number  of  runs  being 

analysed.  This  card  is  required  only  if  L1>1 
and  is  included  after  the  data  of  the  first  run. 
The  output  from  the  program  includes  the  following 

(i)  The  identification  numbers  of  the  orifice,  liquid 
velocity  and  the  number  of  the  run. 

(ii)  The  individual  data  points  X  (I) ,1=1, N. 

(iii)  The  mean  value  of  the  variable  ’X’  from  the  total 
population . 

(iv)  The  standard  deviation  from  the  mean  for  the  data. 
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It  is  defined  as 


n  p 

STDEV  =  /  l  (X(I)-X) 

1=1 

where  X  is  the  mean  of 

i — 1 

i 

the  population. 

(v)  The  sample  size  required. 


The  source  listing  of  the 

three  programs  is 

included . 
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-  PROG  1  - 


-  PROG  1  - 

REAL  FV(50 ) »V(50 ) »R(50 )  » COSTA (50)  ,X(50)  ,Y(50),Q(25»50) 
1  » QT ( 50  »  2  5  )  » Cl ( 25  )  » D2 ( 2  5 , 2  5  )  »  P  (  2  5  )  » D AUG (25*25),  DEV  (  50  )  , 
20R (4 )  » $D  >  AC  »  ACP 1 *A*C *S »VOL » I  NT  1 , I  NT 2  » G  *  U » SP » D » D 1 * XL » 
3XU»X1 (50,25  )  ,  Y 1 ( 5  0  ) 

INTEGER  F  V 1 ( 50 )  *  VI ( 50 )  ,0R 1 ( 4 )  »N>M»NM1 » L  » L 1 » M 1  » M  2 ♦ IM1 » 

1  I  PI »M3 »KP1 *N1 
EXTERNAL  FCT1,FCT2 

C  THE  PROGRAM  EVALUATES  THE  SURFACE  AND  THE  VOLUME  OF  GAS 
C  BUBBLES  THAT  HAVE  BEEN  JUST  DETACHED  FROM  THE  ORIFICE. 

C  THE  TOTAL  NUMBER  OF  BUBBLES  AND  THE  GAS  AND  LIQUID  FLOW 

C  RATES  AND  THE  ORIFICE  SIZE  WILL  BE  READ  IN  AND  PRINTED 

C  OUT. THE  BUBBLES  HAVE  BEEN  PHOTOGRAPHED  BY  A  HIGH  SPEED 
C  CAMERA  AND  DIGITAL  DATA  OF  THE  BUBBLE  PERIPHERY  IS  TAKEN 
C  BY  A  DIGITISING  MACHINE. THE  FOLLOWING  EQUATION  HAS  BEEN 
C  FOUND  TO  FIT  THE  BUBBLE  PER  I FERY  CLOSELY  AND  WILL  BE  USED 
C  IN  THIS  PROGRAM  AS  A  BAS  IS • R= ( A ) # ( K+COS ( THETA ) **P )  WHERE 
CP  IS  AN  EXPONENT  TO  BE  DETERMINED. 

RE AD ( 5  ,94 ) NN  »  Ml 
94  FORMAT ( I  3  *  12 ) 

C  NN  IS  THE  TOTAL  NUMBER  OF  POINTS  MINUS  ONE. 

C  Ml  IS  THE  NUMBER  OF  FUNCTIONS  OF  X  INVOLVED  IN  THE  FIT. 

C  I.E.  IF  Y= A+A 1U 1 ( X )  THEN  M 1 = 1 , I F Y- A  +  A 1U 1 ( X ) + A2U2 ( X )  THEN 
C  M 1  =  2  ETC. 

C  THE  BOTTOM  APEX  POINT  IS  OMITTED  SO  THAT  THE 
C  DATA  MAY  BE  MADE  CONFORMABLE  TO  A  NUMERICAL  PROGRAM 
C  WRITTEN  FOR  A  SIMILAR  PURPOSE . N  WILL  BE  SUBSEQUENTLY 
C  DEFINED  AS  NN+1. 


RE AD (5,96)  XL,  XU 
96  FORMAT (2F12. 8 ) 

READ ( 5,91) NUMRUN 
91  FORMAT (14) 

N  R  U  N  =  1 
90  L=  1 

READ( 5 ,99 ) LI 
99  FORMAT (14) 

READ( 5 ,98  )G,U,SP  ,D 
98  FORMAT (F10.8,F10.5,F9.2,F9.5) 

100  READ ( 5,93  )M 
93  FORMAT (15) 

READ ( 5,101)  ( OR1 (  I  )  ,  1  =  1 ,4) 

101  FORMAT (415) 

READ ( 5 ,95 ) ( VI ( I ) ,FV1 ( I ) » 1=1 ,NN ) 
95  FORMAT ( 16 15 ) 

READ ( 5,9  2 ) NAPEX V  ,  NAPEXF 
92  FORMAT (215) 
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N=NN+1 

VI  (  N ) =NAPEXV 
F  V 1 ( N ) =NAPEXF 

C  THE  ABOVE  STATEMENTS  SUPPLY  ALL  THE  DAT  A • N  =  NO • OF  DATA 
C  POINTS  M=FRAME  NO.G=GAS  FLOW  RATE  U=LIQUID  FLOW  RATE 
C  SP  =  FRAME  SPEED  FV » V=  DATA  COORDINATES  AND  D  =  ORIFICE  SIZE. 

C  THE  DATA  WILL  NOW  BE  TRANSFORMED  THROUGH  AN  AXIS  SHIFT 
C  AND  CORRECTED  FOR  MAGNIFICATION  AND  DISTORSION. 

C  FLOATING  DIGITAL  DATA 

DO  97  I  =  1  »  N 

FV (  I ) =FLOAT ( FV1 (  I  )  ) 

V(  I )=FLOAT( VI ( I  )  ) 

97  CONTINUE 

DO  50  1=1*4 

OR ( I ) =  FLOAT ( OR  1 (  I  )  ) 

50  CONTINUE 

D 1 =OR ( 4 ) -OR ( 2 ) 

102  DO  103  1  =  1, N 

FV (  I ) =  FV (  I )*D/D1 
V ( I ) = V ( I ) *D/Dl*3.48/3.42 

103  CONTINUE 

C  EVALUATING  THE  MEAN  *  Y •  COORDINATE  FOR  TWO  DATA  POINTS  IN 
C  THE  SAME  HORIZONTAL  LEVEL . CONS  I DER I NG  THE  TOTAL  NUMBER  OF 
C  POINTS  TO  BE  ALWAYS  ODD. 


N  M  1  =  N  “  1 

DO  104  I  =  2  *  NM 1 »  2 
I M  1  =  I  - 1 

V ( I ) = ( V( I )+V( IM1  )  )  /2.Q 
V ( I  Ml ) = V (  I  ) 

104  CONTINUE 

C  TRANSFORMING  COORDINATES  BY  AXIS  SHIFT 
SH I FT  =  V ( 1  ) 

SHIFT1=(FV(2)+FV(1)  )/2.0 
DO  10  5  1  =  1  » N 
FV (  I ) =  FV (  I  ) -SH I F  T 1 
V (  I ) = V ( I  ) -SH I  FT 

105  CONTINUE 

C  THE  ABOVE  TRANSFORMATION  GIVES  A  POSITIVE  AND  A  NEGATIVE 
C  VALUE  OF  FV  FOR  EACH  VALUE  OF  V » THE  VALUES  MUST  BE  OF  THE 
C  SAME  ABSOLUTE  MAGNITUDE. 

DO  106  I  =  2  * N M 1 »  2 
I M 1  =  I  “  1 
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-  PROG  1  - 

SUM= ( A3S ( FV (  I  )  ) + ABS ( FV ( I  Ml  )  )  )  / 2  •  0 
FV (  I ) =  +  SUM 
FV ( IM1 )=-SUM 

106  CONTINUE 

C  SINCE  COS ( THE TA )  IS  POSITIVE  IN  THE  FIRST  AND  FOURTH 
C  GUADRANTS  WE  NEED  TO  CONSIDER  ONLY  HALF  THE  GIVEN  POINTS. 

CALL  PR EL  I M  C FV> V  »R  *  COSTA »Y»X*AC»ACP1»A»C*N»NM1) 

Nl= (NM1/2 ) 

I NDEX=2 

DO  120  I  =  1  *  N 1 

XI ( I *1 )=X ( I+INDEX) 

Y 1 ( I ) = Y ( I+INDEX ) 

I NDEX= INDEX+1 

120  CONTINUE 

CALL  LSQF IT ( XI , Yl »Q  »QT  *C1 » D2  * P * N 1  * Ml »M2  ) 

CALL  GAUSS  (XI »  Y 1 » Q  >  QT  » C 1 >  D2  »  P ♦ DAUG ♦ N 1  *  M 1  *  M2 ) 

CALL  STD(P»Y1»QT  »DEV*SD»N1*M2) 

CALL  GAQ10 (XL»XU»FCT1* INTI  * A  * C  * P  ) 

C  THE  REQUIRED  FUNCTION  HAS  NOW  BEEN  INTEGRATED  AND  THE 
C  SURFACE  MAY  BE  EASILY  EVALUATED. 

S-2.0*3.141S926*INT1 

C  AREA  OF  THE  CIRCULAR  FLAT  REGION  AT  THE  BOTTOM  IS  GIVEN  BY 

Sl=3. 1415926*( AC**2 ) 

S=S+S1 

C  EVALUATION  OF  VOLUME  REQUIRES  ANOTHER  INTEGRATION  OF  A 
C  DIFFERENT  FUNCTION. 

CALL  GAQ10 (XL»XU»FCT2  » INT2»A»C»P) 

VCL=3 . 1415926*1 NT2 
IF ( L.GT. 1  )  GO  TO  121 
WRITE (6  » 119  JNRUN 

119  FORMAT ( • 1 • »20X» ' RESULTS  FOR  RUN  NO.  P='»I3) 

121  WRITE (6* 107 )L 

107  FORMAT (' 1 ' »20X »' DATA  AMD  RESULT  FOR  BUBBLE  N0.='»I4) 
WRITE (6 » 108  ) 

108  FORMAT ( 'O' »10X» 'G  CC/SEC ' , 5X , ' U  CM/ SEC ' » 5 X , 1 SP  F/SEC', 
1 5X  » ' D  CM '  » 7X  » ' M 1  ) 

WRITE (6  » 109 )G*U*SP*D*M 

109  FORMAT (  '  '  » 7 X » F 1 0 . 5 » 3X » F 1 0 . 5 » 4X » F9 . 2 > 3 X » F9 . 5 » 2X »  I  6 ) 

WRITE (6. 115 )SD 

115  FORMAT('  '» 10X ♦' STANDARD  DEVIATION  OF  FIT  =  '.E20.8) 
WRITE (6 » 116 )S* VOL 

116  FORMAT('  '» 10X, ' BUBBLE  SURFACE  IN  SO. CMS  =  ' *E16.8»/ 
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-  PROG  1  - 

1 10X  »  '  3UBB  LE  VOLUME  IN  CC  =  SE16.8) 
IF ( L.GE.L1  )GC  TO  117 
L  =  L+1 
GO  TO  100 

117  I F ( NRUN • GE • NUMRUN )  GO  TO  122 
NRUN=NRUN+1 
GO  TO  90 
122  STOP 

END 


' 
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SUBROUTINE  PRELIM 

SUBROUTINE  PRELIM ( FV » V » R *  COST A  * Y ♦ X » AC » ACP 1 » A  » C ♦ N »NM1 ) 
REAL  FV(50)  »  V ( 50 )  »R(50)  > COSTA (50).Y(50)»X(5Q) , AC  > ACPI » 
1 A  » C 

INTEGER  Ni.NMl 

C  THIS  SUBROUTINE  EVALUATES  THE  PARAMETERS  A  AND  C  AND  THEN 
C  COMPUTES  R  AND  COS ( THETA ) . I T  THEN  REPLACES  THE  VALUES  OF 
C  LOG ( R/ A-K )  BY  Y  AND  LOG  COS (THETA )  BY  X. 

ACP 1= V ( N ) -V ( 1 ) 

AC  =  ABS ( F V ( 1 )  ) 

A= ACP 1 -AC 
C= AC/A 

C  R  IS  THE  POSITION  VECTOR  OF  EVERY  POINT  AND  COSTA  IS  THE 
C  COSINE  OF  THE  ANGLE  OF  THE  POSITION  VECTOR  WITH  THE  AXIS. 

DO  10  I  =  3  *  N  ♦  2 

R (  I ) =SQRT (  ( ABS ( F V ( I  )  )  )#*2+(V( I >**2)  > 

COSTA ( I ) = V (  I  )  /R (  I  ) 

Y (  I ) =ALOG (  (R(  I ) /A)-(C)  ) 

X(  I )=ALOG( COSTA (  I  )  ) 

10  CONTINUE 

11  RETURN 
END 
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SUBROUTINE  LSGFIT 

SUBROUTINE  LSQFIT  ( X > Y > Q * QT ♦ C * D » A » N * M » Ml ) 

REAL  X(50»25) ,Y(50) >0(25, 50), QT (50*25) >C(25) *D(25»25) , 
1A ( 25 ) 

INTEGER  N*M*M1» I *<»J* I  Ml* I  PI >M2*KP1 
200  DO  20 1  J  =  1 » N 
Q  (  1 »  J  )  =  1  •  0 
201  CONTINUE 
Ml =  M+ 1 

DO  20  2  I  =  2  * M 1 
IM1= 1-1 
DO  203  J=  1  » N 
Q{ I » J )=X ( j, IM1 ) 

203  CONTINUE 
202  CONTINUE 

C  THE  MATRIX  Q(I»j)  MAY  BE  TERMED  AS  A  COEFFICIENT 
C  MATRIX  IN  WHICH  THE  FIRST  ROW  IS  1  AND  THE  SUBSEQUENT  ROWS 
C  HAVE  Xll ,X21 »X31 ,ETC  THIS  MATRIX  HAS  BEEN  TRANSPOSED  FROM 
C  THE  ORIGINAL  SET  OF  EQUAT I ONS . SEE  DIGITAL  COMPUTATION  FOR 
C  CHEMICAL  ENGINEERS  BY  LAPIDUS. 

C  NOW  WE  SHALL  TRANSPOSE  Q(I*j)  AGAIN  SO  THAT  WE  MAY  EXPRESS 
C  THE  LEAST  SQUARE  FORMULATION. 

204  DO  205  I  =  1  ,  N 
DO  206  J  = 1  * M 1 
QT (  I » J)=Q( J*  I  ) 

206  CONTINUE 

205  CONTINUE 

C  A  COLUMN  VECTOR  C ( I )  IS  DETERMINED  WHICH  CORRESPONDS 
C  TO  SIGMA(1  TO  N)  OF  EACH  COLUMN  OF  Q ( I  * J ) *Y ( J  )  *TH I S  GIVES 
C  N  COLUMN  VECTORS  WHICH  ON  ADDITION  OF  CORRESPONDING 

C  ELEMENTS  YIELDS  C ( I ) • 

DO  207  1  =  1  *  Ml 

C (  I ) =0 .0 

DO  208  J  =  1  » N 

C( I )=C( I )+Q( I »J)*Y( J) 

208  CONTINUE 
207  CONTINUE 

C  A  MATRIX  OF  Q ( I » j ) *QT ( I  * J  )  IS  NOW  DETERMINED. 

DO  209  K= 1 » M 1 
DO  210  J  =  1  » M 1 
D ( K > J  )  =0  •  0 
DO  211  1  =  1  »N 

D(K*J)=D(K>J)+Q(K* I >#QT( I » J ) 

211  CONTINUE 
210  CONTINUE 


« :  =  u  » I )  o 

I— 1 3  c  ■  ■  i 


;  n)  =!i*  )T 


.  u  --  ( y  •  ' 


,*■*1  r xs  -c 


r  n  n  r\  o  n  r 


SUBROUTINE  LSQFIT 


209  CONTINUE 
21?  RETURN 
END 
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SUBROUTINE  GAUSS 

SUBROUTINE  GAUSS  (  X » Y » G » QT  * C > D > A  * DAUG * N  » M » Ml  ) 

REAL  X ( 50  *25 ) > Y ( 50 ) *Q ( 25  »  50 ) »QT ( 50  »  25 ) »C(25) *D ( 25  *25  ) » 
1 A ( 25 )  »  DAUG ( 2  5  »  2  5  ) 

INTEGER  N»M»M1 ,  I  »K»  j,  IM1  ,  IP1  ,M2  >KP1 

C  GAUSSIAN  ELIMINATION  METHOD  IS  USED  TO  EVALUATE  THE 
C  ELEMENTS  OF  THE  RESULTANT  COLUMN  MATRIX. THE  METHOD  HOLDS 
C  PRETTY  WELL  FOR  SMALL  SIZED  *  DENSE  COEFFICIENT  MATRICES. 

M2=M+2 


C  AN  AUGMENTED  MATRIX  WITH  C ( I )  FORMING  THE  M+2TH  COLUMN  OF 

C  D (  I  » j  )  IS  DERIVED. 

DO  310  1  =  1  * M 1 
DO  311  J  = 1  *  M 1 
DAUG { I  * J ) =D ( I » J ) 

311  CONTINUE 

312  DAUG ( I  *  M2 ) =C (  I  ) 

310  CONTINUE 

313  DO  300  1=1 >M1 
DM AX =0.0 

DO  301  K= I » Ml 

IF(ARS(D(K>I) ).LT.DMAX)GO  TO  301 
DMAX  =  ABS  (  DU*  I  )  ) 

I  STORE  =  K 
301  CONTINUE 

C  THE  ABOVE  STATEMENTS  SEEK  THE  MAXIMUM  ELEMENT  IN  A  COLUMN 
C  AND  THE  AIM  IS  TO  BRING  THE  MAXIMUM  ELEMENT  ON  THE 
C  DIAGONAL. THIS  MINIMISES  ROUND  OFF  ERROR. 


I F ( I.EQ. ISTORE)  GO  TO  302 
DO  30  3  L  = I  * M 2 
TEMP=DAUG ( I  *  L ) 

DAUG (  I »  L ) =DAUG ( ISTORE »L) 

DAUG ( I  STORE  »  L )=  TEMP 
303  CONTINUE 
302  IP1=I+1 

DO  304  K= I P 1  *  Ml 

RA T 1 0=  DAUG ( K  » I ) / DAUG ( 1,1) 

DO  305  L=  I  »M2 

I F ( L. EG.  I  )  GO  TO  308 

DAUG ( K  »  L ) =  DAUG ( K  *  L ) -RATI 0*DAUG ( 

GO  TO  305 

308  DAUG(K,L)=0.0 

305  CONTINUE 
304  CONTINUE 
300  CONTINUE 
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SUBROU TINE  GAUSS 

C  THE  AUGMENTED  MATRIX  NOW  CONSISTS  OF  AN  UPPER  DIAGONAL 
C  MATRIX  WITH  A  COLUMN  VECTOR  ATTACHED  AT  ITS  M+2TH  COLUMN 
C  THE  SOLUTION  FOR  THE  RESULTANT  VECTOR  (A)  STARTS  FROM  THE 
C  BOTTOM  TOP. 


A ( Ml ) =  DAUG (Ml »M2 ) /DAUG ( Ml *M1 ) 

DO  306  1  =  1  >M 
SUM=0 • 0 
K  =  M 1  -  I 
KP  1  =  K+1 

DO  307  J=KP1 *  M 1 

SUM  =  SUM+ A ( J ) *DAUG ( K  * J ) 

307  CONTINUE 

A  (  ’<  )  =  (  DAUG  (  K  *  M2  )  -SUM  )  /  DAUG  (  K  *  K  ) 
30b  CONTINUE 
RETURN 
END 


SUBROUTINE  STD 


SUBROUTINE  STD ( A , Y ,QT *  DEV »SD ♦ N »M1 ) 
REAL  A (25  )  *Y ( 50 )  >QT ( 50 ,25  )  »DEV ( 50  )  ,SD 
INTEGER  N*M1 
DO  4001=1 »N 
SUM=-Y ( I ) 

DO  4 0 1 J  = 1  *  M 1 

DEV ( I ) =  SUM+QT ( I »  J ) * A ( J ) 

401  CONTINUE 
400  CONTINUE 

SUMDSQ=0.0 
DO  40  2  1  =  1  *N 

SUMDSQ=SUMDSQ+DEV ( I ) **2 

402  CONTINUE 
SD1=SUMDSQ/FL0AT ( N-l ) 

SD  =  SQRT (SD1  ) 

RETURN 

END 
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SUBROUTINE  GAQ10 

SUBROUTINE  GAQ10 (XL»XU»FCT»Y,A»C»P) 

REAL  X  L  » X  U  » Y  » A » C *P ( 25 ) 

EXTERNAL  FCT 
A 1 "0 • 5* ( X  U+XL ) 

8=XU-XL 

C1=0. 4869533*8 

Y=0*Q3333567*(FCT(  (Al+Cl)  >A»C>P)+FCT(  (A1~C1)  *  A  *  C  »  P  )  ) 
0=0.4325317*8 

Y=Y+0.0 7472567 *( FCT (  (Al  +  Cl)  »A»C»P)+FCT(  (Al-Cl)  »  A  »  C  >  P  )  ) 
C1=0. 3397043*B 

Y=Y+0. 1095432* ( FCT ( (Al+Cl ) >A»C»P)+FCT( (Al-Cl ) »A»C»P) ) 
C1=0. 2166977*B 

Y=Y+0 . 1346334* ( FCT ( (Al+Cl) »A,C,P)+FCT( (Al-Cl) ,A»C»P) ) 
C1=0. 07443717*8 

Y  =  B* ( Y+0 •  1477621* ( FCT (  (Al+Cl)  »A>C»P)+FCT(  (Al-Cl)  >A*C 
1  *  P  )  )  ) 

RETURN 

END 


FUNCTION  FCT 1 


FUNCTION  FCT1 (  Z  >  A  ♦  C  >  P  ) 

REAL  Z  » A  »  C  »P ( 25 ) 

B1=-(A**2 ) *S  I  N (Z  ) 

R2=C  + ( COS ( Z ) **P ( 2  )  ) 

B3= ( SIN (Z ) **2 )*(  (C+( P ( 2 ) +1 )*( COS ( Z ) **P ( 2  )  )  )  **2  ) 
R4=(P(2)*(SIN(Z) ) * ( COS ( Z ) ** ( P ( 2 ) - 1 ) ) -C*COS ( Z ) 
1-COS ( Z  )  ** ( P ( 2 ) +1  ) 

1**2 

FCT1=B1*B2*SQRT ( B3+B4) 

RETURN 

END 
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FUNCTION  FCT2 


FUNCTION  FCT2 ( Z  »  A  ♦  C  » P ) 

REAL  Z  »  A  ♦  C  » P  (  2  5  ) 

B1=-(A**3)*(SIN(Z)**3) 

B2=(C+(COS(Z)**P(2)  )  ) *  *  2 

33  =  C  + (  (P(2)+l)*(COS(Z)**P(2)  )  ) 

FCT2=B1#B2*B3 

RETURN 

END 
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-  FORMING  - 


-  FORMING  - 

REAL  F V ( 50  )  » V ( 50  > » OR ( 4 )  »  PH  I (50)  *  C ( 10)  ♦ ERR ( 50 )  »  T ( 1 00 )  » 
IS ( 100)  » VOL ( 100 )  » DER 1(50)  , FXV(50)  ,FXV1(50)  ,Y(50)  ,G,SP» 
2TIMER #U»D»D1 , INTI » INT2 ,H  »APEX  *PHI 1 ( 50)  » ERROR 
INTEGER  N 1  *  N  » M  *  M 1 »  L 1  *  L  *  F  V 1 ( 50  )  *V1 ( 50 )  »0R1 (4)  ,<  *NAPEX * 
1L2 (50  )  » K 1 (20  ) 

C  THE  PROGRAM  EVALUATES  THE  SURFACE  AND  VOLUME  OF  DESCRETE 
C  BUBBLES  DURING  FORMATION  AT  A  GIVEN  ORIFICE  UNDER 
C  SPECIFIED  CONDITIONS  OF  GAS  AND  LIQUID  FLOW. THE  TECHNIQUE 
C  INVOLVES  NUMERICAL  ANALYSIS. USE  IS  MADE  OF  THE  LAGRANGIAN 
C  POLYNOMIAL  FOR  DER I  VAT  I VES , FORSYTHE  ORTHOGONAL  POLYNOMIALS 
C  AND  NEWTON  COTES  INTEGRATION  FORMULAE. 

READ ( 5,92) ERROR 
92  FORMAT ( E 10  •  2  ) 

RE AD ( 5,95 ) NUMRUN 
95  FORMAT (14) 

NRUN=  1 

94  READ( 5 ,99 )G»U,SP»D, TIMER,  (C(I)»I  =  1,4) 

99  FORMAT (2F9.5 ,F12.2,F9«5 , 17  ,  4F7.2) 

L  =  1 

READ ( 5 , 1 OC  )  L 1 
IOC  FORMAT (14) 

101  READ( 5 ♦ 102 ) Ml 

102  FORMAT (14) 

M=  1 

98  READ( 5 ,103  )N»N1 

103  FORMAT (21 3) 

K  =  N 1  /  4 

K  1  (  M  )  =  K 


C  SO  FAR  THE  TOTAL  NUMBER  OF  BUBBLES  LI, THE  TOTAL  NUMBER  OF 
C  FRAMES  PER  BUBBLES  Ml,  THE  NUMBER  OF  INTERVALS  N 1  *  T  H  E  TOTAL 
C  DATA  POPULATION  N  THE  NUMBER  OF  PANELS  IN  THE  NEWTON  COTES 
C  FORMULA  K  AND  THE  NEWTON  COTES  COEFFICIENTS  HAVE  BEEN 
C  READ. THE  DATA  FROM  EACH  BUBBLE  FRAME  WILL  MOW  BE  READ. 

RE AD ( 5,93) L2 (M) 

93  FORMAT (15) 

READ( 5,104)  ( OR1 (  I  )  , 1  =  1 ,4) 

104  FORMAT (41 5) 

FOR  REASONS  WHICH  WILL  BE  CLEAR  LATER, A  THIN  ALMOST 
SPHERICAL  SECTOR  FROM  THE  BUBBLE  APEX  WILL  BE  CUT  OFF  AND 
TREATED  SEPARATELY  FOR  SURFACE  AMD  VOLUME. 

READ ( 5,97) ( V 1 ( I )  »FV1(I)»I=1,N) 

97  FORMAT (1615) 
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C  THE  APEX  VALUE  WILL  BE  READ  SEPARATELY. NOTE  THAT  THE  VALUE 
C  OF  FV  AT  APEX  IS  BEING  ASSUMED  TO  BE  EQUIVALENT  TO  ZERO* 

C  THE  POINT  BEING  ON  THE  AXIS. 

READ (5*96) NAPEX 
9 fa  FORMAT! 15) 

C  THE  DATA  WILL  NOW  3E  TRANSFORMED  THROUGH  AN  AXIS  SHIFT  AND 
C  CORRECTED  FOR  MAGNIFICATION  AND  DISTORSION. 

C  FLOATING  DIGITAL  DATA. 

DO  105  1=1 >N 

FV (  I  )=  FLOAT ( F V 1 (  I )  ) 

V(  I  )  =FLOAT ( VI ( I )  ) 

105  CONTINUE 

DO  111  1=1*4 
OR (  I  ) =  FLOAT ( OR 1 (  I )  ) 

111  CONTINUE 

APEX  =  FLOAT(NAOEX) 

D 1 =OR ( 4 ) -OR ( 2 ) 

DO  106  1  =  1  iN 
FV ( I ) =  FV (  I  )*D/D1 
V(I)=V( I )*(D/D1)*( 3.48/3.42) 
lCfa  CONTINUE 

AP-EX  =  APEX*(D/D1  )  *  (  3 .48/3 . 42  ) 

EVALUATING  MEAN  Y  COORDINATE  FOR  TWO  DATA  POINTS  IN  THE 
SAME  HORIZONTAL  LEVEL . CONS  I DER I NG  TOTAL  NUMBER  OF  POINTS 
TO  BE  ALWAYS  EVEN. 

DO  107  I  =  2  »  N  *  2 
I M 1=1-1 

V  {  I  )  =  (  V  (  I  )  + V  (  I M 1  )  )  /  2  •  0 
V (  I  Ml ) = V (  I  ) 

107  CONTINUE 

TRANSFORMING  COORDINATES  BY  AXIS  SHIFT. 

SH I FT  =  V ( 1  ) 

SHIFT1= ( F V ( 2 ) +FV ( 1  )  )/2.0 
DO  108  1  =  1  »N 
FV ( I ) =FV(  I ) “SH I F  T 1 
V ( I ) = V ( I ) -SH I  FT 

108  CONTINUE 
APEX=APEX-SHIFT 

THE  ABOVE  TRANSFORMATION  GIVES  A  POSITIVE  AND  A  NEGATIVE 
VALUE  OF  FV  FOR  EACH  VALUE  OF  V » THE  TWO  VALUES  MUST  BE  OF 
THE  SAME  ABSOLUTE  MAGNITUDE. 


. 
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DO  109  I  =  2  *  N  »  2 
I M 1=1-1 

SUM= ( ABS ( FV (  I  )  ) +ABS ( FV ( I  Ml  )  )  ) /2  •  0 
FV ( I ) =SUM 
FV ( I M 1 ) = - SUM 
109  CONTINUE 

C  THE  DATA  WILL  NOW  BE  SMOOTHED  BY  THE  USE  OF  THE  FORSYTHE 
C  ORTHOGONAL  POLYNOMIAL. 


N  2  =  N  /  2 

CALL  FORSYT (FV»V*PHI *  ERR » Y , ERROR ♦ N2 * N » N 1 ) 

C  THE  SMOOTHED  DATA  HAS  PHI(V)  INSTEAD  OF  F ( V )  AND  FURTHER 
C  WORK  IS  DONE  WITH  THE  SMOOTHED  DATA. THE  DERIVATIVE  OF  THE 
C  SMOOTHED  FUNCTION  AT  DESCRETE  POINTS  WILL  NOW  BE  DERIVED 
C  USING  LAGRANGIAN  I NTERPOL AT  I  ON . THE  SPHERICAL  SEGMENT 
C  CONTAINING  THE  BUBBLE  APEX  POINT  WAS  TREATED  SEPARATELY 
C  BECAUSE  THE  DERIVATIVE  AT  THE  APEX  IS  INFINITE  AND  HENCE 
C  NOT  ACCEPTABLE. AN  APPROXIMATION  FOR  THE  SEGMENT  WILL  BE 
C  USED. 


CALL  LADERI (PHI »ERR*Y*DERI »N2*N) 

C  HAVING  EVALUATED  THE  DERIVATIVE  WE  SHALL  EVALUATE  THE 
C  FUNCTION  FXV. 

C  WHERE  FXV=F(V)*SQRT(1+DERI**2) 

DO  110  I  =  2  » N  »  2 

FXV( I )=ABS (PHI ( I ) )*SQRT( 1.0+DERl ( I ) *#2 ) 

110  CONTINUE 

C  THE  FUNCTION  FXV  IS  A  COMPOSITE  FUNCTION  FOR  EVALUATION  OF 
C  BUBBLE  SURFACE.  INTEGRATION  IS  BEING  DONE  BY  THE  USE  OF  A  K 
C  PANEL  NEWTON  COTES  FORMULA  OF  FOURTH  ORDER  I.E.A  FIVE 
C  POINTS  FORMULA  APPLIED  TO  4K+1  EQUIDISTANT  POINTS. 

D I  ST  =  AP EX -V ( 1 ) 

C  AN  ANGLE  OF  FIVE  DEGREES  ON  EITHER  SIDE  OF  THE  MEDIAN  OR 
C  CENTRAL  AXIS  OF  THE  BUBBLE  IS  TAKEN  AND  A  CHORD  JOINING 
C  THE  TWO  POINTS  OBTAINED  CORRESPONDS  TO  THE  PORTION  OF  THE 
C  BUBBLE  SLICED  OFF  FROM  THE  APEX  FOR  SEPARATE  TREATMENT. 

c  this  vertical  height  corresponds  to  about  0.2  percent  of 
c  the  total  bubble  height. 

SLICE=DIST*0. 2/100.0 

REMAIN=DIST -SLICE 

HE  I GHT  =  REMA IN /FLOAT (Nl ) 
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C  'HEIGHT* IS  THE  INTERVAL  LENGTH  OBTAINED  BY  DIVIDING  THE 
C  REMAINING  VERTICAL  DISTANCE  INTO  N1  PARTS. THE  EGUISPACED 
C  DATA  CAN  THUS  BE  GENERATED. 

CALL  LAGI NT (FXV»Y* APEX .HEIGHT  #N2  »N»N1 ) 

CALL  NCOTES (HEIGHT  » F  X  V  »  C  » INTI >  K  >  N  >  N 1 ) 

S  (  M  )  =  2  . 0*3. 141592  6* INTI 

C  THE  ABOVE  IS  THE  SURFACE  OF  THE  BUBBLE  EXCLUDING  THE  SMALL 
C  SEGMENT  AT  THE  APEX. IN  THIS  REGION  THE  DERIVATIVE  IS  VERY 
C  LARGE  AS  COMPARED  TO  1  AND  THE  EQUATION  FOR  THE  SURFACE 
C  MAY  BE  APPROXIMATED  TO  THE  AREA  OF  A  DISC  WITH  RADIUS 
C  F  V  (  N  )  * 

CALL  LAG  I  NT (PH  I » Y.APEX .HEIGHT »N2 .N.N1 ) 

51  =  3.1415926* (PH  I ( N ) **2 ) 

S  (  M  )  =  S  (  M  )  +  S 1 


C  IF  THE  FRAME  CORRESPONDING  TO  THE  FINAL  DETACHED  BUBBLE  IS 
C  UNDER  CONSIDERATION  WE  HAVE  TO  ADD  THE  AREA  OF  THE  BOTTOM 

C  DISC  ALSO. 


IF ( M. LT.M1  )  GO  TO  112 
S2=3.1415926*(PHI (2)**2) 

S ( M ) =S ( M ) +S2 

C  HAVING  EVALUATED  THE  SURFACE  WE  SHALL  NOW  EVALUATE  THE 
C  VOLUME  OF  THE  BUBBLE  AT  THAT  TIME. 

112  DO  114  I  =  2  >  N  ♦  2 

FXV1 (  I  ) =  P  H I  (  I  ) **2 
114  CONTINUE 

CALL  NCOTES (HEIGHT  » FX V 1 » C » I NT2 » K » N » N1  ) 

VOL (M ) =3. 1415926*INT2 

THIS  ONCE  AGAIN  EXCLUDES  THE  VOLUME  OF  THE  SEGMENT. THAT 
VOLUME  IS  BEING  APPROXIMATED  TO  THAT  OF  A  SPHERICAL 
SEGMENT. 

HI =SL I CE 

V2=(3.14l5926*Hl/6.0)*( (3. 0*PHI ( N ) **2 ) + ( HI** 2 ) ) 

V  0  L  (  M  )  =  V  0  L  (  M  )  +  V  2 


C  THE  VOLUME  IS  THUS  EVALUATED. WE  SHALL  EVALUATE  THE  TIME  AT 
C  WHICH  THE  ABOVE  SURFACE  AND  VOLUME  EXISTS. 

C  NOTE  THAT  THE  EXACT  SPEED  MUST  BE  READ  IN. 

T ( M ) =  FLOA T ( L 2 ( M )  ) /SP 
IF(M.GE.Ml)  GO  TO  115 

M  =  |V  +  1 


Y 


. 
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GO  TO  98 

115  IF(L.GT.l)  GO  TO  128 
WRITE (6  » 129 ) NRUN 

129  FORMAT (' 1 ' »20X» ' RESULTS  FOR  RUN  NO.  M='»I4) 

128  WRITE (6  *  116 ) L 

116  FORMAT ( ‘0 » »10X* ' TIME  VS  SURFACE  AND  VOLUME  FOR  BUBBLE 

1  NO.  =  '  »  14  ) 

WRITE(6»117)G»U»SP*D 

117  FORMAT (  '0 •  »10X» 'G  CC/SEC= '  » F10 . 5 » 5X »  ' U  C  M  /  S  E  C  =  '  »  F 1 0 • 5 / 
HOXf'SP  FPS=  '  »F10.2  »5X  »  '  D  CM='»F10.5) 

WRITE (6  » 123 ) 

123  FORMAT ( '0 ' >10X> ' FRAME  N0.'»4X»'T  SEC'»4X»'S  SQ.CM' » 

14X , ' V  CC '  *4X  .  » PANELS  K '  ) 

DO  124  1=1 »M1 

WRITE(6 >125)L2(  I  )  *T( I  )  *S(  I  )  >VOL( I  )  » K1 (  I ) 

124  CONTINUE 

125  FORMAT ('  '  *  1 IX  *  I  4  * 5X » F10 . 6  * 2X » FI 1 . 6  *  IX » F9 . 6 »4X *  I  2 ) 

IF ( L • GE • L 1  )  GO  TO  126 

L  =  L+ 1 
GO  TO  101 

126  I F ( NRUN • GE • NUMRUN )  GO  TO  131 
NRUN=  NRUN  +  1 

GO  TO  94 
131  STOP 
END 
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SUBROUTINE  FORSYT 


SUBROUTINE  FORSY  T ( FV  » V  »  PH  I » ERR » Y » ERROR » N2 >N*N1 ) 

REAL  FV ( 50 ) »V ( 50 ) »PHI ( SO ) >ERR ( 50 ) »Y ( 50 ) » ALPHA ( 20 ) » 
1BETA(20) >B(20) >PHI1(20»50)  »X(50)  » P ( 20  »  50 )  »ER(2G»5Q) » 

2 VAR » ERROR 

INTEGER  N  2  »  M  »  N1 » I  STORE 

C  WE  SHALL  FIRST  EVALUATE  THE  ORTHOGONAL  POLYNOMIALS  AND 
C  THEN  THE  COEFFICIENTS  OF  THE  NORMAL  EQUAT I ONS . FOR  DETAILS 
C  SEE  LAP  IDUS (DIGITAL  COMPUTATION  FOR  CHEM . ENGRS . ) PP333-336 
C  REDEFINING  INDEX  OF  DATA 

DO  9  I  =  1  *  N2 
X  (  I  )  = ASS ( F V ( 2 *  I  )  ) 

Y(  I  )-V(2*I ) 

9  CONTINUE 

DO  10  1  =  1  *N2 
P  (  1  »  I  )  =0.0 
P( 2*  I  )  =  1  -  0 
PHI  1  ( 1 , I ) =0.0 

10  CONTINUE 

C  THE  FIRST  TWO  POLYNOMIALS  HAVE  BEEN  DEFINED. THE  MAXIMUM 
C  DEGREE  OF  POLYNOMIAL  TO  BE  USED  IS  12. 

DO  11  J=  2  » 14 
SI =0.0 

52  =0 • 0 

53  =  0  •  0 

54  =  0  •  0 
S5=0 • 0 

DO  12  1=1 »N2 
S1=S1+Y ( I  ) *P ( J  *  I  ) **2 

52  =  S2  +  P (  j  *  I  ) 

53  =  S3+Y ( I ) *P ( J- 1  *  I ) *P ( J  *  I  ) 

54  =  S4  +  P ( J-l » I  ) **2 

55  =  S5  +  X(  I  ) *P ( j  *  I  ) 

12  CONTINUE 

C  DEFINING  ALPHAS  AND  BETAS 

ALPHA ( J+l )=S1/S2 
IF ( J.EQ.2 )  GO  TO  5 
BETA ( J ) =S3/S4 
GO  TO  6 

5  BETA(2)=G*0 

6  B(J)=S5/S2 

C  DEFINING  P( j+1 ) 


DO  13  1  =  1 »  N 2 


*  ♦ 
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,  ,  U  .  >  Y  +  J  *12 


VI  -v-)  C  I  )  X+  !v^ 


.  =(  ) 


«I*I 


B-  30 


SUBROUTINE  FORSYT 

P ( J+l » I ) = ( Y ( I ) -ALPHA ( j+1 ) ) *P ( J» I ) -BETA ( J ) *P ( J-l , I ) 

13  CONTINUE 

DO  14  1=1 >N2 

PHI1( J»I )=PHI1(J-1*I)+B( J)*P( J*I ) 

ER ( J  » I ) = ABS ( X ( I ) -PH  1 1 ( j» I )  ) 

ER R ( I  ) =E R ( J  » I )**2 

14  CONTINUE 
S  =  0.0 

DO  15  1  =  1 » N2 
S  =  S  +  ERR ( I  ) 

15  CONTINUE 
M=u-2 

I  STORE  =  U 
IN 2  2  =  N  2-M 

I F ( N 2 2 • E Q • 1  )  GO  TO  105 
VmR=S/FLOoT(N2-(M+1 ) ) 

105  I  F ( N2 2 • EQ • 1 )  GO  TO  16 

IF(VAR.LE. ERROR)  GO  TO  16 
11  CONTINUE 

C  WE  SHALL  NOW  REDEFINE  THE  VALUES  THAT  HAVE  TO  RETURN  OUT 
C  OF  THE  SUBROUTINE. 

16  DO  17  1  =  1  *  N 2 

PHI  (  I  )  =P H  1 1  (  I  S T 0 R. E  *  I  ) 

17  CONTINUE 

C  SINCE  THIS  DATA  GOES  DIRECTLY  INTO  THE  NEXT  SUBROUTINE  WE 
C  DO  NOT  NEED  TO  REINDEX  AT  THIS  STAGE, BUT  SHALL  DO  SO  AT 
C  THE  END  OF  THE  SUBROUTINE  LADERI. 

18  RETURN 
END 
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SUBROUT  I  NF  LADERI 

SUBROUTINE  LADERI  (  D  H  I *ERR*Y»DERI » N  2  *  N ) 

REAL  PHI ( 50 )  »ERR ( 50 ) »Y( 50 )  ♦  DER I ( 50 )  >D ( 50  >  50  )  »POLY ( 50 ) 
INTEGER  N 2  » N 

C  THIS  SUBROUTINE  IS  USED  TO  EVALUATE  THE  DERIVATIVE  AT 
C  INDIVIDUAL  POINTS  OF  A  TABULATED  DATA  SET. THE  METHOD  USED 
C  IS  THE  LAGRANGIAN  INTERPOLATION. 

C  EVALUATING  LAGRANGIAN  POLYNOMIALS  AND  THE  D(IK)'S 

DO  201  1  =  1  *N2 

POLY ( I ) = 1 • 0 

DO  202 J= 1 »N2 

I  F (  I.EQ. J)  GO  TO  202 

POLY ( I )=POLY( I )*{ Y( I )-Y( J) ) 

202  CONTINUE 
201  CONTINUE 

DO  203  K= 1  ♦  N2 
DO  20 A  1  =  1 » N 2 
I  F ( I.EQ.K)  GO  TO  204 
D(I»K)  =  (Y(K)“Y(I)) *POLY  (  I  ) 

204  CONTINUE 

203  CONTINUE 

C  EVALUATING  DERIVATIVES  ONE  BY  ONE. 

DO  20  5  K  =  1  *  N  2 
PROD=l .0 
SUM1=0 .0 
DO  206  J= 1 » N2 
IF(K.EQ.J)  GO  TO  206 
PROD=PROD* ( Y ( K ) - Y ( J ) ) 

SUM1=SUM1+ ( 1.0/(Y(K)-Y(J) ) ) 

206  CONTINUE 
SUM 2=0 .0 

DO  207  1  =  1  *  N 2 

IF ( K.EQ. I )  GO  TO  207 

SUM2=SUM2+(PHI ( I ) /D( I *<) )*PROD 

207  CONTINUE 

DERI ( K ) =SUM2+PH I (K)*SUM1 

205  CONTINUE 

C  THE  VALUES  OF  PHI (I)  AND  DERI (I)  AND  ERR ( I )  WILL  NOW  BE 
C  REINDEXED  TO  SUIT  THE  MAIN  LINE  PROGRAM. 

DERI ( N ) =DER I ( N/2 ) 

PHI (N ) =PH I (N/2 ) 

ERR ( N ) =ERR ( N/2 ) 

NM2=N-2 

DO  209  K=2»NM2»2 
I  =N~K 


. 
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SUBROUTINE  LADERI 


PHI ( I  )=PHI  (1/2) 
DERI ( I ) =  D  E  R I ( 1/2 ) 
ERR ( I  )  =ER R  (1/2) 

209  CONTINUE 

210  RETURN 
END 
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SUB ROUT  IMF  LAG  I  NT 

SUBROUTINE  LAG  I N T ( FN » Y , APEX , HE  I GHT * \i2*N>Nl) 

REAL  FN(50) *Y(50)  » APEX > HE  I GHT > FCT ( 50 )  » FX V2 ( 50  )  > Y 1 ( 50  )  ♦ 
1L ( 30 ) 

INTEGER  N 2  ♦  N  » N 1 

C  THIS  SUBROUTINE  GENERATES  EQUISPACED  DATA  FROM  THE 
C  UNEQUI SPACED  DATA  AFTER  THE  COMPOSITE  FUNCTION  FXV  HAS 
C  BEEN  EVALUATED. A  THIRD  ORDER  LAGRANGIAN  INTERPOLATION 
C  POLYNOMIAL  WILL  BE  MADE  TO  FIT  4  POINTS  AND  SOME 
C  EOU I  SPACED  DATA  DETERMINED . THE  FIRST  IS  THEN  DROPPED  AND 
C  THE  5 TH  POINT  IS  ACCEPTED  AND  A  NEW  POLYNOMIAL  TAKEN  TO 
C  GIVE  DATA  FOR  NEXT  SET  OF  INTERVALS  AND  SO  ON. 

C  REINDEX  FN 

DO  10  1=1 >N2 
FXV2 ( I ) =FN ( 2*1 ! 

10  CONTINUE 

N  2  P 1  =  N  2  +  1 
Y ( N2P 1 ) = APEX 
FXV2 ( N2P1 ) =0.0 

C  THIS  POINT  OF  FXV 2  IS  DEFINED  SO  THAT  IN  EVALUATION  OF 
C  THOSE  VALUES  OF  FCT  WHERE  THE  FINAL  APEX  POINT  IS  REQUIRED 
C  WE  MAY  NOT  DRAW  ON  FICTITIOUS  VALUES  FROM  COMPUTER  MEMORY. 
C  GENERATING  NEW  VALUES  OF  Y  WHICH  ARE  EQUISPACED. 

Y 1  (  1 )  =  Y  (  1  ) 

DO  11  1=2 »N2 

I  v!  1  =  I  - 1 

Y 1  (  I  )  = Y 1 (  I  Ml ) +HE I GHT 

11  CONTINUE 
KK=  1 

II  =  1 

13  KK 1=KK+1 
KK  2  =  KK+2 
KK3=KK+3 
DO  12  I  =  I  I  » N 2 
IP1=I+1 

IF ( Y ( KK3 ) . LT. Y1 ( I ) )  GO  TO  16 

THE  ABOVE  IS  A  SAFETY  DEVICE  TO  CHECK  IF  THE  NEW  POINT 
FALLS  WITHIN  RANGE  OF  THE  POLYNOMIAL  BEING  USED. 

FCT ( I ) =0 • 0 
DO  14  K=KK  » KK3 
PROD= 1.0 
PROD1 = 1.0 

C  EVALUATING  LAGRANGIAN  COEFFICIENTS. 
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SUBROUTINE  LAGINT 

DO  15  J=KK»KK3 
IF(K.EQ.J)  GO  TO  15 

C  K  MAY  NOT  BE  EQUAL  TO  j  BECAUSE  THE  LAGRANG I  AN  COEFF.IS 
C  NOT  DEFINED  AT  J.EQ.K 

PROD=PROD*( Y1 ( I ) -Y ( J ) ) 

PRODl=PRODl*( Y ( K ) -Y ( J ) ) 

15  CONTINUE 

L ( K ) =PROD/PRODl 
FCT ( I  )  =FCT ( I ) +L ( K ) *FXV2 ( K ) 

14  CONTINUE 

IF (  I  .EQ.N2  )  GO  TO  12 
IF ( Y ( KK3 ) .GT. Y1 ( I P 1 ) )  GO  TO  12 

THE  ABOVE  STATEMENTS  IMPLY  THAT  IF  THE  NEXT  NEW  POINT  IS 
WITHIN  THE  RANGE  OF  THE  EXISTING  POLYNOMIAL  THEN  USE  THE 
SAME  POLYNOMIAL  INSTEAD  OF  A  NEW  ONE. 

11=1+1 

C  INCREMENTING  KK  TO  ACCEPT  THE  NEXT  POINT  AND  REJECTING  THE 
C  FIRST  ONE. 


16 

K  K  =  K  K  + 1 

GO  TO  13 

12 

CONTINUE 

DO  17  I  =  2  *  N  »  2 

FN (  I ) =FCT  <  I  / 2 ) 

17 

CONTINUE 

18 

RETURN 

END 

' 
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SUBROUTINE  NCOTES 

SUBROUTINE  NCOTES ( HF I GHT  » FXV »C » I  NT . K * N * N1  ) 

REAL  HEIGHT* FXV ( 50)  *C(  10)  »  I  NT  * H *  A 1 » B 1  * C 1 »D2»D3 

INTEGER  K»N#N1»K1*K2*K3 

H=HF I GHT 

A1=2.0*H/45.0 

B 1 = ( FXV ( 2 ) +FXV ( N ) ) 

0=0.0 
K 1  =  8  *  K 

DO  111  I  =  4  » K 1  » 4 
C 1 =C 1  +  FX  V ( I ) 

111  CONTINUE 
D2  =  0 • 0 
'<2=  (  8*K  )  -2 

DO  112  I  =  6  >  ;< 2  *  8 
D2  =D2  +  FX V ( I ) 

112  CONTINUE 
D3  =0  •  0 
K3= ( 8*K  )  -6 

DO  113  I  =  1 0  »  K3  »  8 
D3  =  D3  +  FXV (  I  ) 

113  CONTINUE 

INT=A1*(C ( 1 )*B1+C(2)*C1+C(3 )*D2+C(4)*D3 ) 

114  RETURN 

END 
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-  SAMPLE  “ 


-  SAMPLE  " 


C  THE  VALUE  OF  A  REPRESENTATIVE  SAMPLE  SIZE  FROM  A  DATA 
C  POPULATION  OF  SIZE  N  IS  TO  BE  EVALUATED  USING  STATISTICS. 

REAL  X ( 250  )  » N 1  * S TDEV » MEAN  » UP * X 1 ( 2  50 ) 

INTEGER  N  >  L  »  L 1 »  M  *  M 1  *  I  *  K 1 


C  READ  THE  CODE  WHICH  IDENTIFIES  THE  DATA  SET. 


READ ( 5  » 7 ) Ml >  K 1 

7  FORMAT (213) 

21  L=  1 

19  READ ( 5  *  9 ) N  >  UP 

9  FORMAT ( I4#F10. 4 ) 

C  N  IS  THE  POPULATION  SIZE  AND  UP  THE  SIGMA  LIMIT. 

C  THE  INDIVIDUAL  DATA  POINTS  WILL  BE  READ. 

READ ( 5 ,10 )  (X  (I)  >  1  =  1 >N) 

10  FORMAT (4E20. 8 ) 

IF ( L.GT. 1 )  GO  TO  8 
READ ( 5 ,6 ) LI 
6  FORMAT ( IX  ♦  13 ) 

8  SUM=0.0 

DO  11  1  =  1  *N 

11  SUM=SUM+X ( I ) 

MEAN= SUM /FLOAT ( N ) 

S  U  M 1  =  0 . 0 

DO  12  1  =  1  #iM 
XI (  I  )  =X (  I  ) -MEAN 

12  SUM 1 = SUM 1+ ( X 1 ( I ) )**2 
STDEV=SQRT (SUM 1/FLO AT ( N“1 ) ) 

N 1 = ( (UP*STDEV) /( 0.01*MFAN) )**2 

C  N 1  IS  THE  REQUIRED  SAMPLE  SIZE. 

WRITE (6  *13) Ml *K1 »L 

13  FORMAT (' 1 '» 10X »' ORIFICE  N0.='>12»  5X»'LIQ.FLOW  NO.='» 
112  *5X  * 'RUN  NO.=  '  »  12  ) 

WRITE (6  » 14) 

1*  FORMAT (  'O'  »15X» '  INDIVIDUAL  DATA  POINTS  USED  *  X ( I  )  •  ) 
WRITE (6 >15)  (X(  I  )  »I  =  1*N) 

15  FORMAT ('  ' >10X>F10.6»3X>F10.6»3X»F10.6»3X»F10.6) 

WRITE (6 .16) MEAN 

16  FORMAT (' 0 '» 15X >' THE  MEAN  OF  TOTAL  POPULATION  I S= '  * F12 • 
18  ) 

WRITE (6 > 17 ) STDEV 

17  FORMAT ( '0 ' >15X» ' THE  STANDARD  DEVIATION  IS='»F12.8) 
WRITE (6  > 18 )N1 


(  ,  1,(1)  }  It  » '  )  ' ' 

I  I 


.1  '  ’  IT 

«  y-'.r  .  ix 
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-  SAMPLE  - 

18  FORMAT (  '0 '  »15X» ' SAMPLE  SIZE  REQU I RED= ■  >  F 10  •  3 ) 
IF(L.GE.Ll)  GO  TO  22 
L=  L  + 1 
GO  TO  19 
22  STOP 
END 
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Appendix  C 

The  appendix  is  subdivided  into  the  following 

sections : 

(1)  Rotameter  Calibration. 

(2)  Titration  Analysis  to  standardize  BaCOH)^  and  HC Z 
solutions  and  to  evaluate  C^  the  concentration  of 
carbon  dioxide  in  the  bulk  liquid  phase. 

(3)  Details  of  the  velocity  profile  measurements,  and 
calibration  of  the  probes  used  to  measure  point 
liquid  velocities  in  the  column. 

Calibration  of  Liquid  Flow  Rotameters 

Two  rotameters  R13M-25-3  and  R-10M-25-3  were  cal¬ 
ibrated  to  measure  liquid  velocities  ranging  from  0  to  4.06 
cm/sec  in  a  column  of  20.2  cms  internal  diameter. 

The  calibration  data  is  given  in  Table  7  and  the 
calibration  curves  are  shown  in  Figure  37. 
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Column  I.D.  =  20.2  cms :  Sectional  Area  =  320  cm^ . 
R13M-25-3 


Rotameter  Reading  Flow  Rate  Average  Liquid  Velocity 
_ _ (1/sec)  in  Column  (cm/sec) 


10 

0.135 

0.422 

20 

0.352 

1.10 

30 

0.545 

1.70 

40 

0.731 

2.28 

50 

0.910 

2.84 

60 

1.10 

3.44 

70 

1.30 

4.06 

■25-3 

10 

0 . 0496 

0.155 

20 

0.0878 

0 . 274 

30 

0.1230 

0.384 

40 

0.1610 

0.503 

50 

0.2010 

0.628 

60 

0.2430 

0.759 

70 

0.2860 

0.894 

80 

0.3250 

1.02 

90 

0.3730 

1.17 

100 

0 . 4190 

1.31 

Table-7  ROTAMETER  CALIBRATION  DATA 
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Figure  37 :  Calibration  Curves  for  Rotameters 
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Titration  Analysis 

(a)  Standardization  of  Ba(OH)2  and  HC1  solutions. 

The  standardization  of  Ba(OH)2  and  HC1  solutions  is 
done  through  electrometric  titrations  using  glass  and 
calomel  electrodes.  Bromothymol  blue  is  used  as  an  indicator 
and  the  experiment  is  conducted  in  an  atmosphere  of  nitrogen. 
The  end  point  corresponds  to  the  inflextion  in  a  curve  of 
pH  vs.  acid  concentration. 

The  Ba(OH)2  solution  is  calibrated  against  a  stand¬ 
ard  HC1  solution  of  strength  0.2N.  The  former  is  then 
used  to  standardize  the  HC1  solution  used  for  the  analysis 
of  CO^  in  water  during  the  mass  transfer  runs.  BaCOH)^ 
solution  has  to  be  standardized  frequently  because  its 
strength  changes  due  to  absorption  of  CO^  from  air.  The 
strength  of  the  HC1  solution,  however,  remains  essentially 
constant.  The  data  for  these  calibrations  is  given  in 
tables  8  and  9  and  Figure  38  is  a  typical  pH  vs.  acid  con¬ 


centration  curve. 


s  tc  l  lu.  o  ;  :oh  bna  r  H  ))bH  lo  no  a  b^u  tc. 


t ;  a  d  •  >  3  n  .  n  CO)  >£  ’ 


/ 
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Figure  38 :  A  Typical  pH  Titration  Curve 


h. 


CO 

V  =  Volume  of  HC£  added  (cc) ,  Sample  Volume  =  20  m l  (BaCOH)^), 

Indicator  =  Bromothymol  Blue,  G-Y  =  Greenish-yellow,  Y-G  =  Yellowish-green. 
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Table  8  -  Standardization  of  Ba(OH)0  against  standard  0.02N  HC£ 
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Table  9  -  Standardization  of  HC£  Solution 
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Calculations 

!  — 

Ba( OH) 2 

:  Volume 

=  20  ml. 

Run  1  HC1 

0.02N 

PH  6.95 

HC1  Volume 
at  End  Point 

6.60  cc 

Run  2  ” 

t? 

"  7.00 

n 

6.65  " 

Run  3  " 

Tt 

"  7.10 

I! 

6.65  " 

Average 

6.63  cc 

Ba(OH)2 

:  Volume 

=  20  ml. 

Run  1  HC1 

Unknown 

PH  6.70 

HCL  Volume 
at  End  Point 

19.80  cc 

Run  2  " 

n 

"  7.00 

n 

0 

00 

• 

1 — 1 

Run  3  " 

n 

"  7.00 

1! 

19.75  " 

Average 

19.78  cc 

Normality  of 

Ba(OH)2 

=  6.63  x 

0 . 02N 

20 

Normality  of  HC1  =  6.63  x  0.02N 

20  19.78  =  x  0.02N 

=  0.0067  N. 

"b 

(b)  Titration  Analysis  to  Evaluate 

The  prestandardized  Ba(OH)2  and  HC1  solutions  are 
used  to  determine  the  concentration  of  CO^  in  the  bulk 
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water  phase.  A  known  volume  of  water  containing  dissolved 
CC>2  is  added  to  a  known  volume  of  the  barium  hydroxide 
solution.  The  Ba(OH)2  reacts  with  CO^  in  water  to  give 
BaCO^.  Addition  of  HC1  during  titration  serves  the  fol¬ 
lowing  purposes : 

(i)  To  neutralize  excess  Ba(OH)2  to  BaCl^ 

(ii)  To  convert  the  BaCO^  to  BaCl^  with  evolution  of 


co2. 

The  conversion  of  BaCO  to  BaCl 
namely : 


takes 


place 


in  two  stages 


(i)  conversion  of  BaCO^  to  Ba(HC03)2  and  BaCl^.  This 
corresponds  to  a  change  in  color  of  the  bromothymol  blue 
indicator  from  blue  to  green. 

(ii)  conversion  of  BaCHCO^^  to  BaCl^  and  CO^ .  The 
indicator  changes  color  from  green  to  yellow. 

The  chemical  reactions  involved  in  the  complete  titration 
are  given  as  follows : 


Ba(OH)2  +  2HC 1  ->  BaCl2  +  2H2  0 


(C-l) 


The  above  given  reaction  corresponds  to  the  blank  run  where¬ 
in  no  C02  is  involved.  The  bromothymol  end  point  is  from 
blue  to  yellow. 


+  BaC03  +  H2  0  +  Ba (OH) 2 

(excess ) 


Ba (OH) 2  +  C02 


(C-2) 


V 

‘  riT  .  <  DO  bnfi  I0&3  ol  1  DOH)sS  nols'isvhoo  12) 

■  '  I  ■~iy;  07  H .  £  n:.-  ' 
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This  corresponds  to  CO^  neutralization. 


Ba(OH)  ^  +  2HC1  -+  BaCl2  +  2H20 - (C-3) 

( excess ) 

2BaC03  +  2HC 1  +  Ba(HC03)2  +  BaCl2 - (C-4) 


These  reactions  correspond  to  the  end  of  the  first  stage  of 
the  back  titration.  The  bromothymol  blue  end  point  is  from 
blue  to  green. 

Ba(HCC>3 )  2  +  2HC 1  +  BaCl2  +  2H20  +  2C02 - (C-5) 

This  corresponds  to  the  end  of  the  second  and  final  stage. 
The  bromothymol  blue  end  point  is  from  green  to  yellow. 

If  the  volume  of  HC1  used  to  neutralize  a  given 
volume  of  Ba(OH)2  in  a  blank  titration  is  V^m£  and  the 
volume  used  at  the  end  of  the  first  stage  of  the  back  tit¬ 
ration  is  V2m-£,  then  (V^-V^) m-t  is  a  measure  of  the  C02 
concentration  in  the  water  sample.  This  can  be  shown  as 
follows : 

Say  that  the  initial  volume  of  Ba(OH)2  used  is  =  20  m£. 

Say  volume  of  HC1  required  to  neutralize  20  m t  of  Ba(0H)2 
=  V^m  l 
then 

20  x  N  =  V1  x  H 


( C — 6  ) 
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Where  N  and  are  the  normalities  of  BaCOH)^  and  HC£  res¬ 
pectively  . 

Say  volume  of  HCf.  used  at  the  end  of  the  first  stage 
is  V^mf..  includes  the  volume  for  neutralizing  excess 

Ba(OH)2  and  converting  the  barium  carbonate  to  barium 
bicarbonate.  The  total  volume  of  HCf  used  at  the  end  of 
both  stages  must  equal  because  at  the  end  of 

second  stage  effectively  all  the  BaCOH)^  has  been  consumed. 

In  a  blank  titration  the  conversion  of  BaCOH)^  to  BaCl^ 
takes  place  directly  and  in  the  titration  involving  CO^, 
the  conversion  to  BaCl^  takes  place  through  BaCO^ .  The 
volume  of  HCf  used  in  both  cases  for  the  same  initial 
volume  of  Ba(OH)2  must  be  the  same. 

Thus  the  volume  of  HCf  required  to  convert  barium 
bicarbonate  to  barium  cholride  = 

.*.  Volume  of  HCf  required  to  neutralize  total  BaC0^  =  2  (V-^-V^  )m£ . 
This  is  so  because  equal  volumes  of  HC l  are  used  at  each  step 
of  the  barium  carbonate  conversion  to  BaCl^. 

If  the  absorbed  CO^  is  taken  from  a  20  mZ  Sample 
of  water  then 

20  x  Nco  =  2(Vi"V2)  x  Ni - (C-7) 

Where  NCQ  is  the  normality  of  CO^  in  water. 

If  C1?  is  the  concentration  of  C0o  in  water  in  ?? S 
A  2  litre 
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then 


or 


nco2  2  X  CA 


. * .  20  x  2  x  C^  =  2(V  -V2)  x  N 


,b  _  (V1  V2)  X  N1 


A 


20 


(C-8 ) 


or  =  f[(V1-V2)] 

A  sample  pH  curve  for  a  blank  and  CC>2  titrations  is  shown 
in  Figure  39 . 


Measurement  of  Liquid  Velocity  Profile. 

The  measurement  of  liquid  velocity  profile  in  the 
column  involves  the  following  experiments. 

(i)  Determination  of  the  limiting  current  voltage  range, 

(ii)  Calibration  of  the  platinum  probes. 

(iii)  Measurement  of  point  velocities  in  the  column 
cross  section. 

Determination  of  the  Limiting  Current  Voltage  Range. 

The  effect  of  liquid  velocity  on  the  limiting  current 
and  the  voltage  range  for  which  the  latter  persists  is  given 
by  a  curve  of  polarographic  current  vs.  applied  voltage  at 
different  liquid  velocities. 

A  typical  set  of  data  for  such  curves  is  given  in 
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Figure  39:  Blank  and  C0o  titration  curves 
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Table  10  and  the  curves  for  probes  1  and  2  are  shown  in 
Figures  40  and  4l. 

The  polarographic  current  is  of  the  order  of  yA  and 
measurement  of  a  current  of  this  magnitude  is  difficult. 

To  avoid  inaccuracy,  the  current  was  passed  through  a  high 
resistance  (10.8kfl)  in  series  and  the  voltage  (order  of  mV) 
across  it  was  measured  with  a  digital  voltmeter. 
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Constant  Resistance  in  Calibration  Circuit  =  10.8kft 


C-15 


i — l 

o 

i — i 

on 

VO 

o 

1 - 1 

1 - 1 

CO 

VO 

< 

• 

• 

1 

1 

i — i 

1X3 

ov 

CM 

1X3 

VO 

VO 

C" 

CM 

1 

1 - 1 

i — I 

CM 

CM 

CM 

CM 

CM 

on 

n  un 

E- 

. 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

& 

• 

• 

• 

> 

1 

1 

CM 

-=3" 

1X3 

o 

LP\ 

O 

on 

CO 

o 

CM 

1 

S 

i — 1 

CO 

VO 

1 — 1 

-=3" 

C'- 

oo 

oo 

o 

1X3 

1 — 1 

CM 

CM 

CM 

CM 

CM 

on 

on 

o 

CT\ 

03 

-=r 

o 

on 

1X3 

30 

on 

ov 

c— 

-=r 

03 

CO 

o 

VO 

1 — 1 

< 

• 

• 

+3 

;a 

o 

C\J 

on 

03 

on 

t— 

o 

1 - 1 

CM 

on 

on 

LX3 

o 

cd 

i — 1 

1 — 1 

CM 

CM 

CM 

CM 

CM 

CM 

on 

II  o 

0) 

<=0  • 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

O 

o 

cd 

D  on 

-p 

> 

CO 

ov 

1 — 1 

o 

1 — 1 

i — 1 

e~- 

c~- 

m 

1X3 

1X3 

1 — 1 

£ 

C\J 

on 

o 

1X3 

03 

CM 

on 

m 

[— 

CM 

o 

1 — 1 

1 - 1 

1 - 1 

CM 

CM 

CM 

CM 

CM 

CM 

on 

> 

r\j 

VO 

oo 

1X3 

o 

O 

o 

O 

0) 

o- 

l n 

o- 

VO 

CM 

CM 

e- 

O 

on 

CO 

03 

•H 

< 

• 

1 - 1 

o 

C\J 

l n 

OD 

1 — 1 

VO 

C— 

OO 

oo 

CO 

o 

LX3 

o 

Cu 

1 — 1 

i — i 

1 - 1 

1 — 1 

1 - 1 

i — i 

1 — 1 

CM 

CM 

0) 

a, 

ii  oo 

m 

< 

on 

\ 

<n>  • 

o 

o 

o 

o 

o 

o 

o 

co 

o 

o 

o 

O 

O 

m 

• 

tD  C\J 

• 

• 

£ 

CO 

> 

LT\ 

o 

o 

LT\ 

■=r 

m 

i — 1 

1X3 

oo 

on 

o 

O 

o 

> 

£ 

on 

VO 

03 

CM 

m 

03 

ov 

03 

o 

CM 

CO 

1 - 1 

i — i 

i — i 

1 - 1 

1 — 1 

1 - 1 

CM 

CM 

CM 

£ 

-p 

CO 

CD 

0- 

o 

CM 

VO 

CM 

o 

o 

o 

o 

o 

O 

o 

O 

£ 

C\i 

E- 

oo 

VO 

LT\ 

on 

oo 

vo 

on 

co 

VO 

-P 

£ 

< 

• 

• 

•H 

3 

O- 

o 

CM 

[— 

o 

i — 1 

CM 

on 

on 

on 

-M" 

vo 

1 - 1 

o 

O 

II  o 

1 — 1 

i — 1 

i — 1 

i — I 

i — i 

i — I 

i — 1 

1 — 1 

CM 

o 

in 

Cl) 

•H 

=0  • 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

> 

& 

U 

1 

a 

> 

on 

VO 

VO 

l n 

1X3 

ov 

03 

c— 

ITS 

CM 

Tf 

cd 

6 

CM 

in 

oo 

1 - 1 

CM 

on 

-=r 

in 

CO 

on 

•H 

£ 

1 - 1 

i — 1 

i — i 

i — l 

i — I 

i — I 

i — i 

i — 1 

CM 

£ 

M 

(J1 

o 

•H 

£ 

r- 

-^r 

CM 

LO 

i — 1 

CO 

i — I 

o 

i — i 

0- 

CM 

j 

cd 

on 

03 

o 

f— 

1X3 

o 

o 

i — 1 

t>- 

i — i 

< 

• 

ii 

o 

II  oo 

o 

1 - 1 

-^r 

LT\ 

vo 

VO 

vo 

vo 

vo 

03 

i 

P-I 

CO  O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

i — 1 

> 

1 — 1 

VO 

O 

in 

o 

1X3 

03 

vo 

m 

i 

£ 

CM 

1X3 

vo 

vo 

t- 

VO 

vo 

vo 

CO 

o 

<u 

1 — 1 

1 — 1 

rO 

TO  O  ^ 

cd 

(u  hOra 

EH 

•H  cd  -P 

CM 

on 

l n 

vo 

13- 

oo 

03 

o 

1 — 1 

CM 

CO 

• — 1  -P  i — 1 

O..H  O 

1 - 1 

i — i 

i — 1 

i — i 

i — 1 

1 - 1 

1 — 1 

1 - 1 

CM 

CM 

CM 

CM 

CM 

U  O  > 

<C  > 

Different  liquid  velocities. 


C-l6 
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(Vrt)  q.U0j;j;no 


Applied  Voltage  (volts) 

Figure  40:  Curves  of  Polarographic  Current  vs.  Applied  Voltage  for  Probe 


Current  (yA) 
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Figure  4l:  Curves  of  Polarographic  Current  vs.  Applied 
Voltage  for  Probe  2 
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Calibration  of  the  Platinum  Probes. 

The  apparatus  used  for  the  calibration  of  the  plat¬ 
inum  probes  is  shown  in  Figure  11.  The  technique  involved 
in  this  experiment  is  identical  to  that  of  the  measurement 
of  velocity  profiles  in  the  column,  with  the  exception  that 
the  point  velocity  at  the  probe  is  known  and  the  limiting 
current  at  this  velocity  is  measured.  The  data  for  the 
calibration  of  probes  1  and  2  is  given  in  Table  11  and  the 
calibration  curves  are  shown  in  Figures  42  and  43. 

Measurement  of  Point  Velocities  in  the  Column. 

In  order  to  ascertain  the  radial  position  in  the 
column  at  which  the  liquid  velocity  was  measured,  the  probes 
were  calibrated  in  terms  of  a  characteristic  length  which 
was  a  function  of  the  radial  position.  The  details  of  the 
probe  setting  and  the  definition  of  the  characteristic 
length  are  given  in  Figure  44.  The  data  for  the  calibration 
of  the  radial  position  vs.  characteristic  length  is  given 


in  Table  12. 


. 
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Resistance  in  the  calibration  circuit  =  10.8kfi. 
mV  =  Voltage  on  digital  voltmeter  (mV). 
yA  =  Current  through  probe  (yA). 

=  Liquid  velocity  at  the  probe  (cms/sec). 
r  =  Radius  of  moving  arm  (cms). 


r 

Probe  1 

Probe 

2 

mV 

yA 

V 

mV 

yA 

V 

120.0 

11.1 

0.555  I 

I  135.0 

12.5 

1.180 

100.0 

9.25 

0.395 

101.0 

9.35 

0.549 

5.0 

111.0 

10.3 

0.441 

62.0 

5.74 

0.264 

101.0 

9.35 

0.368 

73.0 

6.75 

0.290 

79.6 

7.37 

0.206 

79.0 

7.31 

0.358 

97.0 

8.98 

0.281 

87.0 

8.06 

0 . 460 

114 . 0 

10.60 

0.440  . 

121.0 

11.2 

0.945 

197.0 

18 . 20 

2.01 

2  30 

21.2 

3.17 

172.0 

15.90 

1.46 

355 

32.8 

6.35 

127 . 0 

11.80 

0 .65 

310 

28.6 

5.06 

122.0 

11.30 

0.63 

265 

24 . 5 

3.92 

138.0 

12 . 80 

0 . 82 

204 

18.9 

2.50 

166.0 

15.40 

1.22 

180 

16.6 

1.79 

201.0 

18.60 

1.95 

159 

14.7 

1.42 

L0 .0 

220.0 

20 .40 

2.50 

146 

13.5 

1.30 

246 . 0 

22 . 80 

3.16 

136.0 

12.60 

0 . 856 

153.0 

14 . 20 

1.11 

322.0 

29 . 80 

5 . 10 

359.0 

33.20 

5 . 86 

299.0 

27.70 

4.39 

280 . 0 

25.90 

4.10 

269.0 

0 

OJ 

3.76 

176 . 0 

16 . 30 

1.52 

149 . 0 

13.80 

0 . 966 

Table  11  -  Calibration  Data  for  Probes  1  and  2. 
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Figure  42:  Calibration  Curve  for  Probe 
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Figure  43:  Calibration  for  Probe 
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Y#  =  Characteristic  Length 
for  Probe  2 


Figure  Details  of  the  position  of  the  probe 

with  respect  to  the  column  and  details 
of  the  characteristic  length 


. 
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With  reference  to  Figure 
Length  of  Probe  1  = 
Length  of  Probe  2  = 
Then  X* 

* 

and  Y  = 


44: 

196.3  mm . 

197 .3  mm. 

196. 3- 23. 0-3. 3-101. 0+r1 

69 . 0+r . 

197. 3- 23. 5-2. 8-101+r1 


=  70+rx 


During  the  process  of  experimentation  the  probe 
was  shortened  by  3*3  mm.  to  restore  a  broken  electrical 
contact;  table  12  thus  gives  both  the  original  and  the 
later  calibrations  of  probe  1. 
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rl/r0 

X* (mm) 
Prior 
to  Short¬ 
ening 

X* (mm) 
After 
Short¬ 
ening 

Y* (mm) 

0.0 

69.0 

65.7 

70.0 

0.1 

79.1 

75.8 

80.1 

0.2 

89.2 

85.9 

90.2 

0.3 

99.3 

96.0 

100.3 

0.4 

109.4 

106.1 

110.4 

0.5 

119.5 

116.2 

120.5 

0 . 6 

129.6 

126 . 3 

130.6 

0.7 

139.7 

136.4 

140.7  ; 

0 . 8 

149.8 

146.5 

150.8 

0.9 

159.9 

156.6 

160.9 

-0.1 

58.9 

55.6 

59.9 

-0.2 

48.8 

45.5 

49.8 

-0.3 

38.7 

35.4 

39.7 

• 

o 

1 

28.6 

25.3 

29.6 

-0.5 

18.5 

15.2 

19.5 

-0.6 

8.4 

5.1 

9.4 

Table  12  -  Radial  position  vs. 

Characteristic  length 
for  probes  1  &  2. 
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The  negative  values  of  r^/r^  in  table  12  indicate 
the  position  past  the  centre  of  the  column. 

At  points  very  close  to  the  orifice  and  the  column 
walls  the  radial  position  is  reported  in  terms  of  the 
characteristic  length  only. 

The  data  for  measurement  of  point  velocities  cor¬ 
responding  to  Figure  16,  17  and  18  is  given  in  Tables  13, 
14  and  15  respectively. 
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Table  13:  Point  Velocity  Measurements  In  different  diametric  directions 


C-27 


in 

i — i 

P 

G 

3 

-=r 

i — i 

i — i 

i — i 

-=r 

in 

CQ 

0 

i — i 

i — i 

P 

0 

i — 1 

i — 1 

0 

0 

o 

p 

i — i 

i — i 

> 

cd 

P 

P 

Eh 

cd 

cd 

Ln 

Eh 

Eh 

CO 

G 

i _ i 

i _ i 

• 

O 

1 — 1 

tin 

G 

O 

O 

o 

II 

0 

•H 

0 

0 

g 

P 

0 

0 

p 

3 

•H 

\ 

\ 

G 

p 

CQ 

0 

0 

0 

3 

0 

O 

a 

a 

G 

•H 

1 — 1 

a 

o 

o 

G 

CQ 

a 

3 

a 

G 

o 

o 

0 

•H 

o 

o 

oo 

S 

G 

• 

• 

faO 

• 

o 

P 

1 — 1 

OJ 

G 

-=r 

o 

p 

0 

•H 

i — l 

o 

a 

II 

11 

P 

g 

o 

p 

cd 

•H 

g 

II 

• 

G 

•H 

>5 

a 

in 

0 

P 

p 

•H 

i — i 

0 

CM 

' — ' 

•H 

i — 1 

• 

i — 1 

0 

o 

CM 

N 

II 

g 

on 

o 

G 

CM 

o 

Q 

1 — 1 

O 

II 

O 

0 

•H 

0 

<P 

G 

G 

P 

G 

> 

G 

3 

■H 

O 

0 

0 

0 

P 

P 

•H 

P 

hO 

P 

P 

3 

G 

P 

•H 

cd 

0 

P 

G 

O 

•H 

3 

P 

a 

0 

O 

0 

CT1 

i — l 

cd 

P 

a 

O 

•H 

o 

•H 

3 

a 

fa£) 

a 

PI 

> 

P 

0 

G 

0 

EH 

•H 

P 

0 

p 

0 

G 

P 

cd 

50 

0 

O 

3 

a 

0 

cd 

•H 

•rH 

CQ 

0 

G 

cd 

G 

i — 1 

a 

CQ 

P 

0 

p 

0 

a 

0 

0 

a 

cd 

> 

a 

G 

G 

ol 

Q 

< 

< 

O 

a 

CO 

n 

a 

i — 1 

30 

OO 

1 — 1 

• 

n 

CM 

MO 

-=r 

r— 

• 

o 

-=r 

i — 1 

o 

n 

m 

a 

CM 

o 

cn 

• 

on 

n 

• 

1 — 1 

o 

n 

1 — 1 

on 

i — 1 

O 

on 

VO 

w 

a 

CM 

crv 

o 

p 

• 

1 — 1 

a 

p- 

• 

< 

o 

a 

G 

0 

X 

CM 

CM 

1 

o 

o 

3 

a 

P 

CQ 

G 

3 

G 

i — 1 

1 - 1 

1 

o 

o 

0 

a 

P 

P 

a 

o 

CQ 

n 

1 - 1 

o 

o 

VO 

G 

a 

on 

• 

0 

n 

0- 

P 

a 

• 

B 

h 

o 

G 

G 

p- 

1 - 1 

o 

on 

ov 

O 

a 

CM 

• 

co 

•rH 

-=r 

on 

p 

1 — 1 

• 

ctf 

r  \ 

o 

•H 

a 

•H 

on 

1 - 1 

O 

a 

CO 

p 

«- 

a 

CM 

• 

o 

G 

i — 1 

n 

1 — 1 

0 

o 

P 

a 

• 

H 

p 

f") 

a 

II 

< 

a 

CM 

CM 

1 

o 

o 

in 

a 

a 

c\ 

-=3- 

a 

#\ 

i — 1 

1 - 1 

1 

o 

o 

on 

a 

P-H 

7 

*\ 

o 

CM 

• 

o 

0 

a 

0 

W  h0 

•  0 

0 

P 

CO  *H 

a  0 

p  \ 

*\ 

O 

0  0  0 

0  \ 

0  o  o 

i — 1 

G 

3  g  a 

rH 

3  i — l  o 

a 

a 

cd  a  w 

►a 

2 . 

l/) 

a> 

/* 

w 

- -7“  /  V" 

{  /  )j,  oc  0 in  -  -&)  _'J  ft  -  0  ' 


C-27A 


-=r 

VO 

on 

in 

in 

on 

in 

X 

on 

crv 

-=r 

1 — 1 

i — 1 

o 

o 

o 

o 

O 

1 — 1 

on 

ov 

co 

• 

• 

o 

1 — 1 

I — 1 

i — 1 

i — 1 

i — 1 

rH 

rH 

i — 1 

rH 

rH 

VO 

o 

CTV 

on 

on 

OJ 

on 

-=d" 

CTV 

i — 1 

on 

in 

cd 

> 

CT\ 

in 

-=r 

-=r 

vo 

CX\ 

s 

l — 1 

i — 1 

i — I 

i — i 

i — i 

i — 1 

i — 1 

1 — 1 

l — 1 

l — 1 

* 

in 

o 

>H 

• 

• 

in 

in 

* 

V0 

0- 

X 

1 - 1 

‘P 

LO 

CVJ 

on 

in 

on 

i — i 

on 

on 

CT\ 

vo 

in 

VO 

OJ 

1 — 1 

o 

o 

o 

o 

o 

o 

1 — 1 

-=r 

• 

• 

• 

O 

1 — 1 

i — 1 

I — 1 

1 — 1 

1 — 1 

l — 1 

1 — 1 

1 — 1 

1 — 1 

rH 

> 

LT\ 

CTv 

on 

CVJ 

1 — 1 

CM 

OJ 

vo 

o 

0- 

1 - 1 

m 

-^r 

-=j- 

-=r 

-=r 

-=r 

-=r 

in 

vo 

cd 

e 

1 - 1 

i — 1 

i — i 

i — i 

rH 

• — I 

rH 

i — l 

i — i 

1 — 1 

rH 

* 

LT\ 

CM 

>H 

• 

c\ 

in 

in 

* 

VO 

IX- 

X 

1 — 1 

<=^ 

CO 

CO 

VO 

cr\ 

ov 

in 

on 

in 

CO 

CO 

0- 

LT\ 

1 — 1 

1 — 1 

o 

o 

o 

o 

o 

o 

o 

o 

Co 

• 

• 

• 

o 

1 — 1 

1 — 1 

1 — 1 

1 — 1 

1 — 1 

rH 

I — 1 

l — 1 

1 — 1 

l — 1 

> 

ov 

OJ 

o 

VO 

vo 

on 

CM 

on 

in 

in 

-^r 

o 

Ln 

in 

-^r 

-=T 

on 

cd 

e 

1 — 1 

l — 1 

1 — 1 

i — i 

rH 

i — l 

rH 

i — l 

rH 

i — l 

i — l 

* 

LT\ 

>H 

• 

• 

in 

on 

* 

vo 

IX 

X 

1 — 1 

oo 

'N1 

oo 

on 

in 

in 

-=r 

o 

cr\ 

ov 

o 

IX 

OO 

CO 

on 

OJ 

l — 1 

rH 

1 — 1 

o 

o 

i — 1 

o 

X 

• 

• 

• 

• 

o 

1 — 1 

1 — 1 

1 — 1 

i — 1 

rH 

rH 

1 — 1 

i — 1 

1 — 1 

o 

> 

o\ 

VO 

OJ 

VO 

in 

on 

CM 

OJ 

on 

1 — 1 

on 

1 — 1 

-^r 

on 

on 

on 

on 

on 

on 

on 

C\J 

cd 

s 

1 — 1 

i — i 

i — l 

1 — 1 

1 — 1 

l — 1 

i — 1 

1 — 1 

1 — 1 

rH 

1 — 1 

* 

LPi 

on 

X 

• 

in 

on 

* 

VO 

X 

X 

1 — 1 

OJ 

-=T 

X 

CO 

CO 

in 

i — 1 

on 

on 

rH 

on 

<p 

CO 

OJ 

1 — 1 

o 

o 

o 

o 

o 

o 

o 

vo 

D 

• 

• 

o 

1 — 1 

1 — 1 

r— 1 

rH 

rH 

rH 

rH 

l — 1 

i — 1 

o 

> 

X 

in 

1 — 1 

in 

in 

on 

I — 1 

OJ 

CM 

i — 1 

on 

CM 

in 

in 

-^r 

-^r 

-=d~ 

-=r 

-=r 

l — 1 

S 

1 — 1 

l — 1 

I — 1 

i — I 

l — 1 

i — i 

i — 1 

rH 

i — i 

rH 

cd 

in 

on 

X 

• 

in 

in 

k 

VO 

vo 

X 

1 - 1 

o 

p 

CTV 

jl 

\ 

• 

CO 

x 

vo 

in 

on 

CM 

i — 1 

( r 

i — 1 

o 

• 

• 

• 

• 

• 

• 

• 

• 

Jr 

u 

o 

o 

o 

o 

o 

o 

O 

o 

o 

> 

0 

£ 

co 

\ 

CO 

Sh 

6 

0 

o 

P 

w 

0 

6 

0 

-p 

o 

i — i 

o 

o 

u 

> 

cx 

i — i 

p 

cd 

cd 

-P 

•H 

X 

hO 

p 

•H 

•H 

Td 

O 

o 

C 

1 — 1 

o 

0 

> 

0 

faO 

cd 

•H 

p 

3 

rH 

O1 

o 

•H 

> 

ii  ii 


X*,Y*  =  Characteristic  length  (mm) 


C-27B 


o 


on 

on 

p 


in 

No 


CM 


on 

p 


w 

o 

o 

on 

on 

H 

o 

on 

on 

i — i 

CM 

CM 

\ 

CA 

3D 

CM 

o 

o 

o 

o 

CD 

o 

i — l 

rH 

in 

s 

* 

• 

• 

o 

1 - 1 

1 - 1 

1 — 1 

l — 1 

1 — 1 

1 — 1 

1 - 1 

i — 1 

i — l 

i — 1 

o 

on 

CM 

CM 

1 — 1 

o 

CM 

CM 

t — 

co 

in 

o 

30 

in 

-^r 

-=r 

-=r 

-=r 

o 

1 - 1 

i — 1 

i — l 

i — l 

i — l 

i — i 

i — l 

1 — 1 

i — I 

1 — 1 

1 — 1 

i — 1 

n 

on 

on 

• 

• 

>a 

on 

P 

VO 

•H 

l — 1 

O 

O 

on 

p 

co 

o 

o 

o 

on 

p 

on 

co 

in 

1 - 1 

d) 

LO 

o 

03 

o 

o 

o 

o 

o 

i — 1 

1 — 1 

0- 

> 

I — 1 

1 — 1 

O 

1 — 1 

1 — 1 

1 — 1 

i — 1 

1 — 1 

i — 1 

1 — 1 

o 

i — 1 

-=r 

00 

o 

o 

o 

CM 

03 

CM 

CM 

•H 

3 

p 

-=r 

on 

■=r 

■=r 

-=r 

-=r 

un 

CM 

o1 

l — 1 

1 — 1 

l — 1 

1 — 1 

1 — 1 

i — i 

l — 1 

i — I 

1 - 1 

1 — 1 

1 - 1 

•H 

on 

0 

• 

• 

hO 

on 

on 

ctf 

0- 

VO 

£ 

1 — 1 

0 

> 

CM 

CO 

1 — 1 

VO 

oo 

o 

on 

LT\ 

on 

o 

o 

CO 

1 - 1 

03 

03 

03 

03 

o 

o 

O 

l — 1 

CM 

co 

£ 

* 

• 

• 

cd 

1 - 1 

o 

O 

o 

o 

1 — 1 

l — 1 

rH 

rH 

1 - 1 

o 

30 

CO 

0- 

oo 

o 

CM 

on 

03 

on 

VO 

M 

o 

on 

on 

on 

on 

-=r 

-^r 

-=r 

in 

CM 

p 

i — l 

1 — 1 

l — 1 

1 — 1 

I — 1 

pH 

i — I 

i — i 

i — I 

l — 1 

rH 

p 

on 

• 

on 

co 


o 

CM 

in 

CO 

m 

o 

-=r 

in 

0 

D — 

o 

03 

o 

i — 1 

i — i 

i — 1 

S 

• 

• 

• 

• 

• 

• 

• 

O 

i — 1 

o 

l — 1 

rH 

i — l 

1 — 1 

>3 

30 

o 

in 

o 

on 

in 

VO 

4-> 

O 

on 

CM 

on 

on 

on 

on 

o 

1 - 1 

1 — 1 

1 - 1 

l — 1 

i — 1 

I — 1 

i — 1 

o 

03 

o 

in 

VO 

03 

o 

in 

in 

VO 

CM 

03 

IP 

03 

co 

03 

03 

o 

o 

o 

1 — 1 

CM 

P 

o 

o 

o 

o 

O 

1 — 1 

l — 1 

rH 

1 — 1 

1 - 1 

o 

in 

CO 

o 

p 

03 

o 

on 

on 

o 

in 

CM 

on 

on 

on 

on 

-=T 

in 

in 

CM 

i— H 

rH 

1 — 1 

rH 

l — 1 

i — i 

rH 

rH 

rH 

( — 1 

1 — 1 

on 

• 

on 

co 


r 

(' 

i — 1 

CM 

on 

in 

vo 

p 

CO 

03 

\  / 

• 

• 

• 

oL 

o 

o 

o 

o 

o 

o 

o 

o 

o 

0 

> 

p 

£ 

•H 

O 

Oh 


0 
i — I 
P 
cti 
Eh 


C-28 


oo 

OJ 

co 

o 

on 

on 

vo 

1 — 1 

on 

rH 

1 — 1 

P 

o 

IP 

on 

OJ 

rH 

I — 1 

o 

o 

rH 

CVJ 

co 

P 

• 

OJ 

OJ 

Ol 

CVJ 

OJ 

OJ 

CVJ 

OJ 

CVJ 

OJ 

CVJ 

p 

in 

rH 

on 

ov 

CO 

vo 

o 

on 

CO 

> 

co 

OJ 

rH 

o 

ov 

CO 

CO 

CO 

o 

o 

CVJ 

in 

PC 

s 

1 — 1 

OJ 

OJ 

OJ 

1 — 1 

1 — 1 

1 — 1 

1 — 1 

CVJ 

CVJ 

CVJ 

st: 

in 

CO 

P 

• 

• 

in 

sK 

VO 

p 

X 

rH 

on 

vo 

i — 1 

on 

o 

on 

on 

p 

on 

o 

vo 

p 

o 

c — 

-=r 

OJ 

Ol 

rH 

l — 1 

CO 

i — 1 

OJ 

D 

• 

OJ 

OJ 

CVJ 

CVJ 

CVJ 

Ol 

CVJ 

rH 

OJ 

CVJ 

CVJ 

in 

p 

CVJ 

-=r 

on 

o 

o 

CVJ 

CO 

on 

-^r 

> 

co 

CVJ 

1 — 1 

o 

o 

o 

o 

CO 

CO 

o 

1 — 1 

CP 

S 

l — 1 

OJ 

OJ 

CVJ 

CVJ 

CVJ 

CVJ 

1 — 1 

1 — 1 

Ol 

CVJ 

st: 

in 

p 

P 

• 

• 

in 

on 

sK 

VO 

p 

X 

l — 1 

o 

p 

CO 

in 

o 

on 

o 

CO 

on 

co 

on 

P 

1 — 1 

vo 

on 

OJ 

CVJ 

1 — 1 

1 — 1 

CO 

i — 1 

o 

1 — 1 

P 

• 

• 

OJ 

OJ 

OJ 

CVJ 

CVJ 

OJ 

OJ 

1 — 1 

OJ 

CVJ 

Ol 

CO 

on 

1 — 1 

m 

on 

o 

OO 

on 

CO 

p 

o 

> 

CO 

CVJ 

rH 

o 

o 

o 

CO 

CO 

CO 

CO 

o 

on 

S 

1 — 1 

CVJ 

OJ 

CVJ 

OJ 

CVJ 

1 — 1 

rH 

1 — 1 

rH 

CVJ 

SK 

m 

p 

>H 

• 

• 

*\ 

in 

on 

sK 

VO 

p 

X 

rH 

o 

p 

p 

p 

o 

on 

■p 

o 

p 

oo 

OJ 

> 

0 

P 

OJ 

CO 

in 

-=d" 

on 

CVJ 

OJ 

OJ 

1 — 1 

o 

6 

w 

D 

• 

— ' 

\ 

OJ 

OJ 

CVJ 

CVJ 

CVJ 

OJ 

CVJ 

OJ 

CVJ 

CVJ 

CVJ 

C/3 

^ — S 

U 

S 

S 

C\J 

o 

vo 

on 

VO 

p 

on 

o 

p 

o 

vo 

co 

0 

o 

s 

PC 

> 

CO 

1 — 1 

o 

CO 

CO 

co 

CO 

co 

CO 

co 

p 

-P 

v'— ' 

S 

1 — 1 

OJ 

OJ 

1 — i 

1 — 1 

rH 

1 — 1 

1 — 1 

1 — 1 

1 — 1 

l — 1 

0 

s 

0 

p 

-p 

P 

p 

SK 

in 

CO 

1 — 1 

O 

hO 

P 

• 

o 

Sh 

c 

in 

on 

> 

a 

0 

SK 

vo 

p 

1 — 1 

X 

i — 1 

aj 

cd 

o 

P 

•H 

p 

CO 

in 

o 

OO 

co 

p 

OJ 

p 

o 

in 

•H 

>> 

p 

p 

1 — 1 

p 

CVJ 

rH 

CVJ 

CVJ 

1 — 1 

1 — 1 

1 — 1 

1 — 1 

co 

hO 

P 

w 

p 

• 

•rH 

•H 

•H 

OJ 

OJ 

OJ 

OJ 

OJ 

OJ 

CVJ 

CVJ 

OJ 

CVJ 

1 — 1 

O 

o 

0 

rH 

CO 

in 

co 

vo 

vo 

1 — 1 

CO 

1 — 1 

co 

vo 

C 

r~H 

P 

i — 1 

> 

o 

OJ 

o 

CO 

o 

o 

o 

CO 

o 

CO 

oo 

o 

0 

O 

PC 

s 

OJ 

OJ 

CVJ 

l — 1 

CVJ 

CVJ 

OJ 

1 — 1 

CVJ 

1 — 1 

1 — 1 

> 

0 

£-1 

hO 

T3 

* 

in 

-=3" 

oj 

•rH 

p 

P 

• 

• 

P 

3 

o 

•V 

in 

o 

f - 1 

o1 

SK 

vo 

vo 

o 

•rH 

II 

X 

1 — 1 

1 — 1 

SK 

o 

it 

n 

I>H 

£h 

r 

\ 

CO 

CO 

p 

vo 

in 

-^r 

on 

CVJ 

1 — 1 

/f 

> 

p 

* 

i — 1 

CL 

e 

p 

X 

U 

o 

o 

o 

o 

o 

o 

o 

o 

o 

£ 

C-28A 


on 

• 

CT\ 

VO 


0- 

CO 

vo 


vo 

vo 


i>- 

co 

VO 


vo 

CO 

oo 

t— 

on 

o 

on 

on 

CM 

OO 

vo 

1 - 1 

o 

ov 

o 

i — 1 

CM 

CM 

on 

in 

on 

CM 

CM 

1 — 1 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

C\J 

CM 

C— 

CM 

m 

CO 

CO 

VO 

CX\ 

vo 

O 

ov 

Cn 

o\ 

cr\ 

o 

o 

O 

CM 

1 — 1 

CM 

CM 

l — 1 

1 - 1 

1 — 1 

1 — 1 

CM 

CM 

CM 

CM 

CM 

0— 


o 

vo 


OO 

VO 

on 

un 

ov 

CO 

CO 

1 - 1 

CO 

CM 

CO 

-=T 

o 

ov 

o\ 

ov 

o 

1 — 1 

CM 

CM 

t"- 

CM 

CM 

I — 1 

l — 1 

I — 1 

CM 

CM 

CM 

CM 

CM 

CM 

m 

VO 

o 

i — 1 

on 

0- 

CM 

on 

VO 

m 

in 

1 — 1 

CTV 

CTV 

cn 

OV 

ov 

o 

o 

o 

CM 

1 — 1 

CM 

1 - 1 

1 - 1 

l — 1 

i — 1 

1 — 1 

CM 

CM 

CM 

CM 

CM 

D — 

• 

o 

vo 


COJ 

CTj 

u 

-=r 

on 

on 

t— 

vo 

CO 

o 

CM 

0- 

crv 

O 

(D 

on 

o 

CTV 

ov 

o 

1 — 1 

on 

on 

on 

0- 

VO 

> 

• 

• 

• 

cd 

CM 

CM 

1 - 1 

! - 1 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

ov 

m 

O 

CM 

VO 

CM 

0- 

CO 

o 

CO 

O 

cti 

o 

OA 

ov 

OV 

CT\ 

O 

o 

o 

1 — 1 

CM 

CM 

CM 

1 — 1 

1 — 1 

1 - 1 

1 - 1 

CM 

CM 

CM 

CM 

CM 

CM 

G 

s'  % 

o 

o 

Cm 

CD 

0- 

0— 

CO 

• 

• 

-p 

\ 

ov 

o 

£ 

£ 

vo 

vo 

CD 

o 

e 

— 

CD 

u 

o 

10- 

oo 

o 

oo 

oo 

VO 

CO 

3 

• 

m 

1 — 1 

o 

o 

1 — 1 

on 

on 

co 

CM 

• 

• 

• 

• 

• 

• 

• 

cd 

l — 1 

CM 

CM 

CM 

CM 

CM 

CM 

CD 

II 

oo 

vo 

OO 

CM 

VO 

on 

LO 

>> 

LO 

CO 

oo 

OO 

ov 

0"\ 

>s 

P 

1 - 1 

1 — 1 

1 — 1 

! - 1 

1 - 1 

l — 1 

1 - 1 

•p 

•H 

•H 

O 

i — 1 

ov 

on 

l — 1 

m 

o 

O 

CM 

1 — 1 

VO 

o 

•H 

ov 

ov 

ov 

o 

o 

CM 

on 

on 

OO 

to- 

o 

• 

PL. 

1 — 1 

1 — 1 

i — 1 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

ov 

on 

o 

-=r 

vo 

on 

to- 

CO 

CM 

1 - 1 

•  . 

CO 

cn 

crv 

ov 

crv 

o 

o 

o 

1 - 1 

on 

CM 

in 

1 — 1 

l — 1 

1 — 1 

1 — 1 

l — 1 

CM 

CM 

CM 

CM 

CM 

CM 

i — 1 

o 

vo 


i — 1 

CM 

on 

^r 

m 

vo 

t— 

CO 

CTV 

• 

.» 

• 

• 

» 

• 

• 

• 

• 

o 

O 

o 

o 

o 

o 

o 

o 

o 

<u 


o 

o 

I — I 
CD 
> 

P 

C 


CD 

I - 1 

P 

cd 

Eh 


CD 

> 

•H 

3 

O' 

•H 


00 

’  . 

• 

ru 

. 

t=.. 

r 

D-l 


Appendix  D 

Evaluation  of  The  Mass  Transfer  Coefficient  from  CO^-Water  Data 

The  amount  of  mass  transfer  during  bubble  formation 
is  obtained  by  the  difference  (V-V*)  where 

V  =  The  equivalent  volume  of  a  gas  bubble  determined 
by  the  equation  G  =  vV. 

V*  =  The  volume  of  the  gas  bubble  at  detachment  as 
determined  from  photographs. 

The  general  equation  for  total  mass  transfer  upto 
a  time  T^,  is  given  by 

Tf 

NA  =  iCt  -  CA>  kLav  /  A(t)  dt  - (D-X) 

o 

In  order  to  determine  k^av  the  following  data  must  be  known: 

-  The  amount  of  mass  transfer  during  bubble  formation. 
This  is  given  by  (V-V*). 

t. 

-  The  concentration  of  the  solute  CO^  in  the  bulk 
liquid  phase . 

C*  -  The  concentration  of  the  solute  CO^  at  the  inter¬ 
face.  If  the  resistance  to  mass  transfer  at  the 
interface  is  assumed  negligible  then  is  given 
by  the  solubility  of  the  gas  in  the  liquid  phase 
at  the  system  temperature  and  gas  partial  pressure. 
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J  A(t)  dt  =  The  value  of  the  area  integral  which  represents 
0 

the  variation  of  the  bubble  surface  area  with 
time  during  bubble  formation. 

The  value  of  is  essentially  zero.  This  can  be  shown  by 
the  following  sample  calculation. 


Size  of  water  sample 
Titre  value  of  a  blank 
titration 

Titre  value  for  a  CO^ 

titration 
Normality  of  HC1 


=  20  cc . 


=  21.2cc  HC1 


=  20.95  cc  HC1 
=  0.0067N 


=  (21.2-20.95)  x  0.0067 

moles 

20 

litre 

=  0.25  x  0.0067  =  8.38  X 

10  moles 

20 

litre 

=  8.38  X 

n  n-8  n 

10  moles 

cc  . 

3.59  x  10 


gm  moles 
cc  soln 


at  the  system  conditions 


.*o  (C*  -  c£)  =  (3.59-0.008)  x  10~5  =  3.582  x  10  5  moles 


The  error  involved  in  assuming  =  0  is  -  0.2%. 

This  is  well  within  the  limitation  of  the  experiment. 
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The  value  of  CA  Is  evaluated  from  solubility  data 


as  follows : 


Temperature  of  investigation  =  25°c 


Assuming  Henry’s  Law  to  hold  for  CO^-water 
The  Henry’s  Constant  =  1.64  x  10^  [73] 


H 


where 

=  moles  CO^ 

moles  Solution 

PA  =  partial  pressure  of  CO^ 

The  gas  CO^  enters  the  orifice  saturated  with  water  vapour. 

The  vapour  pressure  of  water  at  25°c  is  given  by  23.76  mm 

Hg  [74]  =  0.46psia. 

Total  system  pressure  =  14.8  psia 

.’.  pA  =  14.8  -  0.46  =  14.34  psia 

Atmospheric  pressure  =  13.5  psia 

.  '  .  pn  =  14.34  =  1.062  atm. 

13.5 


=  0.6477  x  10-^  gm  moles 

mole  soln. 


XA  =  1.062 

1.64  x  10 


_5 

3.59  x  10  gm  moles 
cc  soln. 


=  0.6477  x  10 
18.04 


The  area  integral  is  evaluated  from  Figures  23,  24  and  25. 


frv 


. 


